Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,067	A→G	Q155Q (CAA→CAG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,067	0	A	G	93.0%	109.5 / 0.3	43	Q155Q (CAA→CAG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (23/17);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

GGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACG  >  NC_000913/4619825‑4620296                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      ggAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAAGAAAGGTCATCCGGTAAAGATACCCCGGCTGGTCACTGGGAAAGCGCCGGGGTACAGGTTCTGTTTGAGTGGGGATATTTCTCCGATAACGAAAACAGATTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAAAGCGATAATCTGCCGGGTTACCTCGGAAACTGCCACTCTTCCGGAACGGTAATTCTGGAACAGCTgggggg                                                                                                                                                                                                                               >  2:59256/1‑248 (MQ=255)      ggACGGCAACGCTGAAGTTATCGGCGCGTACGCCTGGGCGCCCTGAATGGCATCCGATAAAATTAACCCGTCTGGTCACTGGGAAATCGCCGGGGTACCGGCTCTGTTTGAGTGGGGATGTTTCTCCGATCACGAAAACAGCTTCCCGCAAGCGCTGCTGGATAAAATGGT‑CGAACGCGCTACTCTGCCGGGTTACCTCGGAAACTGACAATCTTACGGCACGCTCATTCAGGCTCAGCTGGGCGAAGAGc                                                                                                                                                                                                                          >  2:425158/1‑251 (MQ=255)                    aaggcaTCGGCGCGGACGTATGGGCGCACGAAATGACCTCCGGAAAAGATACCCCGTCTGGTAACTGGGAAATCGCCGGAGTCCCGGTTCAGTTTGAGTGGGGATATTTCTCCGATCACGAAAACCACTTCCAGCACGAGCTGCTGGATAAAATGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGACACTATTCAGGTACGGGAATTCTGGAACAGCGGGGCGTAGAACACATGAAAACCgg                                                                                                                                                                                                             >  2:3585/6‑250 (MQ=255)                              cgcgTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTTTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTG‑ATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAtttt                                                                                                                                                                                                 <  1:219538/251‑1 (MQ=255)                               gcgaccGCATGGGCGCAAGAAATTGCACCCGGTAAAGATACCCCGTCGGGTCCCTGGGAAACCGCCGGTGTCCCGGTTCTGTTTGAGGGGGAATATTTCCCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGG‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAtttt                                                                                                                                                                                                 <  2:228601/246‑1 (MQ=255)                                   ccGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCTGTTTTGTTTGAGTGGGGATATTTCTCCGATCCCGAAAACAGCTTCCCGCAAGAGCTTCGGGATAAACTGGT‑CGAACGCGCTAATCTGCCTGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCta                                                                                                                                                                                              <  1:492891/249‑1 (MQ=255)                                             gcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCtata                                                                                                                                                                                            >  1:371693/1‑242 (MQ=255)                                             gcACGAAATGTCAATCGTAAAAGATACCCCGTCTGTGCACGGGGAAATCGCCGGGGGCCCGGTTCTGTTGGAGTGGGGATATTTCCCCGTGAACAAAAACTGCTTCCCGAAAGAGCTGCTGGATAACCTGGG‑CGAACGGGCTAATCTGCGTGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACAGGAAAACCGGCAAGCCGATTTTCtata                                                                                                                                                                                            <  2:371693/242‑1 (MQ=255)                                                         aTCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTcc                                                                                                                                                                       >  2:99616/1‑251 (MQ=255)                                                             ggTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTcc                                                                                                                                                                       >  1:7576/1‑247 (MQ=255)                                                                      cccccGTCTGTTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCCCCGATCACGAAAACAGCTCCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCa                                                                                                                                                          <  2:136150/250‑1 (MQ=255)                                                                                   cctggaaaaccgCCGGTGTCCCGGTTCTGTTTGAGTGGTGACATTTCTCCAATCACGAAAACAGCTTCCTGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAACCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAccttt                                                                                                                                             <  1:92634/240‑5 (MQ=255)                                                                                       ggAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATTAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTc                                                                                                                                         >  1:100512/1‑251 (MQ=255)                                                                                        aaaaTCGCCGGGGTCCCGGTTTTGTTTGAGTGGGGTTATTTCCCCGATCCCGAAAACAGCTTCCCGCAGGACCTGCTGGATAACCTGGT‑CGAACGCGCTAATCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAcctttggcct                                                                                                                                        <  2:320594/250‑10 (MQ=255)                                                                                         aaaTCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAcctttggcctg                                                                                                                                       >  1:174773/1‑241 (MQ=255)                                                                                                     ccccGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTCCCCGCAAGAGCTGCTGGATAAACTGTC‑CGAACGCGCAAATCTGCCGGTTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTAc                                                                                                                           <  1:267797/250‑1 (MQ=255)                                                                                                       ccGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTACCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACg                                                                                                                          >  1:450035/1‑250 (MQ=255)                                                                                                          gTTCTGTTTGAGAGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACt                                                                                                                      >  1:340907/1‑251 (MQ=255)                                                                                                              tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATCAACTATACGAACTGTGc                                                                                                                  >  1:476691/1‑251 (MQ=255)                                                                                                                         ggtttttcctcagtacACGAAAACGGCTTCCCCCAAGGTCGGCTGGTGAAACGGGT‑CTACCGGCCTAACCTGCCGGTTTCCCTGGGAAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGTGAAGAGCACATGAAAACCGGCAATCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCg                                                                                                       <  2:381693/236‑1 (MQ=255)                                                                                                                           atatTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCAGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATTAAAACCGGCAAGCCGATTTTCTATACCTCAGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGAAACACTCAACGAACTGTGCGAAATCGCCCggg                                                                                                     >  2:314170/1‑249 (MQ=255)                                                                                                                               ttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGa                                                                                                    >  1:231622/1‑248 (MQ=255)                                                                                                                                         aCGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGAGATACTCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCAAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAc                                                                                       >  1:108897/1‑251 (MQ=255)                                                                                                                                           gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                     >  1:126720/1‑251 (MQ=255)                                                                                                                                           gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGAAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                     >  1:474348/1‑251 (MQ=255)                                                                                                                                           aaaaaCATCTTCCCGCAAGTGCTGCTGGAAAAAC‑GTTACGAACCCTCTATTCGGCCGGGTTACCCCGGAAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcg                                                                                      <  2:474348/249‑1 (MQ=255)                                                                                                                                                           aaggccgggtggaaaaccTGGT‑CGAACGCCCTAACCGCCCGGGTTACCCGGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGt                                                                     <  2:370539/234‑1 (MQ=255)                                                                                                                                                           aaGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGACGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATACACTCTACGAACTGTGCGCAATCGCCCGTGAAGCGCTGACCACCGGCGGCTAAAACATCggg                                                                     >  2:87975/1‑250 (MQ=255)                                                                                                                                                            agagCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTc                                                                    >  1:388624/1‑251 (MQ=255)                                                                                                                                                             gagCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAGGGAACCTTTGGCCTGGATAAACTCTACGACCTGTGCGAAATCGCCCGTGAAGAGCTGCCCAACggaggcacacacacccgtcg                                                                   >  2:260362/1‑237 (MQ=255)                                                                                                                                                                     ggATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCg                                                           >  1:17005/1‑251 (MQ=255)                                                                                                                                                                        tAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGCCTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGATAccacctccgtctggatatcaccc                                                        >  2:459015/1‑229 (MQ=255)                                                                                                                                                                                         ccTAACCTGCGGGGTTACCTCGGTAACTGCCACCCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTACTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCTAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTCATCGCTCGTCCGCTTATCggtga                                        <  2:454817/250‑4 (MQ=255)                                                                                                                                                                                               cggCCGGTTTACCCCGGCACCTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCGGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGc                                  <  2:388624/249‑1 (MQ=255)                                                                                                                                                                                                         acctcgtaacctgccaacctTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGTCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCTTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACtt                         <  2:225887/232‑1 (MQ=255)                                                                                                                                                                                                              ggTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCg                   <  1:99616/251‑1 (MQ=255)                                                                                                                                                                                                               gTAACTGCCACTCTGCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGt                  <  1:233809/251‑1 (MQ=255)                                                                                                                                                                                                                      ccACACTTCCGGTACTGTAATTATGGATCAGCTGGGCAAAGAACACATGAAAACCGCCAAGCCGAGTTTCTATACCTCCGCTGACTCCGTGTCCCAGATTGCCTGCCACGAAGGCACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCATCGGCGGCTACAACATCTGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCgtta           <  2:116511/251‑4 (MQ=255)                                                                                                                                                                                                                            ttCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGCTTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGACGGTAACTTCCTGCGTACCGGTaaaccac     >  2:78625/1‑246 (MQ=255)                                                                                                                                                                                                                              ccGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCAc   <  1:258111/251‑1 (MQ=255)                                                                                                                                                                                                                              ccGGTACGGGCATTCTGGAGTAGTGGGGCTAAGAGCCCTTGAAAACCGTCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCCCGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGTCGGTA‑CTTCCAGCGTACCGGTAACCGTCACg  <  2:373706/251‑1 (MQ=255)                                                                                                                                                                                                                               cGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCATCGTACCGGTAACCGTCAc   >  1:373389/1‑250 (MQ=255)                                                                                                                                                                                                                                     ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCggtga                                        >  1:124045/1‑204 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      GGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACG  >  NC_000913/4619825‑4620296

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