Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,289,078	C→A	21.9%	Q89K (CAG→AAG)	yraJ →	putative outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,289,078	0	C	A	21.9%	145.5 / 7.2	67	Q89K (CAG→AAG)	yraJ	putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (40/11);  new base A (2/13);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 8.32e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCACTGTACCGTGATGCAA  >  NC_000913/3288851‑3289315                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 cGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAAAATGCAC‑AACTACCAGGAACTATTCAGGTTGC‑TATCTGGCTGAATAAAAAAAAGGTTTCAAAAAAAAAAA‑tt                                                                                                                                                                                                                                      >  2:52819/1‑239 (MQ=255)           ggCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑Acc                                                                                                                                                                                                                           >  1:777798/1‑238 (MQ=255)           ggCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GGGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGT‑TATCTGGCTGAATAAAAAAAAGGGTTCACAGAAAAAAA‑TTAAA‑ATT‑ACCg                                                                                                                                                                                                                          >  2:944467/1‑239 (MQ=255)                 catcatcaAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGCATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCG‑TTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGCTGA‑TATCTGGCTGAATAAAAAAAAGGTTTACCAGAAAAAAA‑TTAAC‑TTT‑ACTGCCAATga                                                                                                                                                                                                                   >  1:365813/1‑238 (MQ=255)                                  tGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCa                                                                                                                                                                                                    >  1:337395/1‑238 (MQ=255)                                  tGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGCTGA‑TATCCGGCTGAATAAAACGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCAAATGCCGAGCAACTTTTGCa                                                                                                                                                                                                    >  2:523858/1‑238 (MQ=255)                                          gTGCTCACTGGCAGCGGCATTGTTTGCACTCCCGGAAACGCCGAAGGGTATTATTTCGACCCCATTATGCTGGAAACCAAAAAAA‑GTGGTATGA‑AAACAACGGATCTGTCGCGTTTTTCAAAAAAATATGCCA‑AACTACCAGGAATTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTACAAAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGtt                                                                                                                                                                                           <  1:429841/239‑1 (MQ=255)                                                  tGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGCAGAGTATTATTTCGACCCCATTATGCTGGAAACCACCAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAg                                                                                                                                                                                   >  1:97708/1‑239 (MQ=255)                                                     cAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCCATGCAGATCAACTTCTGCAGCCCCAATTTACgggaaaaa                                                                                                                                                                                >  2:962291/1‑233 (MQ=255)                                                          gCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGAAACTTCTGCAGCCCAAGCTTACGGTAGAAAAACt                                                                                                                                                                            >  1:338782/1‑238 (MQ=255)                                                                 ttGCACTACCTCTAACGACGAAAAGTATTATTTCGCCCCCATTATGCTGGAACCCAAAAAAA‑GTGGTATGC‑AAACAACCGTTCTGGCGCGTTTTTCAAAAAAAAATGCAC‑AACTACAAGGACTTTATCAGGTTGA‑TATCTGGCTGAAAAAAAAGAAGTTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAg                                                                                                                                                                     <  1:375817/238‑1 (MQ=255)                                                                  tccACTCCCGCTAACGCAGAAGAGTTTTATTTCGCCCCAATTATGCTGAAACCCAAAAAAA‑GGGGTATGC‑AAACAACCGATCTGTCGGGTTTTTCAAAAAAATAGGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCt                                                                                                                                                                   <  1:48949/237‑1 (MQ=255)                                                                  tGCCCTCCCGCAAACGCCAAAGAGTTTTATTCCGCCCCAATTAGGCTGAAACCCAAAAAAA‑GTGGAATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAGAATATGCAA‑ACCTACCAGAAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCt                                                                                                                                                                   <  1:155895/239‑1 (MQ=255)                                                                            cTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAg                                                                                                                                                         >  1:992353/1‑239 (MQ=255)                                                                            cTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGAAACTTCTGAAGCCACAGTTTACGGTAGAACAAATACGGCAGCTGGGATTTAAg                                                                                                                                                         >  1:44795/1‑239 (MQ=255)                                                                               aaagagAAGGGTATTATTCCCACCCCATTATGCT‑GAAAACACAAAAAGGTGGTAGGC‑AAACAACCGATTTGTCGCGTTTTAAAAAAAAATAGGAAC‑AACTACAAGGAATTTAGCAGTTGGA‑TATTTGGCTGAAAAAAAAGAATGTTCAAAAAAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGt                                                                                                                                                       <  1:293802/233‑1 (MQ=255)                                                                                 gCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTAAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGACCACCTTCTGCAGCCACAGTTTACGGTAAAAAAACAACGTGAACTGGGAATAAAGGGGGa                                                                                                                                                    >  2:490410/1‑239 (MQ=255)                                                                                            ttattTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTCCGTGAGCTGGGAATTAAGGTGGTTGAAATCCCgg                                                                                                                                         >  1:370585/1‑239 (MQ=255)                                                                                              agtcgcaCCCAATTACGCTGAAAACCAAAAAAA‑GGGGTATGA‑AAACAACCGTTCTGGCGGGTTTTTAAAAAAAATAGGCAA‑ACCTACCAGGAACTTATCAGGTTGA‑TATCTGGCTAAATAAAAAGAAGGTTT‑ACAAAAAAAAATTTACT‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGgcg                                                                                                                                       <  1:797066/233‑1 (MQ=255)                                                                                                 tCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGACAAACAACGTGAGCTGGGTATTAAGGGGGATGAAAACCCGGCgctg                                                                                                                                    >  2:467896/1‑239 (MQ=255)                                                                                                  cGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAAAAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCgctg                                                                                                                                    >  1:418738/1‑238 (MQ=255)                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTTTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTAACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGTTCAAACAAACACGTGAGCTGGGGCTTAAGGCGGATGAAATccctgtcccttaa                                                                                                                                 >  2:293802/1‑229 (MQ=255)                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCgctg                                                                                                                                    >  1:713532/1‑236 (MQ=255)                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCcctggct                                                                                                                                 >  1:646135/1‑239 (MQ=255)                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGAGAGCTGGGGATTAAGGGGGAGGAAAACCGGGCgcgggcg                                                                                                                                 >  2:997828/1‑238 (MQ=255)                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTAACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCCCAGTTTACGGTGGAACAACTACGTGAGTTGGGTATTAAGGGGGAGGAAAACCCgagacgggcg                                                                                                                                 >  2:682070/1‑230 (MQ=255)                                                                                                       ccATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGACTTAAGGGGGATGAAATCCCGGGGCTGGCTGaa                                                                                                                              >  2:509067/1‑239 (MQ=255)                                                                                                        cATTATGCTGGAAACCAAAAAAA‑GTGTTATGC‑AAACAACCGATCGGTCGCTTTTTTAAAAAAAATATGCAC‑AACTACCAGGAATTTACAAGGTTGA‑TATCTGGCTGAAAAAAAAGAAGGTTTCACAAAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGaaa                                                                                                                             <  1:509067/239‑1 (MQ=255)                                                                                                             ttgcggAAACCCAAAAAA‑GGGGGAAGGAAAACAACCGATCGTTCGGTTTTTTCAAAAAAATATGACA‑AACTCCCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGTTTTCAAAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGa                                                                                                                         <  1:52819/235‑1 (MQ=255)                                                                                                                 ggAAACCAAAAAAA‑GGGGAAACA‑AAAAAACCGATCTGTCGCTTTTTTAAAAAAAAATGGCAA‑AACTACAAGAAACTTAAAGGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCAAAAAAAAAAA‑TTAAA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACg                                                                                                                    <  1:944467/239‑1 (MQ=255)                                                                                                                  gAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGGAATCCCGGCGCTGGCTGAAAAAGATGACGa                                                                                                                   >  1:72452/1‑239 (MQ=255)                                                                                                                  gAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGa                                                                                                                   >  1:758648/1‑239 (MQ=255)                                                                                                                       cacaAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGAATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCg                                                                                                              >  1:627871/1‑239 (MQ=255)                                                                                                                         aaaaat‑GTGGTAGGA‑AAACAACCGATTTGTCGCTTTTTTCAAAAAAATATGCCC‑AACTACCAGGAACTTATCAGGTGTA‑TATCTGGATAAAAAAAAAAAAGGTTTCAAAGAAAAAAA‑TTACT‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTg                                                                                                            <  1:523858/238‑1 (MQ=255)                                                                                                                              a‑GTTGTATGC‑AAACCCCCTGTCTTGCGCTTTTTTAAAAAAAAGAGTCCA‑AACCACCAGCAAATTAACAGGTTAA‑TTTCTGGCTAAATAAAAAGAAGGTTTCAAAGAAAAAAT‑TTACA‑TTG‑ACCGCCAATGCAGAGCAGCTTCTGCAGCCACAGGTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCaa                                                                                                       <  2:264411/239‑1 (MQ=255)                                                                                                                                          aaaCAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCCCTTGac                                                                                             >  1:399379/1‑238 (MQ=255)                                                                                                                                           aacaccCGATCTGCCGCGTTTTACAAAAAAATATGCAC‑AACTACCAGAAACTTATCAGGTTGA‑TATCTGGCTAAATAAAAAGAAGGTTTCAAAAAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAAc                                                                                            <  1:467896/239‑1 (MQ=255)                                                                                                                                                   gtttttttgcttttttAAAAAAATTAATACAAACTTCCAGAAAGTTAACAGGTTGT‑TATGTGGCTAAAAAAAAAGAAGCTTCAACAGAAAATAA‑TTACT‑TTT‑CCCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCAt                                                                                     <  2:362950/229‑1 (MQ=255)                                                                                                                                                   aTCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAAt                                                                                                                                              >  1:459428/1‑181 (MQ=255)                                                                                                                                                   aTCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCAt                                                                                     >  1:274088/1‑238 (MQ=255)                                                                                                                                                      tGTCGCGTTTTTAAAAAAAATGGGCAA‑AACTTCCAGGACCTTATCAGGTTGA‑TTTTAGGCTGAAAAAAAGGAAGGTTTCAAAAAAAAAAA‑TTACA‑TTA‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTccc                                                                                 <  1:899336/239‑1 (MQ=255)                                                                                                                                                          tataaaaaacaaaaaaTACGAAC‑AA‑TACCCAGAACTTTACTAGTTTATTATTTGGTTAAAAAAAAAAAAGG‑TTCACAAAAAAAAA‑TTTCC‑TTTAACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTa                                                                             <  2:928614/229‑1 (MQ=255)                                                                                                                                                                              cac‑aACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATT‑CCGGTACAGCTGCTGAATTTGATTTCa                                                        >  1:248469/1‑239 (MQ=255)                                                                                                                                                                                   aCTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGa                                                                   >  1:928614/1‑225 (MQ=255)                                                                                                                                                                                    cTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGAGGATGAAATCCCGGCGCTGGCTGAACAAGATGACGATAGCg                                                                                                              >  1:300445/1‑181 (MQ=255)                                                                                                                                                                                       ccAGGAATTTATCAGGTTGA‑TATCTGGCTAAAAAAAAAGAAGTTTTCAAAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAAtcatca                                                 <  1:490410/239‑1 (MQ=255)                                                                                                                                                                                              cTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGAt                                                                                                          >  1:346233/1‑175 (MQ=255)                                                                                                                                                                                              cTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACACCTGCTGAATTTAATTTCAATAATCAGCAAAtt                                          >  2:733807/1‑239 (MQ=255)                                                                                                                                                                                                tATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAACAA‑TTCCA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGCAATCCCGGCGCTGCCTGAAAATAAAGACGATAGCGTGCTCAACACGCTTGAACAAATAATTCCCGGTACAGCTGCTAAATTTGATTCAAACAATAAGCACCTTaa                                        >  2:984837/1‑239 (MQ=255)                                                                                                                                                                                                 aTCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAAt                                       >  1:839878/1‑239 (MQ=255)                                                                                                                                                                                                      gTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGAGGATGAACTCCCGGCGCTGGCTGACAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGGACAGCTGCTGAATTTATTTTCAATCATAAGCACCTACATTTGa                                   >  2:754460/1‑238 (MQ=255)                                                                                                                                                                                                         gA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGa                                                                                                                                                 >  1:986799/1‑125 (MQ=255)                                                                                                                                                                                                         gA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTg                                                                                                            >  2:184336/1‑162 (MQ=255)                                                                                                                                                                                                         gA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCACGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCAt                               >  1:287549/1‑239 (MQ=255)                                                                                                                                                                                                         gA‑TATCTGGCTGAATAAAAAGAAGGTTTCAAAGCAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACACCTACGAGAGCTGGGTATTAAGGCGGATGAAATCCCGGCGCTGACTGAAAAAGATGACGCTCACGTGCTCACATCGCTTAGACAAATCCTTCCCGGTACACCGGCTGAATTTGATTTCAACAATCAGCCACTTAATTTgcaca                                >  2:674308/1‑234 (MQ=255)                                                                                                                                                                                                         aa‑tatCTGGCTGAATAAAAAGAAGGTTCCAAAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTg                                                                                                            <  1:184336/161‑1 (MQ=255)                                                                                                                                                                                                                gtgtaaaaaaaaaaGAAGTTTTACCAAAAAAATA‑TTACT‑TTG‑CCCGCCAATGCGGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCa                         <  2:627871/230‑1 (MQ=255)                                                                                                                                                                                                                  gcaaaaaaaaaaGATGGTTTCAAAAAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCaaa                       <  2:992353/232‑1 (MQ=255)                                                                                                                                                                                                                  gCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGCTAGCGGGATCAACTCGCTTGAACAAATCATTCCCGGTACAGATGCGGAATTTAATTTCAATCATAAGCAACTTAATTTTAACATTCCCCcaa                       >  2:133989/1‑239 (MQ=255)                                                                                                                                                                                                                  gCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGGATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGCTCAACCGGCTTTAACAAATCATCCCCGGTACACCTCCTGAATTGGATTTCAATCACCACAAAATAAATATGAGAATTCCCaaaa                       >  2:546794/1‑239 (MQ=255)                                                                                                                                                                                                                   cgaaaaaaaaaaaaGGTTTTAAAAAAAAAAA‑TTTCACTTT‑CCCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCaaa                       <  2:141714/233‑1 (MQ=255)                                                                                                                                                                                                                       aaaaaaaaaaGTTTCCACAAAAAAAAT‑TTCCA‑TTT‑CCCGCCAATGCAGAGCCACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCa                  <  2:758648/237‑1 (MQ=255)                                                                                                                                                                                                                       aaaaaaaGTAGTTCTCACAAAAAAATA‑TTATT‑TTT‑CCCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACATCTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCa                  <  2:713532/237‑1 (MQ=255)                                                                                                                                                                                                                       aaaaaaaGGAGTTTTCAAAAAAAAAAA‑TTACA‑TTA‑TCCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCa                  <  1:754460/237‑1 (MQ=255)                                                                                                                                                                                                                       aTAAAAAGAAGGTTTCACAAAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCa                  <  1:733807/239‑1 (MQ=255)                                                                                                                                                                                                                               aagGTTTCA‑AAAAAAAAA‑TGTCATTTT‑ACCGCCAATGCAGAGCAACTTTGGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGATGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTTCCCAAATTGCACTGTACCg          <  1:997828/239‑1 (MQ=255)                                                                                                                                                                                                                                agGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTTAGCATTCCCCAAATTGCACTGTa             >  1:141714/1‑235 (MQ=255)                                                                                                                                                                                                                                       acaAAAAAAAA‑TTACA‑TTG‑CCCCCCAAGGCAGAGCAACTTCTGCAGCCAAAGTTTACGGTAGAACAACTACGTGCGCTGGGTATTAAGGTTGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCACTGTACCGTGATGCaa  <  2:777798/239‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 CGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACA‑TTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCACTGTACCGTGATGCAA  >  NC_000913/3288851‑3289315

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