Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,483,254	0	G	A	35.6%	67.9 / 18.3	45	S29L (TCA→TTA)	emrK	tripartite efflux pump membrane fusion protein EmrK
Reads supporting (aligned to +/- strand):  ref base G (24/5);  new base A (0/16);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 3.15e-08
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‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGG  >  NC_000913/2483021‑2483498                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       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)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       GCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGT‑CACCTTGTCGAACGTAGTTCGTATCTTTA‑TGATTAACGACAGTGACACTACCTG‑AGACTTGTGCAGAAATTGG‑ATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAA‑TAGGCATAGGCACCTGAAAAC‑GCAATAAATAAAACT‑ACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGG  >  NC_000913/2483021‑2483498

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