Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,2540GA31.0% 105.6 / 20.8 58S29L (TCA→TTA) emrKtripartite efflux pump membrane fusion protein EmrK
Reads supporting (aligned to +/- strand):  ref base G (29/11);  new base A (0/18);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.54e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTG  >  NC_000913/2483006‑2483504
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                          
ataGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGa                                                                                                                                                                                                                                                           >  1:1137002/1‑250 (MQ=255)
   gTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAAc                                                                                                                                                                                                                                                       <  2:1010192/251‑1 (MQ=255)
   gTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAAc                                                                                                                                                                                                                                                       >  1:454284/1‑251 (MQ=255)
         tCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTGTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑GTGACC‑AATAGGCATAGGCACC‑TGAAAACGCaat                                                                                                                                                                                                                                                  <  2:985906/250‑1 (MQ=255)
                ccGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAAtaaataaa                                                                                                                                                                                                                                           >  1:175363/1‑250 (MQ=255)
                       atatTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGc                                                                                                                                                                                                                                   <  2:741647/251‑1 (MQ=255)
                         atTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCa                                                                                                                                                                                                                                 <  1:15538/251‑1 (MQ=255)
                              ccAGATTATTTTCAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATaa                                                                                                                                                                                                                             >  1:893154/1‑250 (MQ=255)
                                     attTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGt                                                                                                                                                                                                                     >  1:1120063/1‑251 (MQ=255)
                                                 atggAGTGCGATAGTGGCATCATTTTTATCCAGTGAAACAAAAATGTCCCCTTGTCGAAC‑TTAGTTCGTATTTTTATGATAAACGACAGTGACCCTACCTG‑AGACTTGT‑GCAGAAATGGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCCCC‑GGAAACCGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGt                                                                                                                                                                                                         <  1:1184568/248‑1 (MQ=255)
                                                    gAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTAGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                                                                                                                        >  1:661648/1‑249 (MQ=255)
                                                     aGTGCGATAGCGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                                                                                                                        >  2:302885/1‑248 (MQ=255)
                                                              gTGGCTTCATTTTTATCCAGAGAAACAAAAATGTCACCTTGTCGAAC‑TTATTTCTTTTTTTTATAATTAACGCCAGTGACCCTCCCGG‑AGATTTGT‑GCAAAAATTGGATCTGCATCCCCCGTGACATAGCCGTCATCTTTACTAATCTTTTCTTCTA‑TTCCAATCGCCA‑AATAGACATAGCCCCC‑AAAAACCCCAATAAAAAAAACTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttttt                                                                                                                                                                                            <  1:1096298/251‑1 (MQ=255)
                                                                ggaaTCAGTTTTATACAATGAAACAAAAAGGTCACCTTTGCGAAC‑GTAGTTCTTATTTTTATGATAAACGACATTGACATGCCCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATCCCCCGTGACAAAGGCGTAATCTGTACTAATCATTTCTTAAAATTCCA‑TTGCCC‑AAAAGAAAAAGGCACA‑TAAAAACTAAATAAAAAAATCTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAt                                                                                                                                                                                          <  2:12315/248‑1 (MQ=255)
                                                                   aTCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTg                                                                                                                                                                                       >  1:360154/1‑251 (MQ=255)
                                                                      aGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAt                                                                                                                                                                                    >  1:823944/1‑251 (MQ=255)
                                                                           tATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATc                                                                                                                                                                               >  2:550585/1‑251 (MQ=255)
                                                                            aTCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCt                                                                                                                                                                              >  2:60596/1‑251 (MQ=255)
                                                                            aTCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCt                                                                                                                                                                              >  1:60717/1‑251 (MQ=255)
                                                                            aTCCAGTAAAACTAAAATTTCACTTTGTCGAAC‑GTAGTTCGTATCTTTATGTTTACCGACAGTACCACTCCCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATCGGCCTCTTCTGTTCTAACCATTTCTTCTAATCCCA‑TAGCCC‑AATAGGCATAGGCCCC‑TAAAACCGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATc                                                                                                                                                                               <  1:135337/250‑1 (MQ=255)
                                                                             tCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTg                                                                                                                                                                             >  2:775701/1‑251 (MQ=255)
                                                                                aGTGAAACTAAAATTTCACCTTGTGGAAC‑TTGTTTCGTTTTTTTATTATTACCACCAGTGACTCTCCCTG‑ATACTTGT‑GCAGAAATTGGTTCTTCATTCCCCGTGACTAAGGCTTCATTTGTCCAAATCATTTCTT‑TAAATCCA‑ATGCCCCAATAGGTAAAGACACG‑GAAAAACGAAATAAATAAATTTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTc                                                                                                                                                                          <  2:139000/251‑1 (MQ=255)
                                                                                    aaaCTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCATGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATaa                                                                                                                                                                                                                             >  1:1078244/1‑196 (MQ=255)
                                                                                       atAAAATGCCACCTTGTCTACC‑GTAGTTCATTTTTTTATGATTAACTACATGGACTCTCCCT‑‑AAAATTGTTGCAAAAATTGGATCTGCATCCCCTGGGACATGGGCTTCATTTGAACTAAGCATGTTTTCTAATCCCA‑TGGCCC‑AA‑AAGCCTGAGCACCCTAAAAACGAAATAAATAAAACAACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtatt                                                                                                                                                                   <  1:1176947/250‑1 (MQ=255)
                                                                                            aTGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTCCCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtctc                                                                                                                                                              >  1:156794/1‑251 (MQ=255)
                                                                                             tGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtctct                                                                                                                                                             >  2:791820/1‑251 (MQ=255)
                                                                                              gtccccTTGTCGAAT‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTTT‑GAATAAATTGGATCTGCATCCCCCGTAACATAGGCGTCATCTGTACTAATCATGTCTTATAATTCCA‑TGGACC‑AATAGACAAAGGCACC‑TAAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtctct                                                                                                                                                             <  2:893154/246‑1 (MQ=255)
                                                                                               ccACCTTGACGAACTGT‑GTTCGTATTTTAATGATTAACGACAGTACCACTACGTG‑AGACTTGT‑GAAGAAATTGGGTCTGCATTCCCCGTGAAATAGGCGTCATCTGTACTAATCATGTCTTATAATTCCA‑TTGACA‑AATAGCCAAGGGCCCC‑TGAAAACGCAATAAAAAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCa                                                                                                                                                           <  1:178334/250‑1 (MQ=255)
                                                                                                      gTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtctc                                                                                                                                                     >  1:1083183/1‑250 (MQ=255)
                                                                                                         gAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCAt                                                                                                                                                   >  2:1184102/1‑249 (MQ=255)
                                                                                                                 gTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAAt                                                                                                                                          >  1:508376/1‑251 (MQ=255)
                                                                                                                      tATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttt                                                                                                                                                                                              >  2:307897/1‑194 (MQ=255)
                                                                                                                      tATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAAACATGTCTTCTAATTCCA‑TTGACA‑AATGGGCATAGGCCCC‑GAAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttt                                                                                                                                                                                              <  1:307897/194‑1 (MQ=255)
                                                                                                                          tttATGATTAACGACAGTGACACTCCCTG‑AGACTTGT‑GCAAAAATTGGATCTGCTTCCCCCGTTACATAGGCGTCATCTGTACAAATCATGTTTTCTAATTCCA‑TTGACC‑AATAGGCATAGCCACC‑TGAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaat                                                                                                                                  <  1:620288/250‑1 (MQ=255)
                                                                                                                           ttATGATTAACGACAGTGACCCTCCCTG‑AGACTTGT‑GCAAAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCCCC‑TGAAACCGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaatt                                                                                                                                 <  1:791820/250‑1 (MQ=255)
                                                                                                                           ttATGATTAACGACAGTGACACT‑CCTGGAGACTTGT‑GCAGAAATTGAATCGGCATTCCCCGTGACATAGGCGTTTTCTGTACAAATTATGTTTTCTAATGCCA‑TTGACC‑AAAAGATAGGGGCACC‑TGAAAACGCAATAAAAAAATCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaatt                                                                                                                                 <  1:447731/250‑1 (MQ=255)
                                                                                                                            tatgattacgaacattgaccctcccTG‑AGATTTGT‑GCAAAAATTGGATCTGCTTTCCCCGTAACAAGGGTGTCATCTTTACTAATCATGTCTTATAATTCCA‑TAGACC‑AAAAGGAAAAGACACC‑TAAAAACGCAATAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattg                                                                                                                                <  2:360154/227‑1 (MQ=255)
                                                                                                                               aaTTAACGACAGGGACACTCCCTT‑GGTCTTGT‑GCAAAAATGGGATCGGCTTCCCCCGGGTCATGGGCTTTATCTGTACTAATCATGTCTTCAAATCCCA‑TGTCCC‑AAAAGGCATAGCAACC‑TAAAAGCGGAATAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCTCTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattgaat                                                                                                                             <  1:1085399/249‑1 (MQ=255)
                                                                                                                                           tGACACTACCTG‑AGACCTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGAAC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcattc                                                                                                                >  1:384116/1‑251 (MQ=255)
                                                                                                                                                 tACCTG‑AGACTTGT‑GCAGATATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATaaa                                                                                                                                                                                                                            >  1:997143/1‑137 (MQ=255)
                                                                                                                                                  aCCTG‑AGACTGGG‑GCAGAAATTGGATCCGCATCCCCGGTGACATAGGCCTCATCTGTACTAACCATGTCTTTTAATTCTA‑TGGACA‑AAGAGGCAAAGCCCCC‑TAAAAACGCAAAAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCAgtcttc                                                                                                         <  1:775701/251‑1 (MQ=255)
                                                                                                                                                           cTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGtt                                                                                                                                                                           >  1:775873/1‑177 (MQ=255)
                                                                                                                                                             tgt‑gCAAAATTTGGTTCTTCATTCCCCGTGAAATAGGCGTATTCTGTACTAATAATTTCTTCAAATACCA‑TTGACA‑AATAGACAAAGGCAAC‑TAAAACCGAAAAAAATAAACTTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTTAATGAAAAGAGCATTCATTCTTCATATCAAac                                                                                                <  2:661648/250‑1 (MQ=255)
                                                                                                                                                                  cAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGt                                                                                           >  2:1070457/1‑251 (MQ=255)
                                                                                                                                                                            atctccttcccccGTGACAGAGGCGTCATTTGTACTAATCATGTCTTCTAATCCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGtt                                                                                  <  2:60717/241‑1 (MQ=255)
                                                                                                                                                                               tGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGcac                                                                              >  1:452052/1‑251 (MQ=255)
                                                                                                                                                                                gcttcccccGTGACAGAGGCGTCATTTGTACTAATCTTGTTTTCTAATCCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGtt                                                                                  <  1:698658/241‑1 (MQ=255)
                                                                                                                                                                                gCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCCCC‑TAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGc                                                                                <  1:175301/248‑1 (MQ=255)
                                                                                                                                                                                gCATCCCCCGTGACATAGGCGTCATCTGTCAAAATAATTTCTCCAAATTCCA‑TTGACC‑AATAGGAAAAGGCACC‑TAAAAACGAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGc                                                                                <  2:1120063/248‑1 (MQ=255)
                                                                                                                                                                                       ccGTGACAAGGGCTTCATCTGCACTAATCATGTCTTCTATTACCA‑TTGCCC‑AATAGGCATAGTCACC‑TAAAAACGCAAAAAATAAATCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATccc                                                                      <  1:505181/251‑1 (MQ=255)
                                                                                                                                                                                                    gtcatttttacaactcatttcatatAATTCCA‑TTGACA‑AATAGGCATAGGCACT‑TAAAAACGTAATAAAAAAATCTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAg                                                         <  2:175363/227‑1 (MQ=255)
                                                                                                                                                                                                    gTCATCTGAACAAATCATTTCTTTAAATTCCA‑TTGACC‑AATAGACATAGGCACC‑TAAAACCGAAAAAAAAAAATCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAg                                                         <  2:320205/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ccA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGAATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAaac                            >  2:129008/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cc‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACag                       >  1:1126370/1‑249 (MQ=255)
                                                                                                                                                                                                                                            aaTAGGAAAAGGCACC‑TAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGa                   <  1:1192460/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggcaCC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTg         >  2:771189/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        caCC‑TAAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGc       <  1:570374/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           c‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCttt    >  1:899106/1‑251 (MQ=255)
                                                                                                                                                                                                                                                              gAAAACGCAATAAATAAAACTACCGCCAATAAAAAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCACAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTg  >  2:964428/1‑251 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                          
ATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAAC‑GTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTG‑AGACTTGT‑GCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCA‑TTGACC‑AATAGGCATAGGCACC‑TGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTG  >  NC_000913/2483006‑2483504

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