Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,475,493	T→C	F19F (TTT→TTC)	yjgL →	protein YjgL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,475,493	0	T	C	98.3%	190.5 / NA	58	F19F (TTT→TTC)	yjgL	protein YjgL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/38);  minor base A (1/0);  total (20/38)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

ACCGCTGTGCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAG  >  NC_000913/4475248‑4475739                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          aCCGCTGTGCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGGTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAAc                                                                                                                                                                                                                                                       >  2:212770/1‑250 (MQ=255)   cGCTGTGCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTcaaacaaa                                                                                                                                                                                                                                                    <  1:749253/251‑1 (MQ=255)         gcgcAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTg                                                                                                                                                                                                                                              <  1:859790/251‑1 (MQ=255)          cgcAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGa                                                                                                                                                                                                                                             <  1:236585/251‑1 (MQ=255)            cAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAg                                                                                                                                                                                                                                           <  2:172870/251‑1 (MQ=255)              gTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGtt                                                                                                                                                                                                                                         <  1:968599/251‑1 (MQ=255)                gggCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGttt                                                                                                                                                                                                                                        <  1:880140/250‑1 (MQ=255)                 ggCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTa                                                                                                                                                                                                                                      <  1:497490/251‑1 (MQ=255)                  gcgcgtCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGTTGGCGTGATGATTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTaa                                                                                                                                                                                                                                     <  1:251419/251‑1 (MQ=255)                       ccAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATAcc                                                                                                                                                                                                                                <  1:946826/251‑1 (MQ=255)                        cAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGCTGGCGTGATGACTTTTCGCAGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGCTATTTTTTTGTGAGTAAACTTTGTAATGCTTAGACGTTCTTATTCCCTCAAGGAGTTCGTCATGCGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGCAGTTTTCGATACCt                                                                                                                                                                                                                               >  2:111856/1‑251 (MQ=255)                                        aGCAGTAGTGAAGATGCTGCTGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGGTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCTCTCAAGGAGTTCGTCATTAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCCAACAAAAAAGTGTAGTTTTAAATACCTGGCGTGTTGgaatga                                                                                                                                                                                                                >  2:317803/1‑250 (MQ=255)                                                 gAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTg                                                                                                                                                                                                       <  1:10185/250‑1 (MQ=255)                                                  aaGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTGTTTTGTGAGTAAAAGTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGcc                                                                                                                                                                                                     <  2:754197/251‑1 (MQ=255)                                                        gtGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑tt                                                                                                                                                                                              <  1:495590/250‑1 (MQ=255)                                                                                  aaaTAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑a                                                                                                                                                                    <  1:212770/250‑1 (MQ=255)                                                                                                  aCGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑gg                                                                                                                                                   >  1:187098/1‑250 (MQ=255)                                                                                                         gTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCAtttt                                                                                                                                           <  1:144264/251‑1 (MQ=255)                                                                                                             tttaaaaGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGGTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTag                                                                                                                                        <  2:65860/247‑1 (MQ=255)                                                                                                               gAAAAGAGCGTTTTTTGATAATTTTTTGTTAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACGCAAGGAGTTCGTCATGAGCAAAATATCAGATTAAAATTATTCTCAACACATTACATTAGCCGACAAGTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTATAGt                                                                                                                                     <  2:256159/251‑1 (MQ=255)                                                                                                                aaaaGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGACGGGGGGCAGGCTAACAGAAGG‑AATATTATTTCTCCT‑GGTGAATTTTTAGAGtt                                                                                                                                    >  2:916861/1‑251 (MQ=255)                                                                                                                    gagCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTg                                                                                                                                >  1:159509/1‑251 (MQ=255)                                                                                                                      gCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGc                                                                                                                              <  1:366512/251‑1 (MQ=255)                                                                                                                         ttttttGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATTAGCAATATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCAGTC‑TTGCCGGGGGGCAGGATAACGTGAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCaa                                                                                                                            >  2:187751/1‑250 (MQ=255)                                                                                                                           ttttGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAATGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGTGCATTTTTAGTGTTTTTGGCaa                                                                                                                            >  2:733856/1‑248 (MQ=255)                                                                                                                             ttGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCa                                                                                                                                                                                   >  1:227697/1‑193 (MQ=255)                                                                                                                             ttGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCa                                                                                                                                                                                   <  2:227697/193‑1 (MQ=255)                                                                                                                                atagttttttGTGAGTAAAAGTTGTAATGCTTAGCCGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATAttt                                                                                                                    <  2:585720/247‑1 (MQ=255)                                                                                                                                  atTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGTCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTa                                                                                                                   <  1:214113/250‑1 (MQ=255)                                                                                                                                    ttttttGTGAGTAAAATTTGTCATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGACGTTTTAACTACCTGGCGTGTTGGAATGAATGCTTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGA‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGAAAAGATATTTAccc                                                                                                                >  2:499169/1‑251 (MQ=255)                                                                                                                                     tttttGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAccc                                                                                                                <  2:115391/250‑1 (MQ=255)                                                                                                                                          gtgAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGt                                                                                                           <  1:775995/250‑1 (MQ=255)                                                                                                                                            gAGTAAAAGTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGTCAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAt                                                                                                        <  2:775728/251‑1 (MQ=255)                                                                                                                                                aaaattgttaagtcctAGACGTTCTTAATCACTCAAGGATTTTGTCATGATCAAAATATCATATTTAAATTATTCTCAACACATTACATTAGCCTACAAGTTCAAACAAAAAAGTGAAGTTTTAAAGACCTGGCGTGGTGGAATGTATGATTTTGCCCGTA‑TTGCCGTGGGTCATGATAACAGAAGG‑AATATTCTTTCGCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTCCTCTGGGTTAtgttg                                                                                                    <  2:49585/236‑3 (MQ=255)                                                                                                                                                   aTTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGAtt                                                                                                 <  2:274638/251‑1 (MQ=255)                                                                                                                                                   aTTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGAtt                                                                                                 <  2:230118/251‑1 (MQ=255)                                                                                                                                                        tAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTAAATTAGCCGACAATTTCAAACAAAAAAGAGAAGTTTTAAATACCTGGCGTGTTGGAATGA‑TTAATTTGCCCCTAATTCCCGGGGGGAAGG‑TTACAAAAGGGATTATT‑TTTTTCCTTGTGGCTTTTTTGGAGTTTTTGGCAAAAATTTTTCCCCTGGGtttgggggtttttagc                                                                                            >  1:627633/1‑237 (MQ=255)                                                                                                                                                                aGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCt                                                                                     <  2:979807/250‑1 (MQ=255)                                                                                                                                                                 gACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTcc                                                                                   <  1:436871/251‑1 (MQ=255)                                                                                                                                                                     ttCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACg                                                                               >  2:83297/1‑251 (MQ=255)                                                                                                                                                                     ttCTTATTCACTCAAGGAGTTCGACATGAGCAAAATTTCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGACTGAATGATTTTGCGCGTA‑TTGCCGGGGGGCTGGATAACCGAGGG‑CATCTTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAACCGCTCCAACg                                                                               >  2:878427/1‑251 (MQ=255)                                                                                                                                                                           ttCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCgg                                                                         >  1:307796/1‑251 (MQ=255)                                                                                                                                                                                      aGTTCGTCTTGCGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTAAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAGGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCATGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAatttg                                                              >  2:590963/1‑248 (MQ=255)                                                                                                                                                                                          cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTCCCAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg                                                          >  1:503347/1‑251 (MQ=255)                                                                                                                                                                                             cATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCt                                                        <  2:405098/250‑1 (MQ=255)                                                                                                                                                                                             cATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCt                                                        <  1:766399/250‑1 (MQ=255)                                                                                                                                                                                                gAGCAAAATAACAGATGTAAATTATTCGCTACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGTGGGCAGGATAACAGAAGG‑AATAATCTTTCTCCG‑GTAGCATTTTTAGAGTTTTTTGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCat                                                     <  2:756955/250‑1 (MQ=255)                                                                                                                                                                                                       atatCAGATTTAAATTAATCGCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATTACTGGCGTGTTGGAATGAATGATTGTTCCCGTA‑TTGCCGGTGGGCAGGATAACATAAGG‑AATATTCTTTCTCCT‑GGAGTATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTTGATTTTAGCAAACGCTCCAACGAATCGGGTAGAAATATGATGGCTCATATTAAGtt                                             <  2:71829/251‑2 (MQ=255)                                                                                                                                                                                                          tCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTc                                          >  1:565918/1‑251 (MQ=255)                                                                                                                                                                                                                ttAAATTATTCTCAACACATTACATTAGCCGGAAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCtt                                     <  2:267561/250‑1 (MQ=255)                                                                                                                                                                                                                    attattCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAGATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCt                                <  1:111856/251‑1 (MQ=255)                                                                                                                                                                                                                           ccAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATGCCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATATTAGCAAACGCTCCAACGAAGCGGGTAGAAATAGGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATa                         <  2:259034/250‑1 (MQ=255)                                                                                                                                                                                                                            cAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGAt                          <  1:727523/249‑1 (MQ=255)                                                                                                                                                                                                                                acattacatgAGCCGACAATTTCACACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGAAAACAGAAGG‑AATATTCTTTCTCCT‑GGTGCATTTTTAGAGTTTTTGGCAAGGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCAAACGACGGGGGAAGAAATATGATGGCTCATATTAAGTCCTAATTTTATTCTAAAGATACTAAt                    >  2:989329/1‑251 (MQ=255)                                                                                                                                                                                                                                       attaGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATg             <  1:317803/251‑1 (MQ=255)                                                                                                                                                                                                                                          aGCCGACAATTTCAAACAAAAAAGTGACGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTTCCGGGGGGCAGGTTAACAGAAGG‑AATATTCTTTCGCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGaaa          <  2:654754/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aaTTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAACTGa    >  1:493653/1‑250 (MQ=255)                                                                                                                                                                                                                                                 aaTTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGa    >  1:706667/1‑250 (MQ=255)                                                                                                                                                                                                                                                  aTTTCAAACAAAAAAGTGAAGTTGTAAAGACCGGGCGTGTTGTAATGAATGATTTTGCCCGTA‑TTGCCGTGGGGCAGGCTAACAGAAGT‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGTCAATGAAAAAATGaat  <  2:19798/251‑2 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          ACCGCTGTGCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTA‑TTGCCGGGGGGCAGGATAACAGAAGG‑AATATTCTTTCTCCT‑GGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAG  >  NC_000913/4475248‑4475739

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