Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,243,185	C→T	T156T (ACG→ACA)	folE ←	GTP cyclohydrolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,243,185	0	C	T	98.6%	269.3 / NA	73	T156T (ACG→ACA)	folE	GTP cyclohydrolase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/35);  minor base G (0/1);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGCGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACC  >  NC_000913/2242941‑2243421                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             ataaTCGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGCCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGCTTTGAACAATCCACCAAGCGAGGTCGTTGTCGTGGCCCTGGTTGCATCGCGGATTCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACCGCCACGTTATTGGTGCCCACCAGCGTTTGTAGCGCCCTACGAATTTGCTGTgtccg                                                                                                                                                                                                                                         >  2:303‑M2/5‑248 (MQ=255) ataaTCGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCACTAAGAATTTGCTGTGTc                                                                                                                                                                                                                                           >  2:400567‑M2/5‑248 (MQ=255)   aaTCGCGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCCCAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACg                                                                                                                                                                                                                                      >  2:586517/1‑251 (MQ=255)      cGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGtt                                                                                                                                                                                                                                    <  1:479939/250‑1 (MQ=255)      cGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTc                                                                                                                                                                                                                                   <  2:191923/251‑1 (MQ=255)      cGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTc                                                                                                                                                                                                                                   >  2:412387/1‑251 (MQ=255)          cGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGCAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATCAGAATTTGCTGTGTCAGACGTTCCTGc                                                                                                                                                                                                                               >  2:534673/1‑251 (MQ=255)              aCGTTGCGCTCCATTGTTCCTGCCGTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTTTGTCAGACGTTCCTGCACCt                                                                                                                                                                                                                           <  1:443883/251‑1 (MQ=255)                gTTGCGCTACATTGTTCCTGCCTTTTACTCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGc                                                                                                                                                                                                                         >  2:224682/1‑251 (MQ=255)                   gcgcTCCATTGTTCCTGCCTTTTAATCAGTGGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTAGTCGGACTGGATTTGAACAATCCACCAAGAGATGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTTCTGTGTCAGACGTTCCTGCACCTGCgg                                                                                                                                                                                                                       <  2:389060/250‑1 (MQ=255)                             gTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGACAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCa                                                                                                                                                                                                             >  1:74567/1‑250 (MQ=255)                               tCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTTGTTGCATCGCGGATGCCACGCGCCTGCACGCAGTAATGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGGTCCGGCGCCTGCGGACGCTGGGCaaa                                                                                                                                                                                                           <  2:391246/250‑1 (MQ=255)                                            tCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGc                                                                                                                                                                                             >  1:532523/1‑251 (MQ=255)                                               gTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGt                                                                                                                                                                                          <  2:576189/251‑1 (MQ=255)                                                        aCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAAAAATACACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGCCGCTGGGCAAAGAACTGCACAATGCGGTTACTTTTTg                                                                                                                                                                                 >  2:193667/1‑251 (MQ=255)                                                          gCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGCACTGCACAATGCGGTTAATTTTTGAc                                                                                                                                                                               >  2:306019/1‑251 (MQ=255)                                                          gCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGAc                                                                                                                                                                               >  1:140750/1‑251 (MQ=255)                                                           cacaGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAACGAACTGCACAATGCGGTTAATTTTTGAc                                                                                                                                                                               >  2:456649/1‑250 (MQ=255)                                                            ccaGCGCGAAGAAACTTGTGGCGCGTATTCTGGCTCGATGTGAACAATCCCACGAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCTTGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCGGCGTTTGTAGCGCAATAAGAATTCTCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAAGGCGGTTAATTTTTGACAg                                                                                                                                                                             <  2:525682/250‑1 (MQ=255)                                                              aGCGCGCAGAAACTCGTGGCGCGTATTCGGCATGGATTTCAACAATCAACCAAGAGCGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACCGCCCCGTTATTGGGTCCCAGCAGCGTTTGTAGCGCAATCAGACTTTGCTGTGACAGACGTTCCTGAACCTGCGCAAGCTGGGCAAAGAACTGCACCCTGCGGTTAATTTTTGACAgcc                                                                                                                                                                           >  2:499276/1‑249 (MQ=255)                                                                     aGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGGTTGTAGCGCAATAAGAATTTGCTGTGTCATACGTTCCTGCACCTGCGGACGCTGGGCAATGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCAc                                                                                                                                                                    <  1:351828/251‑1 (MQ=255)                                                                       aaaCTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGGCGTTGTCGTGTCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTGTGACATACCGATCACCg                                                                                                                                                                  <  2:341215/251‑1 (MQ=255)                                                                        aaCTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCg                                                                                                                                                                  >  2:202487/1‑250 (MQ=255)                                                                           tCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGCTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATc                                                                                                                                                              >  1:549567/1‑251 (MQ=255)                                                                            cGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCt                                                                                                                                                             >  2:380296/1‑251 (MQ=255)                                                                              gggCGCGTATGCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCttt                                                                                                                                                           <  2:5567/250‑1 (MQ=255)                                                                                  gcgTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCggg                                                                                                                                                       >  2:389352/1‑251 (MQ=255)                                                                                  gcgTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCggg                                                                                                                                                       >  2:212582/1‑251 (MQ=255)                                                                                    gTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGat                                                                                                                                                     >  1:562423/1‑251 (MQ=255)                                                                                           gACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAg                                                                                                                                                 <  2:317552/248‑1 (MQ=255)                                                                                               ggATTTGAACAATCCACCAAGAGGGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCAc                                                                                                                                            <  2:306049/249‑1 (MQ=255)                                                                                                       aCAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTg                                                                                                                                  >  2:495520/1‑251 (MQ=255)                                                                                                        cAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCttt                                                                                                                                   >  1:296366/1‑249 (MQ=255)                                                                                                        cAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCttt                                                                                                                                   <  2:296366/249‑1 (MQ=255)                                                                                                         aaTCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGcc                                                                                                                                >  1:169410/1‑251 (MQ=255)                                                                                                         aaTCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGcc                                                                                                                                <  2:444480/251‑1 (MQ=255)                                                                                                          aTCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGcc                                                                                                                                >  2:451774/1‑250 (MQ=255)                                                                                                          aTCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGGTAATTTTTGACAGACCGATCACCGAATCTTTCTGGATAGATGCCACCGTCGCTTTGCCa                                                                                                                               <  2:76316/251‑1 (MQ=255)                                                                                                          aTCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGATGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGcc                                                                                                                                >  1:524134/1‑250 (MQ=255)                                                                                                                caAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATg                                                                                                                         >  1:247870/1‑251 (MQ=255)                                                                                                                 aaGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATgg                                                                                                                        >  1:90610/1‑251 (MQ=255)                                                                                                                     gagGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAAc                                                                                                                    <  1:586517/251‑1 (MQ=255)                                                                                                                      agGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTc                                                                                                                                        <  1:553650/230‑1 (MQ=255)                                                                                                                      agGTCGTTGTCGTGGCACAGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGAACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGCCGTTCCTGCAACTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTc                                                                                                                                        >  2:553650/1‑230 (MQ=255)                                                                                                                                                  gcgGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATGTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCagag                                                                                       <  1:540015/251‑1 (MQ=255)                                                                                                                                                         ccACGCGCCTGCAATCAGTAATGCACCGCTTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGTGTTTGTAGCGCACTAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATAACCGAATATTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGTGGGTCAGAGCGATATc                                                                                <  2:74567/251‑1 (MQ=255)                                                                                                                                                               gcCTTCACGCAGTAAGGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCAc                                                                          <  2:140750/251‑1 (MQ=255)                                                                                                                                                                 cTTCACGCAGTAATGCACCCCGTCGATCGGGCCAGCCACGTTATTGGTACCCTGCAGCGTTTGTAGCGAAATAAGATTTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGTACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGCCGCTTTGCCATCGTTGGTAACAACATGGTGTTCACAGGTGCTGGTCAGTGTGATATCGCGCACgg                                                                        >  2:575127/1‑251 (MQ=255)                                                                                                                                                                    cACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTTCTGGTCAGAGTTATATCGCGCGCGGTg                                                                      <  1:193667/250‑1 (MQ=255)                                                                                                                                                                      cGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGAcc                                                                   <  1:534673/251‑1 (MQ=255)                                                                                                                                                                         aGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCAtt                                                                >  1:346461/1‑251 (MQ=255)                                                                                                                                                                            aaTGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTc                                                              <  2:562423/250‑1 (MQ=255)                                                                                                                                                                            aaTGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTc                                                              >  1:538797/1‑250 (MQ=255)                                                                                                                                                                                  ccGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTTGGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGACATCGTCGCTTTGCCAGCGATGGTAACAAAATGGTGTTCTGAGGTGGTGGCCAGAGTGATATCGCGCACGGTGACCATTTCATCGAcc                                                       <  1:575127/251‑1 (MQ=255)                                                                                                                                                                                     cgTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGGTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTc                                                    <  2:524134/251‑1 (MQ=255)                                                                                                                                                                                        cgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCAt                                                  <  2:181515/250‑1 (MQ=255)                                                                                                                                                                                            cgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTGTTTACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGGCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTg                                              <  1:212582/250‑1 (MQ=255)                                                                                                                                                                                               gaCAGCCACGTTATTGGTGCCAAGCAGCGTTTGTAGAGCAATAAGAATTTGCTGTGTCAGAAGTTCCTGCACCTGCGGACGCTGGGCAAAGAGCTGCACAATGCGGTTACTTTTTGACAGACCGATCACAGACTCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGGCACAAAATGGTGTTCCCAGGTGCTGGTCAGAGTGCTATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGtttt                                          >  1:565020/1‑251 (MQ=255)                                                                                                                                                                                                aCAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCtt                                                                                                                                    >  2:465965/1‑160 (MQ=255)                                                                                                                                                                                                aCAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGGCGCtt                                                                                                                                    <  1:465965/160‑1 (MQ=255)                                                                                                                                                                                                     cACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGACAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATgg                                    >  2:274883/1‑250 (MQ=255)                                                                                                                                                                                                              ggTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGAttt                           >  1:349111/1‑251 (MQ=255)                                                                                                                                                                                                                  cccAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTTATTTTCgg                       <  1:456649/251‑1 (MQ=255)                                                                                                                                                                                                                    cagcagCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGa                     <  1:306019/251‑1 (MQ=255)                                                                                                                                                                                                                           gTTTGTAGCGCAATAAGCCTTTGCTGTGTCAGCCGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATCGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGc              >  2:272391/1‑251 (MQ=255)                                                                                                                                                                                                                             ttGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCg             <  1:400567/250‑1 (MQ=255)                                                                                                                                                                                                                               gTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTaa          <  1:389352/251‑1 (MQ=255)                                                                                                                                                                                                                                   cgcAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCa      <  1:380296/251‑1 (MQ=255)                                                                                                                                                                                                                                     cAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGa    <  2:532523/251‑1 (MQ=255)                                                                                                                                                                                                                                      aaTAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGAc   <  1:473691/251‑1 (MQ=255)                                                                                                                                                                                                                                       aTAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGCCGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCTTCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTg                               >  2:267648/1‑222 (MQ=255)                                                                                                                                                                                                                                       aTAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCGCCGAATCGTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAAGGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTg                               <  1:267648/222‑1 (MQ=255)                                                                                                                                                                                                                                       aTAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGAcc  >  1:351273/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             AAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGCGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACC  >  NC_000913/2242941‑2243421

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