Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6610TG35.9% 34.4 / 10.7 42pseudogene (943/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (6/19);  new base G (14/0);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 3.34e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGCAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTC  >  NC_000913/3862430‑3862894
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑AGGTGCgggg                                                                                                                                                                                                                                             >  1:2641317/1‑238 (MQ=255)
     aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑AGTTGCggggaata                                                                                                                                                                                                                                         >  1:502812/1‑233 (MQ=255)
     aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGGTTA‑AGATGCggggaaata                                                                                                                                                                                                                                        >  1:2090977/1‑233 (MQ=255)
      tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑AGATGCGGGAAACAGc                                                                                                                                                                                                                                       >  1:852776/1‑239 (MQ=255)
                          gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGTTGCGGGGATCCGCAAG‑CCCGCCCCC‑TTTTA‑cc                                                                                                                                                                                                                 >  1:1554624/1‑238 (MQ=255)
                          gcgTGATTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTGA‑TGTTTCGGGGAAAAGCAAA‑GCCC‑CCAC‑CTTTA‑CCCg                                                                                                                                                                                                               >  2:1709366/1‑239 (MQ=255)
                           cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGGTGGGGGGATCAGCAAG‑CCCC‑CCCC‑ATTTT‑cccca                                                                                                                                                                                                              >  1:2253117/1‑237 (MQ=255)
                                ttGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGGTTCGGGGAACAACAAA‑CCCC‑CCCC‑ATTTT‑CCCGAtttt                                                                                                                                                                                                          >  1:1411540/1‑238 (MQ=255)
                                ttGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGGAATTCCCA‑CCAGCATGGCGGGTAA‑AGGTGGGGGGAGCAGCAAG‑GCCC‑CCAC‑TTTTT‑CCCgtattt                                                                                                                                                                                                          >  1:2185470/1‑233 (MQ=255)
                                                     ttttGCGGTAAAACCTAGTCAATCAGGCCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGTTGCg                                                                                                                                                                                                                                                <  1:1303400/183‑1 (MQ=255)
                                                        tGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGGTGCGGGGGACAGCCAG‑CCCC‑CCCC‑CTTTT‑CCCCACTTTCCGGTGCGGGC‑GTGAAGTACATCg                                                                                                                                                                                 >  1:1132809/1‑238 (MQ=255)
                                                           ggTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGc                                                                                                                                                                                                                                       <  1:830021/186‑1 (MQ=255)
                                                                     aGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCC‑CCAGCATGGCGGTTAA‑TGGTGCGGGGAGCATCAAG‑CCCG‑ccccc                                                                                                                                                                                                                         >  2:830021/1‑185 (MQ=255)
                                                                     aGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑AGGTGCGGGGAACAGCAAA‑CCCC‑CCCC‑TTTTT‑CCCCATTTTTCCGGGCGGCA‑GGGAAGGACATACCTaaa                                                                                                                                                                            >  1:2191410/1‑230 (MQ=255)
                                                                                                       acgaaaGAGATAAAATCGGGTGAACATGGGGGCCGCCGGCACGGGGTGCACCGAAAAAAGCAAGGGTAAAAGAAAAAGGAAGGTAATTC‑CCGGCGGTTCGGTAATCCCCA‑CCAGCGTGGCGTTTAA‑GGTTGGGGAAATCAACACC‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTAc                                                                                                                                                                                     <  2:1529024/184‑1 (MQ=255)
                                                                                                         aaaaGAGAAACATTACGGGCACAAATGACGCCGACAGGACGGCGTGACACACAAACACCAACGGAAAAAAAAAAAGGAGGAAAAATT‑CCGGCGTTTCGGTAACCCCCA‑CCTGCAGGGCGGTGAA‑GGTTGCGGGAGTCAACAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCa                                                                                                                             <  2:207745/238‑1 (MQ=255)
                                                                                                                 tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGGTGCGGGGATCAGCAAA‑CCCG‑CCAC‑CTTTA‑CCCCATTTTTCGGTGCGGCA‑GTGGAGTACAACGCTaaaa                                                                                                                                                                           >  1:1529024/1‑186 (MQ=255)
                                                                                                                 tACATTACGGTCGAAATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGGTTA‑GGGTGTGGGGATCAGCAAGACCCC‑CCCA‑TTTTC‑CCCGTTTTTCGGGGGCGG‑G‑GGGGAGGACAACGCCAAAAAAAAGCACCCCGCACCACAGACTTTTTTAGTT‑TCG‑TGGCGGGCAAAAACGc                                                                                                                      >  2:1797338/1‑238 (MQ=255)
                                                                                                                          gTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGGTAA‑AGGTTCGGGAAACAACCAA‑CCCC‑CCACTTTTTA‑CCCCTTTTTCC‑GTGCGGGAGGGGAGGAACTCCCCTAAGAAAATCCCCCCGCCACAAAAGA‑TTTTTTGTTTACAG‑T‑CAGGGGAAAACCGCCCTCACGGGa                                                                                                            >  2:179519/1‑239 (MQ=255)
                                                                                                                              aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGGTGCGGGGAACAGCAAA‑CCCG‑CCCC‑CTTTA‑CCCGAATTTTCGGTGCGGCA‑GGGAAAAAAATAGCCAAAGAAAAACACACAGGACCACAAA‑GTTTTTTGTTACCGCGGCAGGGAAAAACGCACTTC‑CGGGaacag                                                                                                        >  2:2094094/1‑238 (MQ=255)
                                                                                                                                    aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGTTGCGGGGATCAGCAAG‑CCCG‑CCCC‑TTTTA‑CCCGATTTTTCGGGGCGGCA‑GTGAAGTACATCg                                                                                                                                                                                 >  1:1548389/1‑162 (MQ=255)
                                                                                                                                         gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGCAATTCCCA‑CCAGCATGGCGGGTTA‑GGCTGCGGGAGACAGCAAA‑CCCC‑CCCC‑CTTTA‑CCCCATTTTTCCGTGCGGCA‑GGGGAGTAGAACGGTACAAAAAAACCCCCGGCACCCAAAA‑TTTTTTAGGTATCG‑TGCAGGGGAAAACCACCTTT‑CGGGGACACACACTTCAGCgg                                                                                            >  2:1232236/1‑239 (MQ=255)
                                                                                                                                                     tgaaCCGCAAACAGCAAGGGAAAAAGGACCAGGAAGGAAAATC‑CCAGCGGTTCGATAATCCCCC‑CCAGCAGGGGGTTTAG‑GGTTGGGGAAATCAGCAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTg                                                                                                                                                                                                    <  2:2484755/124‑1 (MQ=255)
                                                                                                                                                           gaaaaCAGCAACGGTGAAATGAAAAGGAAGGTAAATT‑CCAGCGGTTCGGTAATCCCCC‑CCAGATGGGGGGTTAA‑GGTGGGGGAAATCAGCAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGaa                                                                         <  1:1793387/237‑1 (MQ=255)
                                                                                                                                                               aCAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑AGTTTCGGGGATCAGCAAA‑CCCC‑CCAC‑TTTTA‑CCCCCTTTTCCCGTGCGGGA‑GTgga                                                                                                                                                                                         >  1:2484755/1‑125 (MQ=255)
                                                                                                                                                                       ggggAAAGGAAGAGGAAGGTAAATC‑CCAGGGGTTCGGTAATCCCCC‑CCAGCAGGGGGTTAAA‑GGTTGCGGAAACAAGCAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAAt                                                                                                                                                                       <  1:1576263/134‑1 (MQ=255)
                                                                                                                                                                        tgaaaaaGAACAGGAAGTAAAATC‑CCGGGGGTTGGGTAATCCCCC‑CCAGCTGGGCGGTTAA‑GGTGGGGGAAACAACCAGG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGAActct                                                                     <  1:1797338/227‑1 (MQ=255)
                                                                                                                                                                              agaacgatgaGGTAAGTA‑CCGGCGGTTCGTTACTCCCCC‑CCA‑CAAGGGGGTTAAGGGTGGCGAAAATCAGAGAC‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                      <  2:2185470/230‑1 (MQ=255)
                                                                                                                                                                                acaaggaaggaGACTC‑CCAGCGG‑TCGGTTAATCCCACCCAGCGGGGGGGTGAG‑GGTGGGGGAAACCAGCAAC‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                      <  2:2104407/234‑1 (MQ=255)
                                                                                                                                                                                 acgggaggtaaaaTC‑CCGGCGGTCCGTAAATCCCCC‑CCAGCAGGGCGGTTAA‑GGTGGCGGGAACCAGCAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                   <  1:1228077/229‑1 (MQ=255)
                                                                                                                                                                                 aaaaggatgTAAATTCCCGGCGGTTCGGAACTCCCCA‑CCAGCAGGGGGGTAAA‑GGTGGCGGAAACCAGCAAC‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                      <  1:1292254/229‑1 (MQ=255)
                                                                                                                                                                                 aCAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑AGGTGCGGGGAACAGCAAG‑CCCC‑CCAC‑ATTTACCCCGATTTT‑CGGGGGGGCA‑GGGAAGGAAATCCCTAAAAAAAAACCCACCGCa                                                                                                                                                             >  2:1576263/1‑137 (MQ=255)
                                                                                                                                                                                   gagagggTAAATC‑CCAGCGTTGCGGTAATCCCCC‑CCAGCATGGCGGTAAA‑GGTGGGGGGAAACAGCAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGgcgc                                                                                   <  1:2171979/200‑1 (MQ=255)
                                                                                                                                                                                            aaTT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGTTGCggg                                                                                                                                                                                                                                              >  1:41363/1‑50 (MQ=255)
                                                                                                                                                                                            aaTT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGTTGCggg                                                                                                                                                                                                                                              <  2:41363/50‑1 (MQ=255)
                                                                                                                                                                                               tcccgggggTTCGGTAATTCCCC‑CCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGCAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGa                                                                                                            <  1:2704102/163‑1 (MQ=255)
                                                                                                                                                                                               t‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGGTAA‑GGGTTGGGGAATCAACAAG‑CCCC‑CCAC‑TTTTA‑ACCGATTTTTCGGGGCGGCC‑GTGGAGGACAACGCTAAAAAAAAGCCCACCGGACCCAAGA‑TTTTTTAGTTACCC‑TGCAGGGCAAAACCGCCCTC‑CCGGGAAAACGATTTTAGCGGGTCCGGCGCTCAACCTAACCCCTTTCAGCTGCT‑CGccccagaccat                                             >  2:649339/1‑223 (MQ=255)
                                                                                                                                                                                                        tcgtgTAATCCCCA‑CCAGCAGGGCGGTTAA‑GGTGGGGGGAACCGGCAAC‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCgg                                                                                                                                                                                                 <  1:1131587/76‑1 (MQ=255)
                                                                                                                                                                                                         tCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGGTGCGGGGATCAGCCAA‑CCCC‑CCCC‑CTTTT‑CCCCATTTTTCGGTGCGGGC‑GGGGAGAAAAACGCCAAAAAAAAACCCACGGCACCAAAAA‑TTTTTTAGGTAACG‑TGGAGGGGAAAAACCCCCTC‑CGGGGAAAACGCttt                                                                                                  >  2:2704102/1‑170 (MQ=255)
                                                                                                                                                                                                                   cca‑ccaGCATGGGGTTTAG‑GGTGGGGGGAACGGGAAAA‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTa                                                                                                                                                                                      <  1:2166797/79‑1 (MQ=255)
                                                                                                                                                                                                                                tggTTAA‑GGTTGCGGGAATCAGCACC‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa       <  1:2094094/237‑1 (MQ=255)
                                                                                                                                                                                                                                     tagggTTGGGGGAATCAGCAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACTAAAGATAAACTGTc  <  1:1709366/236‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCA‑CCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGCAAG‑CCCG‑CCAC‑TTTTA‑CCCGATTTTCCGGTGCGGCA‑GTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGA‑TTTTTGAGTTACCG‑TGCAGGGCAAAACCGCCTTC‑CGGGAACAACGATTTCAGCGGCT‑CGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTC  >  NC_000913/3862430‑3862894

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