Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,2680TA31.2% 83.4 / 18.6 48L24F (TTA→TTT) emrKtripartite efflux pump membrane fusion protein EmrK
Reads supporting (aligned to +/- strand):  ref base T (26/7);  new base A (0/15);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.56e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTG  >  NC_000913/2483021‑2483504
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
gCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATATGCGTCATCTGTACTAATCATGTCTGCTAAT‑TTCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAataaataa                                                                                                                                                                                                                                             <  1:394241/250‑1 (MQ=255)
   gAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACt                                                                                                                                                                                                                                         >  1:793341/1‑251 (MQ=255)
    aaCAATATTTGCCAGATTATTTTTAGCTTTAATGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTGCTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑CCAATAAATAAAACT‑a                                                                                                                                                                                                                                       <  2:600024/251‑1 (MQ=255)
     aCAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCATTTTTATCCAGTGAAACTAAAATGT‑CGCCTTGTCGAACGTAGTTCGTATCTGTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGGGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATGTACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑Ac                                                                                                                                                                                                                                      <  1:913255/251‑1 (MQ=255)
     aCAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTCCAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGCCTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑Ac                                                                                                                                                                                                                                      >  1:876269/1‑251 (MQ=255)
                      attTTTAGCTTTATTGGGTGCGATAGTGGCATCAGTTTTATCCAGTGAACCTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATATGCATTCCCCGTGACATCGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGt                                                                                                                                                                                                                     >  2:576068/1‑251 (MQ=255)
                          ttAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTtttt                                                                                                                                                                                                                 >  2:635522/1‑251 (MQ=255)
                              cTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTttctt                                                                                                                                                                                                             >  2:249129/1‑251 (MQ=255)
                                 tattgattcgaatattggaaTCAGTTTTATCCATTAAAACTAAAATGT‑CCCTTTGTCGACCGTAGTTCTTATTTTTA‑TGATTAACGACAGGGACACT‑CCTGGAGCCTTGTGAAGAAATTGT‑ATCTGCATCCCCCGTGCCAGAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCTTTGACCAA‑TGGGCATAGCCCCCTAAAAAC‑GAAATAAATAAACCT‑ACCGCCAATAAAGAAAAGTATTTtcttct                                                                                                                                                                                                           <  1:668819/231‑1 (MQ=255)
                                  aTTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGACACTAAACTGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGt                                                                                                                                                                                                         >  2:470719/1‑251 (MQ=255)
                                     tccagccgtgtgctctccctatttACCCAGTGAAACAAAAATGT‑CCCTTTGCCGAACTTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTCCCGG‑AGACTTGGGAAAAAATGGG‑ATCTGCATCCCCCGTGACAAGGGTGTCATCTGTACTAATATTGTCTTCAAAT‑TCCATGCACCAA‑TAGGCAAAGGCACCTAAAAAC‑TCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTa                                                                                                                                                                                                       <  1:948118/229‑1 (MQ=255)
                                       gTGCGATAGGGGCATCAGTTTAACCCAGTGAAACAAAAATTC‑CCCCTGGTCGAACGTGGTTCGTTTCTTTA‑GAATAAACGACAGTGACCCTCCCGG‑AGATTTGTGCAAAAATTGG‑ATCGGCATCCCCCGGGACAAAGGCTTGATCTTTACAAATGATGTCTTC‑AATGTCCATTGACCAA‑GGGGAAAAGGCCCCTAAAACC‑GCAAAAAAAAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGaa                                                                                                                                                                                                    <  1:470719/251‑1 (MQ=255)
                                              aGTGGCATCAGTTTTACCCAGTGAAACAAAAATGT‑CACTTTGTCGACCTTAGTTCGTATCTTTA‑TGATTAACGACAGTGACTCTCCCGG‑AGACTTGTGCAAAAATTGG‑ATCTGCATTCCCCGTAACAGAGGCGTCATCTGTACTAATCTTGTTTTCTAAT‑CCCATTGCCCAA‑TAGGCATAGGCACCTGAAAAC‑GAAAAAAAAAAACCT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttt                                                                                                                                                                                              <  1:442041/250‑1 (MQ=255)
                                                 ggCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTa                                                                                                                                                                                                       >  1:754499/1‑238 (MQ=255)
                                                  gCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAtt                                                                                                                                                                                         >  2:656569/1‑251 (MQ=255)
                                                    aTCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTg                                                                                                                                                                                       >  2:434074/1‑251 (MQ=255)
                                                         ttttACCCAGTGAAAATAAAA‑GTGCCCCTTGTCGACATTAGTTC‑TATTTTTAATGATTACCGACAGTGACACAACCGG‑AGACTTGTGAAAAATTTTG‑ATCTGTATCCCCCTTGACAGAGGCTTTATGTGTACTAATCTTGTCTCCTAAT‑CAAATTGCCCAA‑TAGGCAAAGGCAC‑TAAAAACCGCAA‑AAAAAAAACTTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                                                                                                                                   <  2:801160/250‑1 (MQ=255)
                                                         ttttACCCAGTAAAACAAAAATGT‑CCCCTTGTAGAACGTAGTTCGTTTTTTTA‑TGATAAACGACATTGCCACTACCTG‑AGACTTGGGAAGAAATTGG‑ATCGGCATCCCCGGTGACGTGGGCGTCATCTGAATTAATCATGTTTTTTAAT‑CCAATTCACCAA‑AGGGAAGAGGCACCTGAAAAC‑CAAATAAATAAAACT‑TCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                                                                                                                                   <  1:434074/250‑1 (MQ=255)
                                                           ttATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAAt                                                                                                                                                                                >  2:30794/1‑251 (MQ=255)
                                                             aTCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTa                                                                                                                                                                                                       >  2:948118/1‑226 (MQ=255)
                                                                    tgaaataaaaagctccccTTGTCGAACGTAGTTCGATTCTTTA‑TGATAAACGACAGTGCCACTCCGTG‑AGCTTTGGACAGAAATGGG‑ATCGGAATCCCCCTTAAAATAGGCTCCATCTATACTAATCATGTCTTCTAAT‑TCCATT‑ACCAAATAGGCATAGACACCTAAAAAC‑TCAAAAAAAAAAACT‑TCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTcc                                                                                                                                                                         <  2:159909/234‑1 (MQ=255)
                                                                    gAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCAc                                                                                                                                                                       >  1:541365/1‑251 (MQ=255)
                                                                    aaaaaaaaaaatgtccctttgcagAACGTAGTTCGTATCTTTG‑TGATTAACGACAGTGCCTCTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTCCATTCCCCGTGACACAGTCGTCATCTGTAATAATAATTTCTTCTAAT‑CCCATTCACCAA‑GAGACAAAGGCACTTAAAAAC‑TCAATAAAAAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCa                                                                                                                                                                        <  2:155289/228‑1 (MQ=255)
                                                                      aaacaaaatggcctctttGTCAAACGTAGTTCGTTTCTTTA‑TGATTAACGACAGTGCCTCTACCTG‑AGGTTTGTGCAAAAATGGG‑ATCTGCATTCCCCGTGACATAGGCGTCTTCTTTTCAAATCATGTCTTCTATT‑CCCTTTGCCAAA‑GAGGCATAGGCACCTAAAAAC‑GAAATAAATAAACCT‑GCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACt                                                                                                                                                                      <  1:635522/235‑1 (MQ=255)
                                                                       aCTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGAATTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGAACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtat                                                                                                                                                                    >  1:499339/1‑251 (MQ=255)
                                                                                                             aa‑TGATTAACGACATTCACACTACCTT‑AGTCTTGAGAAGAAATTGT‑ACCTCCATCCCCCGTGAAATAGGTGTCATCTGTTCAAATCATGTCTTCTAAT‑TCCATTCACCAA‑TAGTCAAAGCCACCTAAAACC‑GCAATAAAAAAACCT‑ACCGCCAATAAAGAAAAGTATTTTCTACTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattg                                                                                                                                <  2:279268/249‑1 (MQ=255)
                                                                                                                 gATTAACGAAAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGACTGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATtgaatg                                                                                                                            >  2:585423/1‑251 (MQ=255)
                                                                                                                     accgaaagggaccctcccTG‑AGCCTTGTGCATAAATTGG‑ATCTCAATCCCCCGTGACATGGCCGCCTTTTGTTTAAATCATCTTTACTAAT‑TCCATACACCAA‑TAGGAATAGCCACCTAAAATC‑GTAATAAAAAAATCA‑ACCGCCAATAAATAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAGTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAAGTTATGAATGGAATGaaaa                                                                                                                        <  1:656569/235‑1 (MQ=255)
                                                                                                                     aaCGACAGTGACACTACCTG‑CGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGGACTACTCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAGATAAAACT‑ACAGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaa                                                                                                                         >  2:422010/1‑250 (MQ=255)
                                                                                                                            gTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcat                                                                                                                  >  1:926310/1‑250 (MQ=255)
                                                                                                                                      ccggagctttGTGCA‑AAATTGGGATCTCCTTTCCCCGTGACATAGGCGTCATCTGTACTAATCATTTTTTCAAAT‑CCCATTAACAAA‑TAGGAATAGGCACCTAAAAAC‑AAAAAAAAAAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCa                                                                                                        <  1:349930/243‑1 (MQ=255)
                                                                                                                                                      aaaTTGG‑ATCTGCATCCCCCGGGACATAGGCGTCATCTGTACAAATCATTTCTTCAAAT‑CCAATTGCCCAA‑TAGGCATAGGCACCTAAAAAC‑GAAAAAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGGAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCa                                                                                        <  1:30794/251‑1 (MQ=255)
                                                                                                                                                            g‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGtt                                                                                  >  1:115860/1‑251 (MQ=255)
                                                                                                                                                                      cccccGTGACATAGGCGTCTTTTTTACTAATCATTTCTTCAAAT‑TCCATTAACAAA‑TAGCAATAGCCACCTAAAACC‑GAAAAAAAAAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTa                                                                          <  1:295768/249‑1 (MQ=255)
                                                                                                                                                                                    gCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTgg                                                           >  1:168926/1‑251 (MQ=255)
                                                                                                                                                                                      gTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATCAAACT‑ACCGCCAATCAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAgg                                                                   >  1:399863/1‑241 (MQ=255)
                                                                                                                                                                                         aTCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGcc                                                      >  2:837505/1‑251 (MQ=255)
                                                                                                                                                                                         aTCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGCAGGcc                                                      >  2:835181/1‑251 (MQ=255)
                                                                                                                                                                                              aacaaaacATGCCTTCGAAT‑CCCAATGACCAA‑GAGGAATAGCCACGAAAAAAA‑GCAAAAAAAAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCGCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAAtt                                                 <  1:835181/244‑1 (MQ=255)
                                                                                                                                                                                                 tAATCATGTCTTCTAAT‑TCAATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATACAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGCGGCCTACTTTa                                               >  2:583890/1‑250 (MQ=255)
                                                                                                                                                                                                         tcttctAAT‑TCAATTGACCAA‑TAGGCAAAGGCACCTAAAAAC‑GCAATAAAAAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGt                                       <  1:249129/250‑1 (MQ=255)
                                                                                                                                                                                                          cttctAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GGAATAAATAACCCT‑ACCGCAAATCAAAACACGTATTTTCATCTGTTAGAATGTTTATTCTTTGCATTACTCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAACGAGCATTCATTCTTCATATCCAACACAGTGCATGTGTTGCACATTATCCACGGAACGTTGGAGGCCTAATTTAAATGAGGTGc                                     >  1:887570/1‑251 (MQ=255)
                                                                                                                                                                                                              aaaatcccATTAACCAA‑TAGGCATAGGCACCGAAAAAC‑CAAAAAAAAAAAACT‑CCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaaa                                  <  2:168926/245‑1 (MQ=255)
                                                                                                                                                                                                                                    cataggcaCCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAAc             >  1:610196/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gcaCCTAAAAAC‑GAAAAAAAAAAATCT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTgg        <  1:253413/251‑1 (MQ=255)
                                                                                                                                                                                                                                              caaaaaC‑GAAAAAAAAAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCttt    <  1:287549/248‑1 (MQ=255)
                                                                                                                                                                                                                                               agaaaccgaaataaaaaAACCT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACTCAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCtt     <  1:585423/240‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaaaC‑GCAAAAAAAAAACCT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATTAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTAGCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTg  <  1:583890/250‑1 (MQ=255)
                                                                                                                                                                                                                                                aaaaC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCttt    >  2:287549/1‑248 (MQ=255)
                                                                                                                                                                                                                                                       cAATAAATAAAACT‑ACCGCCAATAAAGCAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCt                                                     >  1:815572/1‑193 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
GCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTG  >  NC_000913/2483021‑2483504

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