Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	11,427	C→T	12.7%	V120V (GTG→GTA)	yaaI ←	UPF0412 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	11,427	0	C	T	12.7%	156.1 / 3.7	72	V120V (GTG→GTA)	yaaI	UPF0412 family protein
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (8/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 8.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTT  >  NC_000913/11200‑11653                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      gcgtCGGTTTAAACTTTTTTTTGCGCATCGTGACGCTGGAGAGGCTGGTTTTGCTTTTTTCATTATGGGGTAGCAAACGGGCGAAATTTGCTAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTc                                                                                                                                                                                                                                       <  1:680535/225‑1 (MQ=255)                     ttGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTGCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGt                                                                                                                                                                                                                  <  2:82155/228‑1 (MQ=255)                     ttGCGCAATAGTACGCTGGGGAGGGTAGTTTTTCCTTTTTCATTTAGGGTTAGCAATCTGGCGAAATTTGCCAACTTGTCCTCACTACAATGCTTAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGGGGCTTCTATATTCTGGCGTTAGTCGGCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGt                                                                                                                                                                                                                  <  1:691553/228‑1 (MQ=255)                                  cGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAACAATTGGGCGAAATTTGCCAACTGTTCCTCATTACAATGCCGAAGAAAATCCAGATCTGAATCATTCAGGGAATTAACATTCATTTTTTGTGGCTTATCTGTTTTGGCGCTTGTTTCCGCCGATACTTTTTTGTGTTGCCCTAACCTACGAGTTTCCCGTTTG‑ACCCGAAAAAGTGATTTTTtagaac                                                                                                                                                                                                     >  1:361268/1‑226 (MQ=255)                                   gCTGGGAAGGCTAGTTTTGCCTTTTTCATTATGGGGGAGCAATCGGGCGAAATTTGCCAACTGTTCATCACTACAATGCTGAAGAAAAGCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGGGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGc                                                                                                                                                                                                    <  1:386042/228‑1 (MQ=255)                                    cTGGAGAGGAGTGTTTTGCCTTTTTCATTAGGGGTGAGAAATCGGGCTAAATTTGCCAACTGTTCCTCACTACAATGCTTAAGAAAATCCATATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCATCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCa                                                                                                                                                                                                   <  1:496409/228‑1 (MQ=255)                                     tGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGGGTTTGTCGTGGCCCACAATTTTTCGCGTGGCCATATCCATTGGGGTTACCGCATT‑GACCGTAAAGGGGAGTTTTGTCACGCAg                                                                                                                                                                                                  >  1:533450/1‑228 (MQ=255)                                     tGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGGGAGCAATCGGGCGAAATTTGCCAACTTTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTTTATTCTGGCATTAGTCGTCCCCGATAATTTTCAGCGGTGCTATATTCGATGAGTTTCCCCTATT‑ACCCCAAAAGGTGATTTTTGGAACGc                                                                                                                                                                                                    >  2:684568/1‑226 (MQ=255)                                       gagaGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTTACATTCATTTTTTGTGGCTTCTATATTCTGGCGATAGTTGACCCCGATAATTTTCAGCGGGTCCCATTCCCATTAGTTTACCGTTTG‑ACCCGGACAggcgtgtttttggtccc                                                                                                                                                                                                    >  2:309087/1‑209 (MQ=255)                                       gagaGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTCTATTCTGGCGTTAGACGTCCACCTTAATTTTCAGCTTGGGCTTATCTGTTTAGTTTCCCGTGTA‑CCACACAAAGGGGAttttttaaagccaacg                                                                                                                                                                                                >  1:227818/1‑217 (MQ=255)                                       gagaGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAACTTTCAGCGTGGCCCTATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGGGATTTTTTAGACGcgccg                                                                                                                                                                                                >  2:315623/1‑224 (MQ=255)                                       gagaGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGATCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTTTGTCGTTAGTAGTCACCGTTAATTTTTACGGTGCGCCTATCCCGTGGGTT‑ACCCTGTGTACCCCAAAAGGTGATTTTTTTGAGGCAgag                                                                                                                                                                                                >  1:622370/1‑226 (MQ=255)                                            gCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTa                                                                                                                                                                                            <  2:31109/227‑1 (MQ=255)                                               aGTTTTGCATTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTTTATTCTGGCCATGGTCGTCCCCGATAATTTTTAGCGGTGCCCTAACCGATGAGTTCCCCTTTTG‑ACCCGCAAAAGGGATTTTTGTGAGGCGCGGTTTATTgc                                                                                                                                                                                        >  1:404238/1‑227 (MQ=255)                                                       cTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTGGTGGCTTCTATATTCTGGCGTTAGGCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATc                                                                                                                                                                                 <  2:247025/227‑1 (MQ=255)                                                           ttCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGTTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCg                                                                                                                                                                                <  2:697388/224‑1 (MQ=255)                                                            tCATTTGGGGTGAGCAAAGGGGGTAAATTTGCCATATTGTCCTCACTAAAATGCTGAATAAAAGCCAGATCTTAATCATGCAGGTAATTCACATTCATTTTTTGTGGCTTCTATTTTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGATTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGt                                                                                                                                                                            <  1:159781/227‑1 (MQ=255)                                                                tATGGGTGCGAAATCGGGCGAAATTTCCAAACTGTTCCTAACTACAATGCGTAAGAAAAGCCAGATCTAAATTATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTATTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATg                                                                                                                                                                       <  1:359461/228‑1 (MQ=255)                                                                     gTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCTGCGTTGCCATATCCGATGCGTTCACCCTATG‑ACCCGAAAAGGTGCTTTTTTAGTCGCAGCGTTTATTGTCCTTTTCGCTGTTAATGtatat                                                                                                                                                                  >  2:680818/1‑224 (MQ=255)                                                                     gTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGAt                                                                                                                                                                  <  2:361268/228‑1 (MQ=255)                                                                        aGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTTTATTCTGGCGTTAATCCTCGCCGACTATTTTCAGCGTTGCCATTTCCGATGATTTCCCCGTTTG‑ATCGGCTAAGGTCATTTTTGAGCCCCAGCGTTTATTGTCGTTATATCTTTTaaatcttaccca                                                                                                                                                               >  2:333136/1‑218 (MQ=255)                                                                        aGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGGGTGGCCCTATCCTATGAGTTCACCGTCTG‑ACCCGAAAAGGGGATTTTTTGTACGCAGCGGTTATTTTCGTTATCCGTGTTAATTTTTATc                                                                                                                                                                 >  1:617628/1‑226 (MQ=255)                                                                         gCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTg                                                                                                                                                                    <  2:617628/222‑1 (MQ=255)                                                                                  cGAAATTTGCCAACTGTTCCTAATTACAATGCTGAATGAAATCCAGATCTGAATCATTCAGGTAATTAACCTTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGt                                                                                                                                                     <  2:97660/227‑1 (MQ=255)                                                                                            tgcgggTTCGTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTACCATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCACCGATACTTTTTAGCGTGTCCATATCCGCTGATTTCACCGGATG‑ACCCGAAATGGTGATTTTTGGGTCGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGTTTCTATCC‑GTGGGTTTC‑Ccttttt                                                                                                                                         >  1:271563/6‑228 (MQ=255)                                                                                                     cTCACTACAATGTTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGACGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTg                                                                                                                                <  2:636205/228‑1 (MQ=255)                                                                                                          tACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTTGTCGCCGTTAATTTTCAGGCTGGCCCTATACGATGAGATTACCGTCTG‑CCACGAAAAGGTGGTTTTTGGGCAGCCGGGTTTATTGCAGTTATCGCAGTTAATGTGTATGCAGGAA‑CGGGTTTTC‑CCTTCTTTTATTTCTGAAgg                                                                                                                            >  2:535056/1‑227 (MQ=255)                                                                                                          tACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGGGGTTTTTGAGACGCAGCGGTTATTTGCGTTATCGCTGTTAATGTTTGTCCATGCA‑GTGGTTTGC‑CCTTCTTTTATTTTTGAAgg                                                                                                                            >  2:458543/1‑227 (MQ=255)                                                                                                           aaaaTGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGaa                                                                                                                          <  2:533450/226‑1 (MQ=255)                                                                                                           aCAATTCTTAATAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTGTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGaa                                                                                                                          <  1:458543/228‑1 (MQ=255)                                                                                                               tGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGTGGTTAGTCGTTGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAATGTGATTTTTGTGACGCATCGGTTATTTCCGTTTTCGCTGTTAAGGTTGATCCAGTCA‑GGGGTTTTCTCCTTTTTTTTTTTTTTAA‑GAatata                                                                                                                      >  1:63929/1‑227 (MQ=255)                                                                                                               tGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGCTAGACGTCGCCGATAATTTTCAGCGTGGCCATATCCGTTGAGTTTACCCTATG‑ACCCGAAAAGTTGATTTTTTAAACCGCGCGTTTATTGTCGTTATCCCTGTTAATGTTTATCCAGTCA‑TGGGTTTGC‑CCTTCTTTTTTTTTTTAAAGAatat                                                                                                                       >  1:165465/1‑227 (MQ=255)                                                                                                               ggCTGAAGAAAATCCAGATCTGAATTATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTAAATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATAtt                                                                                                                      <  2:165465/227‑1 (MQ=255)                                                                                                                    aagaaAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTTACCGTATG‑AACCGAAAAGGGGATTTTTGAGACGCAGGGTTTATTGTTGTTATCGCTGTTAATGTTGATCCAGTCA‑ATGGTTTGC‑CCTTCTTTTATTTTTTAAGGAATATTTAgg                                                                                                                  >  2:302447/1‑227 (MQ=255)                                                                                                                    aagaaAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGGTTGC‑CCTTCTTTTATTTTTGAAGGAATATTCAgg                                                                                                                  >  2:529089/1‑227 (MQ=255)                                                                                                                     taaaaaTCCATATCTGAATCGTTCAGGTTATTAACCTTCTGTTTTTTTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAAATTTCAGCGTTGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCATTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGct                                                                                                                <  1:434342/226‑1 (MQ=255)                                                                                                                         aaTCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGGGGCTTCTATATGCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGa                                                                                                            <  2:404238/228‑1 (MQ=255)                                                                                                                                     aaTCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATCTTCTGCCGTTAGTCTTCGCCGATATTTTTCAGCGTGGCCCTATCCGATGAGTTCACCGTATG‑ACCCGACAAGGTGCTTTTTGGGACGCAGCGTTTTTTGTCGTTATACGTGTTAATGTTGTTCCATTCG‑GTGGGTTGC‑CCTTCTTTTTTTTTTTAAGGAAATTTTAGGcccggcacttttccttt                                                                                                 >  1:348040/1‑211 (MQ=255)                                                                                                                                      tttttttttgTATTTACCATTTATTTTTTGGTGCTTCTATATGTTGGCGTTAGTGGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACggg                                                                                               <  1:302447/219‑1 (MQ=255)                                                                                                                                      attcttCAGGTTATTAAATTTCATTTTTTGGGGCTTCTATTTGCTGGCGTTGGTTGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTTTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACggg                                                                                               <  1:450496/224‑1 (MQ=255)                                                                                                                                      aTCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACggg                                                                                               <  2:22603/228‑1 (MQ=255)                                                                                                                                      aTCATGCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATAtt                                                                                                                      <  2:348040/205‑1 (MQ=255)                                                                                                                                           tatggTTATTTATATTTATTTTTTGTGGCTTGTATATTCTGGCGTTAGTCGTTGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCt                                                                                          <  1:309087/225‑1 (MQ=255)                                                                                                                                            cAGGTAATTGACATTCTTTTTTTGGGGCTTGTATATTCTGGCGTTAGTCGTTGCCGATAATTGTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCt                                                                                          <  1:4465/227‑1 (MQ=255)                                                                                                                                                  aTTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGGTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTTAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGATCTGACTGGCGCTCCGGTCGGCTTTGgaa                                                                                    >  2:59030/1‑227 (MQ=255)                                                                                                                                                  aTTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGCTAGTCGTCACCGATAATTTTTAGCGTGGCCATGTCCGCTGAGTTCAACGTATG‑ACACGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTaac                                                                                                                                                                        >  2:132516/1‑144 (MQ=255)                                                                                                                                                   ttATCAGTCATTTTTGGGGGCTTCTATATTCTGGCGGTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCTTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGGTTATTGTCGTTAGCGCTGTTAATg                                                                                                                                                                       <  1:132516/145‑1 (MQ=255)                                                                                                                                                   ttAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATAGCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATa                                                                                  <  2:365872/228‑1 (MQ=255)                                                                                                                                                   ttAACATTCATTTTTTGTGGCTTCTATATTATGGCGTTAGTCGTCGCCGATAATTTTCCGCGTGGCCATATACGATGAGTTCACCGTATG‑ACCCCAAAAGGTGATTTTTTAGACGCACCGTTTATTGTCGGTATCGCTGTTAATGTTTATCCA‑TCACGAGGTTTGC‑CCTTTTTTTATTTCTTAAGGAATATTCAGGCTCTGACTGGGCCCACCGGGGGGTTTGaaaa                                                                                   >  1:571302/1‑226 (MQ=255)                                                                                                                                                         ttcattTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTTAGCGTGGCCATATCCGATGAGTTCACCGTTTG‑ACCCGAAAAGGTGATTTTTGAGACGCAGGGTTTATTGTCGttta                                                                                                                                                                                  >  1:234427/1‑126 (MQ=255)                                                                                                                                                         ttcattTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTTAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAAGTGATTTTTGAGACGCAGCGTTTTTTTGCGTTATCGCTGTTAATGTTTATCCAGGCA‑GGGGGTTGC‑CCTTTTTTTTTTTCTTAAGGAATATTCAGGCCCTGGCTGGGGCAACGGGCGGCTTTGGAATaaaaaaa                                                                            >  1:15718/1‑224 (MQ=255)                                                                                                                                                         ttcattTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTAACCGTCTG‑ACCCGGAAAGGGGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGATTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTATGGGCGGCTTTTGAATAAACCg                                                                             >  2:435918/1‑227 (MQ=255)                                                                                                                                                         ttcattTTTTGTGGCTTCTATATTCTGGCGGTAGTCGTCGCCGATAATTTTCAGCGTTGCCATATCCGATGAGTTCACCGTATG‑GACCGAAAAGGCGATTTTTGAGACGCAGTGTTTACTTTCTTTATCGCTGTTAATGTTGATACAGTCA‑GTGGTTTGC‑CCTTCTTTTCTTTCTGAAGGTATATTCAGGATCTGACTGGCGCTCTGGCCGGGTTTGCAATAAAccct                                                                            >  2:405041/1‑226 (MQ=255)                                                                                                                                                          tcattTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATc                                                                                                                                                                                 <  2:234427/128‑1 (MQ=255)                                                                                                                                                            attTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATg                                                                          <  2:309736/227‑1 (MQ=255)                                                                                                                                                                                tCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCGACGGGCGGCTTTGAAATAAACCGATGCCCCGCTTAACTGTAAAACGc                                                     >  1:13579/1‑228 (MQ=255)                                                                                                                                                                                  tGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCa                                                   >  2:179276/1‑228 (MQ=255)                                                                                                                                                                                  gtgggtttttgctcgCCGATAATTTTCAGGTTTTGCATATCCGATGAGTTCACCGTATG‑ATCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTTTCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGATTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGcc                                                    <  1:684568/215‑1 (MQ=255)                                                                                                                                                                                    gCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCGGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGACTATTCAGCCTCTGACTTGCGCTAGGGGCGGCTTTGAACTAAACCGATGCACCGCTTACCTGTAAATCGCCATg                                                 >  2:243195/1‑228 (MQ=255)                                                                                                                                                                                    gCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACgg                                                                                                >  2:294999/1‑181 (MQ=255)                                                                                                                                                                                    gCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATg                                                 >  2:663773/1‑228 (MQ=255)                                                                                                                                                                                         aGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTTATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGTATTAAACCGTTTCACCCCTTAACTGGAAATCGCCCGGGTCgg                                            >  1:302677/1‑228 (MQ=255)                                                                                                                                                                                         aGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑TGGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGTCTTTGAATTAAACAGATGGACCGCTTACCGTGAAATCGCCATGTTCgc                                            >  2:699948/1‑227 (MQ=255)                                                                                                                                                                                         aGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACAGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CTTTCTTTTTTTTCTGAAGGAATTTTCAGGCCCTGCCTGGCGTTCCGGGCGGCTTTTAAATAAACCGATGCACCCGCTAACTGTAAATCTCCAAggctgg                                            >  1:400206/1‑224 (MQ=255)                                                                                                                                                                                              tcgCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGAGTTTATTGTCGTTATCGCAGTTAATGATGTTCCAGACA‑GTGGTTTTC‑CCTTCTTTTATTTCTGAAGGATTATTCGGGATCTTCATGTCGCTACCGGCGGCTTTTAAATAAACAGTGTCACACCTTAACTTTAAACCGCCGTgggcagaaaaa                                       >  1:253185/1‑219 (MQ=255)                                                                                                                                                                                                      taTTTTCAGTGTTGCCATATCCGGTGATTTCCCCGTATG‑ACCCGAAAAGGTGATTTGTGAGACGGAGCGTTGATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTAt                               <  1:59030/227‑1 (MQ=255)                                                                                                                                                                                                                tgtgtcATATCCTATGGGTTTACTGTGTG‑ACCCGAAAAGGTTATTTTTGAGCCGCTGGGTTTTTTGTCGTTTTCGCTGTTAATGTGGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTcaca                     <  1:315623/223‑1 (MQ=255)                                                                                                                                                                                                                gTGGACATATCCGATGATTTCACGTTATG‑ACCCGAAAAGGTGATTTTTGAGACGGAGGGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTcaca                     <  1:680818/228‑1 (MQ=255)                                                                                                                                                                                                                      atatCCGATGAGTTCACCGTATG‑ACCCTAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGa               <  2:400206/228‑1 (MQ=255)                                                                                                                                                                                                                        ttCCGATGGGTTTACCGTTTG‑TACCGAAAAGGTTATTTTTGAGACGGAGGGTTTATTTTCGTTAACTCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACa             <  1:179276/227‑1 (MQ=255)                                                                                                                                                                                                                             aTGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAAttggccaacca         >  2:417452/1‑217 (MQ=255)                                                                                                                                                                                                                              tGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCTCGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCCATGCACCGCTGAACTGTAAATCGCCATGGTCGGCAGAGGGTTGTATGCGTTTCACAATACGCCAACCAg        >  2:435896/1‑227 (MQ=255)                                                                                                                                                                                                                                  ttCACCGCATG‑ACCCGAAAAGGAGATTTTTGAGACGCAGAGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CATTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTATGAAATAAACCGATGCACCGCTTACCTGTAAATCGCCATGGTAGGCAGAGAGTTGTATGCGATTCACAATTCGACATACATGAAgc   >  2:703626/1‑227 (MQ=255)                                                                                                                                                                                                                                    cACCGTATG‑ACCCGAAAAGGAGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAGCCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGTCAAACAGGAAGtt  >  2:446883/1‑227 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      CCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATG‑ACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCA‑GTGGTTTGC‑CCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTT  >  NC_000913/11200‑11653

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