Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 530,939 +G 2.6% intergenic (+17/+124) valS → / ← RCAP_rcc00497 valyl‑tRNA synthetase/UbiA prenyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034530,9392.→G 2.6% 2.5 0.0 ‑0.1115intergenic (+17/+124)valS/RCAP_rcc00497valyl‑tRNA synthetase/UbiA prenyltransferase
Reads supporting (aligned to +/- strand):  new base (1/2):  ref base (43/69):  total (44/71)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.20e-01

Only 100 of 139 total aligned reads displayed.

GAAGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTC..CCGCGCTTGTCGGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCAA  >  NC_014034/530839‑531022
                                                                                                     |                                                                                    
gAAGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTc                                                                                       >  2:2575407/1‑101 (MQ=255)
 aaGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTct                                                                                      <  2:2232954/101‑2 (MQ=255)
 aaGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTct                                                                                      <  2:2015408/101‑2 (MQ=255)
 aaGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTct                                                                                      <  2:1718474/101‑2 (MQ=255)
  aGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctg                                                                                     <  1:103595/101‑3 (MQ=255)
   gCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgc                                                                                    <  2:1785140/101‑4 (MQ=255)
   gCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctcc                                                                                    >  1:381454/1‑98 (MQ=255)
      cgcgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgcccc                                                                                 >  1:1126552/1‑95 (MQ=255)
       gcgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccc                                                                                >  2:3692051/1‑94 (MQ=255)
       gcgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccc                                                                                <  2:2199773/101‑8 (MQ=255)
       gcgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccc                                                                                <  1:3586054/101‑8 (MQ=255)
        cgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgcccccc                                                                               <  2:3176491/101‑9 (MQ=255)
         gAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccccg                                                                              >  2:1601768/1‑92 (MQ=255)
         gAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccccg                                                                              <  1:753047/101‑10 (MQ=255)
          aaaaCCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccccgc                                                                             >  1:2344928/1‑91 (MQ=255)
                               aaGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATAttt                                                       >  1:1974813/1‑101 (MQ=25)
                                aggggcggcGAAGATGGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATAtttt                                                      <  2:2959473/98‑1 (MQ=21)
                                agagGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATAtttt                                                      >  1:4093463/1‑101 (MQ=25)
                                agagGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATAtttt                                                      <  2:811852/101‑1 (MQ=25)
                                agagGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATAtttt                                                      <  2:607751/101‑1 (MQ=25)
                                agagGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATAtttt                                                      >  1:788967/1‑101 (MQ=25)
                                  agGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCg                                                    <  1:501055/101‑1 (MQ=25)
                                  agGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCg                                                    <  1:1119140/101‑1 (MQ=25)
                                  agGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCg                                                    >  1:2456427/1‑101 (MQ=25)
                                    ggggcgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGc                                                  <  1:3014581/99‑1 (MQ=21)
                                     cggcgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCa                                                 <  2:3052313/101‑1 (MQ=25)
                                       gcgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAg                                               <  1:557607/101‑1 (MQ=25)
                                        cgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGa                                              <  2:3823317/101‑1 (MQ=21)
                                        cgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGa                                              <  2:3681422/101‑1 (MQ=21)
                                         gAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGAt                                             >  1:1221764/1‑101 (MQ=21)
                                          aaGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATa                                            <  1:3047707/101‑1 (MQ=21)
                                          aaGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATa                                            >  1:3605514/1‑101 (MQ=21)
                                           aGATCGGCGCGGGCCGGAAGCGTCTGCCCGAATTAGGATAAAGGAAGAAGGGGGGCTCG.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATaa                                           <  2:768917/101‑1 (MQ=12)
                                           aGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATaa                                           <  1:1200835/101‑1 (MQ=21)
                                           aGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATaa                                           >  1:4141491/1‑101 (MQ=21)
                                            gATCGGCGCGGCGCTGAGGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGCCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATaaa                                          <  2:1431599/101‑1 (MQ=21)
                                            gATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATaaa                                          >  1:486256/1‑101 (MQ=21)
                                             aTCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATaaac                                         >  1:1796730/1‑100 (MQ=21)
                                             aTCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAg                                         <  1:3485493/101‑1 (MQ=21)
                                             aTCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAg                                         <  1:3877528/101‑1 (MQ=21)
                                              tCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAgg                                        <  1:1274122/101‑1 (MQ=21)
                                               cggcgcggcgCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGCTATTTTGGGCAAGATAAAGGc                                       >  2:1247673/1‑101 (MQ=21)
                                                 gcgcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGccc                                     <  2:1544582/101‑1 (MQ=21)
                                                  cgcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCg                                    <  1:4052966/101‑1 (MQ=21)
                                                  cgcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCg                                    <  2:2557665/101‑1 (MQ=21)
                                                   gcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCgg                                   >  1:4128181/1‑101 (MQ=21)
                                                   gcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCgg                                   <  2:3781012/101‑1 (MQ=21)
                                                   gcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCgg                                   >  1:1049718/1‑101 (MQ=21)
                                                    cggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGa                                  <  1:2770012/101‑1 (MQ=21)
                                                    cggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCAGCAAGATAAAGGCCCGGa                                  <  1:3444402/101‑1 (MQ=21)
                                                     ggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGaa                                 <  1:2541088/101‑1 (MQ=21)
                                                      gcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAg                                <  1:532589/101‑1 (MQ=21)
                                                       cgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGc                               >  1:3647200/1‑101 (MQ=21)
                                                        gcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGcc                              <  1:2819209/101‑1 (MQ=21)
                                                         cTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCt                             >  1:203511/1‑101 (MQ=21)
                                                         cTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCt                             <  1:232449/101‑1 (MQ=21)
                                                            aaGCGTCTGGGCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTccc                          <  2:1784789/101‑1 (MQ=21)
                                                             aGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCt                         <  2:921711/101‑1 (MQ=21)
                                                             aGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCt                         <  2:3976566/101‑1 (MQ=21)
                                                              gCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCtt                        <  1:1079300/101‑1 (MQ=21)
                                                              gCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCtt                        <  1:1367046/101‑1 (MQ=21)
                                                              gCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCtt                        <  1:3308880/101‑1 (MQ=21)
                                                               cGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCttt                       <  2:2407748/101‑1 (MQ=21)
                                                               cGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCttt                       >  1:180405/1‑101 (MQ=21)
                                                                 tCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTAt                     <  2:3529722/101‑1 (MQ=21)
                                                                 tCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTAt                     >  1:656153/1‑101 (MQ=21)
                                                                 tCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTAt                     <  2:777106/101‑1 (MQ=21)
                                                                  cTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATc                    <  1:23208/101‑1 (MQ=21)
                                                                   tGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCt                   <  2:3990292/101‑1 (MQ=21)
                                                                   tGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCt                   >  1:4024707/1‑101 (MQ=21)
                                                                    ggCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCtt                  <  2:4180641/101‑1 (MQ=21)
                                                                    ggCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCtt                  <  2:928703/101‑1 (MQ=21)
                                                                    ggCCGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCtt                  <  2:2252934/101‑1 (MQ=21)
                                                                       cGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGcc               >  2:2283685/1‑101 (MQ=21)
                                                                       cGAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGcc               <  1:1576435/101‑1 (MQ=21)
                                                                        gAACTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCg              <  2:317916/101‑1 (MQ=21)
                                                                         aaCTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCgg             <  1:1824468/101‑1 (MQ=21)
                                                                         aaCTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCgg             <  2:818523/101‑1 (MQ=21)
                                                                         aaCTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCgg             >  1:998422/1‑101 (MQ=21)
                                                                         aaCTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCgg             >  2:1283362/1‑101 (MQ=21)
                                                                         aaCTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCgg             >  2:1331295/1‑101 (MQ=21)
                                                                          aCTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGa            <  2:3588917/101‑1 (MQ=21)
                                                                           cTCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGaa           <  1:2539619/101‑1 (MQ=21)
                                                                            tCGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGaaa          <  1:3111348/101‑1 (MQ=21)
                                                                             cGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAat         >  1:299735/1‑101 (MQ=21)
                                                                             cGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAat         >  2:820012/1‑101 (MQ=21)
                                                                             cGGATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAat         <  1:1574579/101‑1 (MQ=21)
                                                                              ggATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAata        >  2:1254971/1‑101 (MQ=21)
                                                                               gATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAatat       >  2:3005797/1‑101 (MQ=21)
                                                                               gATAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAatat       <  2:464236/101‑1 (MQ=21)
                                                                                aTAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATc      <  2:2074130/101‑1 (MQ=21)
                                                                                aTAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATc      <  2:3557905/101‑1 (MQ=21)
                                                                                 tAAAGGAAGAAGGGGTGCTCT.CCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATcc     <  2:3727721/101‑1 (MQ=25)
                                                                                 tAAAGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAACATcc     >  2:3396319/1‑101 (MQ=21)
                                                                                  aaaGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATccc    >  1:997909/1‑101 (MQ=21)
                                                                                  aaaGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATccc    <  1:1490817/101‑1 (MQ=21)
                                                                                   aaGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGGTATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCa   >  1:3001598/1‑101 (MQ=21)
                                                                                   aaGGAAGAAGGGGGGCTCT.GCCCCCCGCTTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCa   >  1:431742/1‑101 (MQ=21)
                                                                                    aGGAAGAAGGGGGGCTCTGCCCCCCGCTTC.GCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCaa  <  2:1221668/101‑1 (MQ=21)
                                                                                    aGGAAGAAGGGGGGCTCTGCCCCCCGCTTC.GCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCaa  >  1:1100706/1‑101 (MQ=21)
                                                                                                     |                                                                                    
GAAGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTC..CCGCGCTTGTCGGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCAA  >  NC_014034/530839‑531022

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG