Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,917	T→C	I105I (ATT→ATC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,917	0	T	C	95.6%	124.6 / ‑5.4	46	I105I (ATT→ATC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (18/25);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGA  >  NC_000913/4619671‑4620150                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            ttgttaCGTCGGGGCTGACACCCTGGTTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGcc                                                                                                                                                                                                                                          <  2:581922/245‑1 (MQ=255)                 gACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGt                                                                                                                                                                                                                          <  1:317203/250‑1 (MQ=255)                             cATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGata                                                                                                                                                                                                           <  1:593645/251‑1 (MQ=255)                                atCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGAAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAATGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATTTGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATAttt                                                                                                                                                                                                        <  1:446623/251‑1 (MQ=255)                                       aaGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGAt                                                                                                                                                                                                 >  1:662486/1‑251 (MQ=255)                                            ggtgCAAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAAGTTGCCAAATCTGCCCCGGCTGGGGCTGGCGAAAGCACACGAAGGTTCGACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTGCTCCGATCACGa                                                                                                                                                                                            <  2:545350/250‑1 (MQ=255)                                                  aaaGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCCCTGGGAAATCGCCGGGGTCCCGGGTCTTTTTT‑TG‑GGGGGGTGTTTCTCCGATCAAGAAAAc                                                                                                                                                                                        >  2:36325/1‑249 (MQ=255)                                                        gcaGCTGATAACGGTCGTAAATGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGGAGTTATCGGCGCGTACGCCTGGGCG‑ACAGCAATGTCAACCGGGAAAGAAAACCCGGCTGGTCACTGGGAAAACGCCGGGGGCCAGGTTCTTGTTG‑AGGTGGGG‑TAATTCTCCGAACAAGCAAACAACTTccc                                                                                                                                                                                >  2:351673/3‑251 (MQ=255)                                                         acgCTGATAACGGCAGCAAAGGCCCGCTCAATCTGACAAATCTGACCCGTCTGGGGCTGTCGAAAGCATACGTAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCTGTAAAGATCCCCCGTCTGGTTACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TTGGGATATTTCTTCGATCTCGAAAACAGCTTCCCg                                                                                                                                                                               <  2:654888/249‑1 (MQ=255)                                                         aaGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCg                                                                                                                                                                               <  1:234032/251‑1 (MQ=255)                                                            cggATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAAGCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGAt                                                                                                                                                                                                 <  2:387453/228‑1 (MQ=255)                                                            cTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGGCGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGAt                                                                                                                                                                                                 >  1:387453/1‑230 (MQ=255)                                                                        gTAAAGACCCGCTCAATCTGCCAAATCTGAACCGGCTGGGGCTGGGGAAAGAACACGTAGGTGCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTAGCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGTCTTTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAATAGCTTCTGGAt                                                                                                                                                                <  1:407674/251‑1 (MQ=255)                                                                            aGGCCCGCTCAATCTGCCAAATCTGACCAGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTAACGGTTCTGGTTG‑AG‑TGGGGATATTTCTCCGATAACGAAAACAGCTTCCCGCAATAGCTGCTTGATaaa                                                                                                                                                             >  1:587112/1‑250 (MQ=255)                                                                              gCCCGCTCAATCTGCCAAAGCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCAGCCGGTAAAGATACCCCGTCTGGTCACTGGGAAAGCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTg                                                                                                                                                          <  1:458887/251‑1 (MQ=255)                                                                                     cAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGGTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAc                                                                                                                                                   <  2:541935/251‑1 (MQ=255)                                                                                            ccAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGGAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTccc                                                                                                                                                                                                                                  <  1:565989/167‑1 (MQ=255)                                                                                               aaTCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTTATTCCGGCGGGAATGGACGGCAACGCTGAAGTTAGCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGTTAAAGATACCCCGTCTGGTCACTGGGAAATCTCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCt                                                                                                                                         <  2:638132/251‑1 (MQ=255)                                                                                               aaTCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATAGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGGCTGGTCACTGGGAAATTGCCGGGGTCCCGGGTCTGTTTGGAG‑‑GGGGATATTTCTCCGATCACGAGAACAGCTTCCCGCAAGAGCTGCTGGCTAAACTGGGCGAACGCGCTAAtat                                                                                                                                         >  1:307361/1‑249 (MQ=255)                                                                                                 tCTGACCCGTATGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GCAATGTCATCCGGTAAAGATACCCCGTCTGGTCACCGGGAAATTGACGGGTTCCCGGTTCTGCTTG‑TG‑TGGGGTTATTTATCCGATCACGAAAACAGCTTCCCGCCAGAGCCTCTGGATAAACTgtgcgaaccacccaccctcc                                                                                                                                       >  2:626593/1‑232 (MQ=255)                                                                                                                gCTGGAGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCCTACGCATGGGCGCAC‑GAAATGTCATACGGTAAAAATAACCCGTCTGGTCACTGGGAAATAGCCGGTGTCCCGGTTCAGGTTG‑AG‑TGGGGCTAGTTCGCCGATCACGAAAACAGCTTCCCGCGAGGGCTGCTGGATGAACTGGTCGAACGCGCAAATCTGCCGGGTTAACTCggg                                                                                                                         >  2:26666/1‑249 (MQ=255)                                                                                                                    gggAAAGCACCCGAAGGTTCGTCCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCTC‑GAAATGGCATCCGGTAAAGATACCCCGGCTGGTCACTGGGAAATCGCCGGTTTTCCGGTTCTGTTTG‑AG‑TTGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTg                                                                                                                    <  2:580551/249‑1 (MQ=255)                                                                                                                        aaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgcg                                                                                                                                                <  1:325299/219‑1 (MQ=255)                                                                                                                            aaccggagggtCTACCGGTTTCATTCGGGGGGGAAGGGACGGCAACGCTGAAGTTATCGGCGCGGACGCGTGGGCGCAC‑GAAATGTCATCCTGGAAAGATACCCCATCTGGTCACTGGGAAACCGCCGGGGTCCCGGTTCGGTTTG‑AG‑TGGGGATATTTCGCCGATCCCGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACGGGGCGAACGCGCTAAGCTGCCGGTTTGCCTCGGTACCTGCCTCTCtt                                                                                                            <  2:338733/241‑1 (MQ=255)                                                                                                                              acGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTcc                                                                                                          <  1:661680/251‑1 (MQ=255)                                                                                                                              acGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTcc                                                                                                          >  1:316872/1‑251 (MQ=255)                                                                                                                                   ggTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GGCATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGACGGTGTCCCGGTTCTGTTTG‑GG‑TGGGGATATTTCTCCGATCACGAAAACCGCTTCCCGCAAGAGCTGCTGGCTAAAATGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGCAAATGCCACTCTTCCGGTAc                                                                                                     >  1:638746/1‑251 (MQ=255)                                                                                                                                                           ggcaTGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATATCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCCTTCTGGATCAGCTGGGCgg                                                                             >  1:265870/4‑250 (MQ=255)                                                                                                                                                             aaTGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAAGGTCATCCGGTAAAGATACCCCGTCTGTTCACTGGGAAATCGCCGGTGTCCCTGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAag                                                                           <  1:450264/251‑1 (MQ=255)                                                                                                                                                                 gACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCAATGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑GG‑TGGGGGTATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcc                                                                       >  1:292921/1‑250 (MQ=255)                                                                                                                                                                    ggCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACAt                                                                    <  1:213379/251‑1 (MQ=255)                                                                                                                                                                          gCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GGAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTACTCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGCAGAGCACATGaaacc                                                              >  2:130413/1‑249 (MQ=255)                                                                                                                                                                                   aTCGGCGCGTACGCATGGGCGCCC‑GAACTGTCCTCCGGTAAAGATACCCCGGCTGGGAACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGACGAACGCGCTAATATGCCGGGGTACATCGGTAACTGCCACTCTTCCGGGACGGTCATTCTGGATCAGATGGGCGCAGTGCACATGAAAACCgggcagcc                                                     >  2:301163/1‑245 (MQ=255)                                                                                                                                                                                        ggcgTACGCATGGGCGCCC‑GAAATGTCATCCGGTAAAGATCCCCCGTCGGGTCTCTGGAAAATCGCCGGTGTCCGGTTTCGGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGCGTTGCTGGATAACCTGGTCGAACGCGCTAGTCTGCCGGGTTACACCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAttt                                                <  2:662486/250‑1 (MQ=255)                                                                                                                                                                                           gTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCt                                             >  1:345062/1‑251 (MQ=255)                                                                                                                                                                                             aCGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGTAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGTTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCta                                            <  1:486503/250‑1 (MQ=255)                                                                                                                                                                                              cGCATGGGCGCAC‑GAAAGGCCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTTTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCtata                                          <  1:351673/251‑1 (MQ=255)                                                                                                                                                                                                              aaTGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCg                           >  1:350011/1‑251 (MQ=255)                                                                                                                                                                                                                tGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCgtg                         >  1:508139/1‑251 (MQ=255)                                                                                                                                                                                                                  tCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGtt                       >  1:648744/1‑251 (MQ=255)                                                                                                                                                                                                                   cATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGtt                       <  1:36325/250‑1 (MQ=255)                                                                                                                                                                                                                   cATCCGGTAAAGATACCCCGTCTGGTCACTGG‑AAATCGCCGGTGTCCCGGTTCTGTTTG‑CG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGATAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGTCGAAGAGCACATGAAAACCGGCAAGCCGATTTTATATACCTCCGCTGACTCCGTGTTcc                     >  2:25639/1‑251 (MQ=255)                                                                                                                                                                                                                          taaagctgccccgacaggTCACTGGCAAATCGCCGCTGGCCCGGTTCTGTTTG‑AG‑AGGGGCCATTTCTCCGATCACGAAAACAGCCTCCCGCAAGAGCTCCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCAATCTGGATCAGCTGGGCGAAGCGCACATAAAACCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCAGTGTTCCAGAtt                >  2:277375/17‑250 (MQ=255)                                                                                                                                                                                                                            aaGATACCCCGTGTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGGTTG‑AG‑TGGGGATATTTCCCCGATCACGAAAACAGCTTCCCGCAAGCGCTGCTGGATAAACTGGTCGAACGCGTTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcc             <  2:153044/251‑1 (MQ=255)                                                                                                                                                                                                                                  ccccGTCTGGTCACTGGGAAATCGCAGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATCTTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAATGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAACTTGCCTGCCac       >  2:501122/1‑250 (MQ=255)                                                                                                                                                                                                                                          ggccACTGGGAAATCGCTGGTGTCCCGTTTTTGTTGG‑AG‑TGTGGATATTTCGCCGATCACTAAAACGGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCGCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACgaagg  <  2:236002/245‑2 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            TGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTG‑AG‑TGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGA  >  NC_000913/4619671‑4620150

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