Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,909,6020GA59.6% 34.1 / 39.6 49A22V (GCA→GTA) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base G (8/11);  new base A (15/13);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTCTTTACGCAGAACGTCAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAAC‑GCTCGTAGTGCCCCAGGTCCAGGT‑CGGTTTCAGCGCCGTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAA  >  NC_000913/2909356‑2909842
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
cgTTCTTTACGCAGAACGTCAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAAC‑GCTCGTAGTGCCCCAGGTCCAGGT‑CGGTTTCAGCGTCGTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGc                                                                                                                                                                                                                                              >  1:126891/1‑251 (MQ=255)
           cAGAACGTCAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGTTCATTTTGGTACGAATGAAAC‑GCTCGTAGTGCCCCAGGTCCAGGT‑CTGTTGCAGCGCCGTCTTCAGTAACGAACGC‑TTCCCCGTGGGGGATGGGGCTCACAGTGCCTGGATCGACTTTGATGGAC‑GGATCCAGTTTCATGATGGTCACATTGAGATCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTAc                                                                                                                                                                                                                                   <  1:574439/251‑1 (MQ=255)
             gAACGTGAGAGTAGATCCGACCCGTGGTGAAGTTGTTGCGTCGGCTCATTTTGGTACGAATGGAGC‑GCTCGTAGGTCCCCAGGGCCAGGT‑CGGTTTCAGCGCAGTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGACGCATAGTACTTGGATCGACGTTTATGTAC‑GGATCCAGTTTCATGATGGGCACATTTAGGCCACGGGCTTCAAGAACGTCGACGAGTGAGAGTGCTGCAATGCCTTTAccc                                                                                                                                                                                                                                 <  2:511147/251‑1 (MQ=255)
                   cAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAAC‑GCTCGTAGTGCCCCAGGTCCAGGT‑CGGTTTCAGCGCCGTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCagaga                                                                                                                                                                                                                            <  1:160298/250‑1 (MQ=255)
                   cAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAAC‑GCTCGTAGTGCCCCAGGTCCAGGT‑CGGTTTCAGCGCCGTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCagaga                                                                                                                                                                                                                            <  2:402055/250‑1 (MQ=255)
                                     gggtgAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAAC‑GCTCGTAGTGCCCAAGGTCCAGGT‑GGGTTTCAGCGCCGTCTTCAGTAACGAACCC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCGGGATCGACGTTGATGTAC‑GGATCCATTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGCCCCCGCCGGTc                                                                                                                                                                                                         <  1:366145/250‑1 (MQ=255)
                                                  gcggcgGCTCATTTTGGTACGAATGAACC‑GCTCGTAGTGCCCCAGGGCCAGGT‑CGGTTGCAGCCCCGTCTTCAGTAACGAACGC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATTTAC‑GGATCCATTTTCATGGGGGTCACATTGAGTCCACGTGCTTCAAGATTGGCTGCGAGGGAGGCTACTGTAGTGCGTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAg                                                                                                                                                                                             <  1:44152/250‑1 (MQ=255)
                                                   cggcgcttCATTTTGGTACGAATGAAAC‑GCTCGTAGTGCCCCAGGTCCAGGT‑CGGTTTCAGCGCCGTCTTCAGTAACGAACCC‑TTCCCCGTGTTGGATTGGGCTCATAGTCCCTGAATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAAAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGtt                                                                                                                                                                                           <  1:137870/244‑1 (MQ=255)
                                                    ggcggcTAATTTTGTTACAAATGAAAC‑CCTCGTAGTGCCCC‑GGTCCCGGTTGGGTTTAAGCGCCTTCTTCAGTAACGAACCC‑TCCCCCGTGTGGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GAATCCATTTTCAGGATGTTCACATTGAGGCCCCGGGCTTCAAGAAGGGCTGCGAGGGAGGCTACGTAAAGCCCTTTCCCCAGAGAGAATACGACCCCGCCGGTCACAAAAATATAGTTc                                                                                                                                                                                          <  2:126891/251‑1 (MQ=255)
                                                        gcTCATTTTGGTACGAA‑GAAACCGCGCGTAGTGCCCCAGGTCCAGGT‑CGGTTTCAGCGCCGTCTTCAGAAACGAACCC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCGGGATCGCCGTTGATGTAC‑GGATCCAGTTTCAGGATGGTCACATGGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGAAATGCTTTTACCCAGAGAGGATACGCCCCCGCCGGTCACAAAAATATAGTTCGTTg                                                                                                                                                                                      <  2:268621/251‑1 (MQ=255)
                                                                      gAATGAAAC‑GCTCGTAGTGCCCCAGGTCCAGGT‑CGGTTTCAGCGCCGTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTg                                                                                                                                                                        >  1:71163/1‑251 (MQ=255)
                                                                             cc‑GCTCGTAGTGCCCCAGGTCCAGGT‑CGGTTTCAGCGCCGTCTTCAGTAACAAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTCCCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGCCCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGtt                                                                                                                                                                 <  1:596681/250‑1 (MQ=255)
                                                                              ccgcgcGTAGTGCCCGAGGCCCAGGT‑CGGTTTCCGCGCCGTCTTCGTTAACGAACGC‑GTCCCCGGGTGGGATTGGGCTCAAAGTACCTGGATCGAGGTGGA‑GTCCCGGATCCAGTTTCATGATGGTCACATTGAGGCCGCGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGAAATGCCTTTCCCCAGATAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGtt                                                                                                                                                                 <  2:165051/246‑1 (MQ=255)
                                                                                    gTAGTGCCCCAGGCCCAGGT‑CGGTTTCAGCGCCGTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCccgcc                                                                                                                                                                                                             <  2:112319/201‑1 (MQ=255)
                                                                                      aGTGCCCCAGGTCCAGGG‑CGGTTTCAGCGCCGTCTTCAGTAACTAACAC‑TTCCCCGTGTGGGATTGGGCTCATAGTACCTGGAGCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCGGCGGCAATGCCTTTACCCAGAGAGGATCCAACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGa                                                                                                                                                          <  1:60674/250‑1 (MQ=255)
                                                                                              gggtccaggt‑cGGTTTCAGCCCCGTCTTCAGTACCGAACAC‑TTCCCCGGGGTGGATGGGGCTAATAGTACCTGGATCGGCGTGGATGTAC‑GGATCCTGTTTCATAAGGGTAACATGGAGGCCACGGGCTACAAGAATGGCTGGGAGGGGGGCTACGGAAATGCCTTTACCAAGAGAGGAGACGACCCCGCCGGTCACAAAAATATAGGTCGTTGTCATGCTGAACCTGAGAAGTTAGGTGGAAAAGACGAt                                                                                                                                                 <  2:391766/250‑1 (MQ=255)
                                                                                              aggtccaggt‑cGGTTTCAGCGCCGTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCACAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGAt                                                                                                                                                 >  1:591558/1‑251 (MQ=255)
                                                                                                gtccaggt‑cGGTTTCAGCGCCGTCTTCAGTAACGAACAC‑TTACCCGTGTTGGATTGGGCTCATAGTAACTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATgg                                                                                                                                               >  1:407549/1‑251 (MQ=255)
                                                                                                     ggt‑cGGTTTCAGCGCCGGCGTCAGTAACGACCAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGTTGAACCTGAGAAGTTAGGTTGCAAAGACGATGGAAtg                                                                                                                                           >  2:415534/1‑249 (MQ=255)
                                                                                                        ggggTTTAAGCGCCGTCTTCAGTAACGAACCC‑TTCCCCGTGTTGGATTGGTCTCATAGTACCTGGTTCGACGTTGAGGTGC‑GGAGCCAGTTTCAGGAGGGTCACATTGGGGCCACGGGCTCCAAGAAGGGTGGCGGGGGGGGCTACGGAAAGGCCGTTCCCCAGAGAGGTTGCCACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAAcc                                                                                                                                        <  2:400607/249‑1 (MQ=255)
                                                                                                          ggTTTCAGCGCCGTCTTCAGTAACGAACAC‑TTCCCCGTGTCGGATTGGGCCCCTAGTACCTGGATCGACGCTGATGTAC‑GGATCCCGTTTCATGATGGTCACCTTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCACTGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAACAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCa                                                                                                                                       >  1:576025/1‑250 (MQ=255)
                                                                                                           gttccggccccctCTTCAGTAACGA‑CACCTCCCCCGGTTTGGATTGGGCTAAAAGTACCTGGTCCGGCGTTGATGTAC‑GGATCCAGTTTAATGATGGTCACATGGAGGCCCCGGGCTTCAAGAAGGGCTGCGAGGGAGGCTACGGAAATGCTTTTCCCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAg                                                                                                                                     <  2:9338/239‑1 (MQ=255)
                                                                                                               cAGCGCCGTTTTCAGTAACGAACAC‑GTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTCC‑GGATCCAGTTTCATGATGGTCACATTGAGCCCACGGGCTTCAAAAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTCCCCATAGAGGATACGCCCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACgg                                                                                                                                 <  1:89209/251‑1 (MQ=255)
                                                                                                                      gTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTCC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGc                                                                                                                           >  2:193892/1‑250 (MQ=255)
                                                                                                                        cTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGt                                                                                                                        >  2:64714/1‑251 (MQ=255)
                                                                                                                        cTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCAAGTTTCATGATGGTCACATTGAGGCCCCGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGCAGTTTGGTTGAAAAGACGATGGAATAACTAAGCAGGGAAAGCAGt                                                                                                                        >  2:628/1‑251 (MQ=255)
                                                                                                                           cAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGGGCAGAGTACCTGGATCGACGTTGAAGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTAGCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGtat                                                                                                                      >  1:91448/1‑250 (MQ=255)
                                                                                                                                    acac‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACAt                                                                                                            >  1:297616/1‑251 (MQ=255)
                                                                                                                                      ac‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATg                                                                                                           >  1:117984/1‑250 (MQ=255)
                                                                                                                                             ccGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCtgtg                                                                                                    >  1:2728/1‑251 (MQ=255)
                                                                                                                                                   tGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCcaca                                                                                               >  1:172335/1‑250 (MQ=255)
                                                                                                                                                     gATTGGGCTAATAGTACCTGGATCGGCGTTGGTGTCC‑GGATCCAGTTTCAGGAGGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTCCCCAGAGAGGATACGCCCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAAc                                                                                            <  2:117984/251‑1 (MQ=255)
                                                                                                                                                         gggCTCATAGTACCTGGATCGCCGTTGATGCAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCACGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCAAAAACTTTc                                                                                        >  2:563792/1‑251 (MQ=255)
                                                                                                                                                                                 ttGATGTCC‑GGATCCAGTTTCATGAGGGTCACATTGAGGCCACGGGCTTCAAGAATGGCGGCGGGGGAGGCTCCGGCAATGCCTTTACCCAGAGAGGATACGCCCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAAGAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAgc                                                                  <  1:513838/249‑1 (MQ=255)
                                                                                                                                                                                  tGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTACTCTTCGCCAgcgc                                                                >  2:160836/1‑250 (MQ=255)
                                                                                                                                                                                       tAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACCAAAATATAGTTCGTTGTCATGCTGAACCTGAGACGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAACCTTTCATTATCCCTCCTCTTCGCCAGCTCACTAtt                                                          >  2:130918/1‑251 (MQ=255)
                                                                                                                                                                                       gaC‑GGATCAGGTTTCATGATGGCCACAGGGGCGCCACGGGCTTAAAGAATGGCCGGGAGGGAGGCTGCGGCAATGCCTTTACCCAAAGAGGAGACACCCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAATCTTTCATTATCCCTCCTCTTCGCCAGCGCACTAtt                                                          <  2:256069/250‑1 (MQ=255)
                                                                                                                                                                                          cgacccgaGTTTCATGAGGGTCACATTGAGGCCACGGGCTTCAAGAATGGCGGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCCCTTCGCCAGCGCACGATTGa                                                        <  1:193892/244‑1 (MQ=255)
                                                                                                                                                                                           ggATCCAGTTTCAGGAGGGCCACATGGAGGCCACGGGCTACAAGAATGGCTGCGGGGGGGGCTGCGGCAATGCCTTTCCCCAGAGAGAATACGCCCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGaa                                                       <  2:297616/251‑1 (MQ=255)
                                                                                                                                                                                            aatccattttaaggatgGCCACATTGAGGCCACGGGCTTCAAGAAGGGCGGCGAGGGAGGCTACGCCAATGCCTTTGCCCAGAGAGAATCCGCCCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGaaa                                                      <  2:336849/238‑1 (MQ=255)
                                                                                                                                                                                             aTCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAg                                                                                                                                                      >  1:331166/1‑154 (MQ=255)
                                                                                                                                                                                                  gTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAgg                                                 >  1:45923/1‑250 (MQ=255)
                                                                                                                                                                                                    ttCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAgg                                              >  1:23065/1‑250 (MQ=255)
                                                                                                                                                                                                           gggTCATATAGAGGCGACGGGCGTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCTTTTCCCCAGAGAGGTGACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCGGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCTCCAGCGCACTATTGAAATCAGGAGTGAGaaaa                                       <  2:574005/250‑1 (MQ=255)
                                                                                                                                                                                                               cacaTTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGc                                   >  1:259904/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tCAAGAATGACTGCGAGGGAGGTTCTGGCAATGCCTTTACCCCGAGAGGATACGCCCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGtt                 <  1:37276/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gAATGGCTGCGAGGGAGGCTACGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTg             >  1:47603/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gCTGCGAGGGAGGATGCGGCAATGCATTTACCCAGAGAGGATACGAACCAGCCGGTCACACAAATATAGTTCGTTGTCATGCTGAACCTGAGACGTTAGGTTGACAAGACGATGGACTAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCa        >  2:268857/1‑250 (MQ=255)
                                                                                                                                                                                                                                                aagggagggTGCGGCATGGCCTTGACCAAGAGGGGAAACGCCCCCGCCGGTCACAAAAATATAGTACGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTACATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTGTGACGCAAATCaa  <  2:246195/250‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
CGTTCTTTACGCAGAACGTCAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAAC‑GCTCGTAGTGCCCCAGGTCCAGGT‑CGGTTTCAGCGCCGTCTTCAGTAACGAACAC‑TTCCCCGTGTTGGATTGGGCTCATAGTACCTGGATCGACGTTGATGTAC‑GGATCCAGTTTCATGATGGTCACATTGAGGCCACGGGCTTCAAGAATGGCTGCGAGGGAGGCTGCGGCAATGCCTTTACCCAGAGAGGATACGACCCCGCCGGTCACAAAAATATAGTTCGTTGTCATGCTGAACCTGAGAAGTTAGGTTGAAAAGACGATGGAATAACCAAGACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAA  >  NC_000913/2909356‑2909842

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑