Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 432,414 G→T 3.7% intergenic (+17/+131) RCAP_rcc00403 → / ← RCAP_rcc00404 hypothetical protein/Sua5/YciO/YrdC/YwlC family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034432,4140G→T 3.7% 8.1 ‑4.2 ‑0.3162intergenic (+17/+131)RCAP_rcc00403/RCAP_rcc00404hypothetical protein/Sua5/YciO/YrdC/YwlC family protein
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (48/108):  total (53/109)
Fisher's exact test strand distribution p-value = 1.47e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.68e-01

Only 100 of 179 total aligned reads displayed.

GATCACGCGACGAAATGGACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.T.CTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCGGGGGTGAATTTTG  >  NC_014034/432314‑432516
                                                                                                    |                                                                                                        
gATCACGCGACGAAATGGACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGgcg                                                                                                          <  2:2927167/101‑1 (MQ=255)
gATCACGCGACGAAATGGACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGgcg                                                                                                          <  1:1813557/101‑1 (MQ=255)
  tCACGCGACGAAATGGACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCt                                                                                                        <  2:3143095/101‑1 (MQ=255)
      gcgACGAAATGGACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGccc                                                                                                    <  2:957450/101‑1 (MQ=255)
       cgacgaAATGGACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGgccc                                                                                                   >  2:2284049/1‑101 (MQ=255)
          cgaAATGGACCCGGGGGCTGGGCACGATGGGGATCGTTCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCg                                                                                                <  1:3612676/101‑1 (MQ=255)
           gaAATGGACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCGGCCCCCCgc                                                                                               >  2:3065124/1‑100 (MQ=255)
               tGGACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.T.cctg                                                                                         <  2:1569581/101‑4 (MQ=255)
                ggACCCGGGCGCTGGGCACGATGGGGATCGGTCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.T.cctgc                                                                                        >  1:3966807/1‑97 (MQ=255)
                ggACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.T.cctgc                                                                                        <  1:624619/101‑5 (MQ=255)
                 gACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.T.cctgcg                                                                                       <  1:4369590/101‑6 (MQ=255)
                    ccGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGAc                                                                                     <  1:3289210/101‑1 (MQ=255)
                     gcagcgcTGGGCACGATGGGGATAGATCCGCGCACGCTTTCGGCGGCGGCGGGCCAAGACTCAGGCCGGGCAGGAGAGAGCTGACCCCCGCC.CTGCGGACt                                                                                    >  2:726375/4‑101 (MQ=255)
                      gggCGGTGGGCACGATGGGGATGGATCCGCGCGCGCTTGCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTc                                                                                   <  2:1454362/101‑1 (MQ=37)
                      gggCGCGGGGCACGAGGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCTGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTc                                                                                   <  2:10693/101‑1 (MQ=37)
                       ggCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTcc                                                                                  <  1:3893579/101‑1 (MQ=255)
                         cgcTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTcccc                                                                                <  1:2240108/101‑1 (MQ=255)
                          gatGGGCACGATGGGGAGCGATCCGCGCACGCTTGCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTccccc                                                                               <  2:3861814/99‑1 (MQ=37)
                            tGGGCACGATGGGGATCGATCCGCGCCCGCTTTCGGCGGGGGCGGGCCACGCCTGAGGCCGGGCGGGGGGGGGCGGCCCCCCG.TCCTGCGGACTCCCCCCg                                                                             >  2:3587263/1‑101 (MQ=37)
                            tGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCg                                                                             >  1:3674532/1‑101 (MQ=255)
                            tGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCg                                                                             <  1:4064384/101‑1 (MQ=255)
                            tGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCg                                                                             >  2:3640517/1‑101 (MQ=255)
                            tGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCg                                                                             >  2:2162566/1‑101 (MQ=255)
                            tGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCg                                                                             <  1:1265373/101‑1 (MQ=255)
                             gggCACGATGGGGATCGATCCGCGCACGCTTTCGGAGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCgg                                                                            <  1:2234599/101‑1 (MQ=255)
                               gCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGa                                                                          >  2:3309019/1‑101 (MQ=255)
                               gCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGa                                                                          <  1:141001/101‑1 (MQ=255)
                                cACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGat                                                                         <  2:3438372/101‑1 (MQ=255)
                                 aCGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGata                                                                        <  1:558876/101‑1 (MQ=255)
                                 aCGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCCGGgta                                                                        >  2:793445/1‑101 (MQ=37)
                                  cGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGatat                                                                       <  1:1729243/101‑1 (MQ=255)
                                  cGAGGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGatat                                                                       <  1:1642365/101‑1 (MQ=255)
                                     tGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTg                                                                    <  1:262315/101‑1 (MQ=255)
                                     tGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTg                                                                    <  2:2005975/101‑1 (MQ=255)
                                      ggggATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCCGGGTAtttca                                                                   >  1:863520/1‑99 (MQ=25)
                                       gggATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAg                                                                  >  2:478899/1‑101 (MQ=255)
                                         gatcgatcCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGc                                                                <  1:1030510/101‑1 (MQ=255)
                                          atcgatcCGCGCACGCTTTCGGCGGGGGCGGGCCACGCCTTAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGGCTCCCCCCGGGATATTTGAGGCa                                                               <  1:351498/101‑1 (MQ=25)
                                          atcgatcCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCa                                                               <  1:3332201/101‑1 (MQ=255)
                                           tcgatcCGCGCACGCTTTCGGCGGCGGCGGNCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCaa                                                              <  1:1170312/101‑1 (MQ=255)
                                            cgatcCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAg                                                             <  1:1530272/101‑1 (MQ=255)
                                             gatcCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGa                                                            <  1:4144742/101‑1 (MQ=255)
                                             gatcCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGa                                                            >  1:481258/1‑101 (MQ=255)
                                              atcCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGAt                                                           <  2:3698868/101‑1 (MQ=255)
                                               tnNGCNNACGCTTTCGGNGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGTCAAGATg                                                          <  2:1127254/101‑1 (MQ=21)
                                                 cgcgCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGaa                                                        <  2:3168342/101‑1 (MQ=255)
                                                 cgcgCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGACATTTTCGGGCAGATGaa                                                        >  2:857126/1‑101 (MQ=21)
                                                  gcgcACTCTTTCGGCGGCGGGGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGaaa                                                       <  1:1102042/101‑1 (MQ=18)
                                                  gcgcACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGaaa                                                       >  1:2826810/1‑101 (MQ=255)
                                                     cACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGcc                                                    <  1:4299844/101‑1 (MQ=255)
                                                     cACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGcc                                                    <  1:354675/101‑1 (MQ=255)
                                                      aCGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCg                                                   <  1:3722350/101‑1 (MQ=255)
                                                       cGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGa                                                  <  2:4331726/101‑1 (MQ=255)
                                                         cTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATc                                                <  2:2507141/101‑1 (MQ=255)
                                                           ttCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCg                                              <  2:1370921/101‑1 (MQ=255)
                                                           ttCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCg                                              <  1:1968169/101‑1 (MQ=255)
                                                              ggcggcggcggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGc                                           >  1:2776396/1‑101 (MQ=255)
                                                               gcggcggcggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCt                                          <  2:941398/101‑1 (MQ=255)
                                                                cggcggcggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCtt                                         <  2:2928521/101‑1 (MQ=255)
                                                                 ggcggcggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTc                                        <  2:2773287/101‑1 (MQ=255)
                                                                 ggcggcggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTc                                        <  1:1214384/101‑1 (MQ=255)
                                                                 ggcggcggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTc                                        >  2:2595712/1‑101 (MQ=255)
                                                                 ggcggcggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTc                                        <  2:1791387/101‑1 (MQ=255)
                                                                 ggcggcggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCCGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCCTc                                        >  2:3170009/1‑101 (MQ=21)
                                                                    ggcggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGtt                                     <  1:2047272/101‑1 (MQ=255)
                                                                     gggggCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGttt                                    <  1:3136909/99‑1 (MQ=255)
                                                                      cggGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGGTTc                                   >  2:4038378/1‑101 (MQ=255)
                                                                         gCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATc                                <  2:2054554/101‑1 (MQ=255)
                                                                          ccACGCCTGAGGCCGGGGAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCt                               <  1:1267872/101‑1 (MQ=255)
                                                                           cACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCtt                              <  2:173876/101‑1 (MQ=255)
                                                                           cACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCtt                              >  1:2659957/1‑101 (MQ=255)
                                                                           cACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCtt                              >  2:3922964/1‑101 (MQ=255)
                                                                           cACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCCTCGTTTCATcct                              >  2:797111/1‑99 (MQ=255)
                                                                           cACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCAGATCCGGGCTTCGTTTTACCtt                              >  2:59682/1‑101 (MQ=18)
                                                                            aCGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTg                             <  1:1916388/101‑1 (MQ=255)
                                                                            aCGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTg                             <  1:1877662/101‑1 (MQ=255)
                                                                             cGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGc                            <  1:2370166/101‑1 (MQ=255)
                                                                              gCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGcc                           <  1:2575464/101‑1 (MQ=255)
                                                                              gCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGcc                           >  2:4364611/1‑101 (MQ=255)
                                                                                cTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGcccc                         <  1:719163/101‑1 (MQ=255)
                                                                                      ccGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAatat                   <  2:2311039/101‑1 (MQ=255)
                                                                                       cGGGCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATc                  <  2:3175109/101‑1 (MQ=255)
                                                                                        gggCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATcc                 <  1:2715247/101‑1 (MQ=255)
                                                                                         ggacggggtgCTCTGCCCCCCG.G.CTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAA....G....GGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAtt     <  2:1564473/97‑1 (MQ=2)
                                                                                         ggacggggggCTCTGCCCCCCG.G.CTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAA....G....GGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAtt     >  1:1895554/5‑101 (MQ=2)
                                                                                          gacggggggCTCTGCCCCCCG.G.CTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAA....G....GGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttg    >  1:3903003/4‑100 (MQ=2)
                                                                                          gCAGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCg               <  1:1391283/101‑1 (MQ=255)
                                                                                           acggggggCTCTGCCCCCCG.G.CTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAA....G....GGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttgg   >  2:3465243/3‑99 (MQ=2)
                                                                                            cggggggCTCTGCCCCCCG.G.CTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAA....G....GGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggc  <  1:2878129/100‑4 (MQ=2)
                                                                                            aGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCggg             <  2:1915454/101‑1 (MQ=255)
                                                                                            aGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCggg             >  2:1951916/1‑101 (MQ=255)
                                                                                            aGGGGGGCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCggg             >  1:2153747/1‑101 (MQ=255)
                                                                                             ggggggCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCgggg            <  2:458048/101‑1 (MQ=255)
                                                                                             ggggggCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCgggg            >  1:1315167/1‑101 (MQ=255)
                                                                                               ggggCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGCTATTTGAGGCAATACGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCTAAATATCCCGGGGGt          >  2:3000739/1‑101 (MQ=25)
                                                                                               ggggCGCTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCGGGGGt          <  2:1267855/101‑1 (MQ=255)
                                                                                                   cgcTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGGAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATGCCGGGGGCGAAt      >  2:645146/1‑101 (MQ=25)
                                                                                                   cgcTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCGGGGGTGAAt      >  2:127379/1‑101 (MQ=255)
                                                                                                   cgcTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCGGGGGTGAAt      <  2:1397837/101‑1 (MQ=255)
                                                                                                   cgcTGCCCCCCG.TCCTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCGGGGGTGAAt      >  1:337882/1‑101 (MQ=255)
                                                                                                    |                                                                                                        
GATCACGCGACGAAATGGACCCGGGCGCTGGGCACGATGGGGATCGATCCGCGCACGCTTTCGGCGGCGGCGGGCCACGCCTGAGGCCGGGCAGGGGGGCGCTGCCCCCCG.T.CTGCGGACTCCCCCCGGGATATTTGAGGCAAGATGAAAGCCGATCCGGGCTTCGTTTCATCTTGCCCCAAATATCCCGGGGGTGAATTTTG  >  NC_014034/432314‑432516

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG