Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,400,207	G→A	T294I (ACC→ATC)	mreB ←	dynamic cytoskeletal protein MreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,400,207	0	G	A	94.7%	99.5 / ‑2.3	39	T294I (ACC→ATC)	mreB	dynamic cytoskeletal protein MreB
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (16/20);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

TTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGG  >  NC_000913/3399974‑3400447                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 ttcttAAGTTATGAGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCATATTACTCTTCGCTGAACAGGTCGCCGCAGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCCCCGCGCGCCACACAGGCCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGATTCTACCATTACCAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGACGAGCCCCATgacaccc                                                                                                                                                                                                                                 >  2:189636/5‑245 (MQ=255)    aTAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCg                                                                                                                                                                                                                              <  1:150255/251‑1 (MQ=255)               gTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCgg                                                                                                                                                                                                                    >  1:317379/1‑250 (MQ=255)                         aTCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTcc                                                                                                                                                                                                         <  1:222494/251‑1 (MQ=255)                              aCCAGGCAGGGTAAACAGACACTTCCCCTGCCGGCAGCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCcggcgg                                                                                                                                                                                                    <  1:337644/251‑1 (MQ=255)                                    caggGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGTCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTg                                                                                                                                                                                              <  1:249278/251‑1 (MQ=255)                                      gggTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGtt                                                                                                                                                                                            >  1:274216/1‑251 (MQ=255)                                      gggTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGtt                                                                                                                                                                                            >  1:267805/1‑251 (MQ=255)                                            aCAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTg                                                                                                                                                                                      <  1:254675/251‑1 (MQ=255)                                                ccacTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAAc                                                                                                                                                                                  <  1:70588/250‑1 (MQ=255)                                                    ttCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGCTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAAccctt                                                                                                                                                                              >  1:243923/1‑248 (MQ=255)                                                                cATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTTCAACCCTTACCGCGCTcac                                                                                                                                                                   >  1:265478/1‑250 (MQ=255)                                                                   ccgagttctctTCCCTGAACAGGCTGCCGCCGTGTATGTCAATTATTTCCAGCGCTTTGCCGCCTCCGCGCGCCACACAGGTCATCGGGTCTGCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCTCCGATTAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTGCCAGTGCAACCATTACCGCGCTCACAAt                                                                                                                                                                <  2:337738/243‑1 (MQ=255)                                                                         aCTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGt                                                                                                                                                          >  1:54636/1‑250 (MQ=255)                                                                            ctTCGCTGAACAGGACGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGGCAGc                                                                                                                                                      >  1:62513/1‑251 (MQ=255)                                                                              tCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAg                                                                                                                                                       <  1:9755/248‑1 (MQ=255)                                                                                    aaCAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTg                                                                                                                                              >  1:195268/1‑251 (MQ=255)                                                                                               gccgTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGa                                                                                                                                   <  1:374224/251‑1 (MQ=255)                                                                                                  ttGCATGTCGATCATTTCCAGCGCTTTGCCGACACCGCGCGCCACACCGGTCAGCGGGTCTTCAGCAACAACGACTGGACTGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGAAGCAGTGCGCCACAACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCAATGTTCCTGGCCAAACAATACCGCGCTCACACTACCGGCCAGCGGGTCCTGCAggccttggggg                                                                                                                                 >  2:130536/2‑240 (MQ=255)                                                                                                    gCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTATTCCATTAACAAACGATCAAGGTTACGAAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCAGCAATGCTTCGAGGAtt                                                                                                                              >  1:51175/1‑251 (MQ=255)                                                                                                          cggtcttttCCAGCGCTTGGCCGCCCCCGCGCGCCACACAGGTCAGCGGGTCCTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGAGGTCGGAAGCCAGTTCCGGCGGGCTCTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCAtt                                                                                                                         <  2:195268/244‑1 (MQ=255)                                                                                                                                          gcCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTc                                                                                        >  1:319263/1‑251 (MQ=255)                                                                                                                                           ccACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGAGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCt                                                                                       >  1:251191/1‑251 (MQ=255)                                                                                                                                                catgtctgCGGGTCCTCAGCAACACCCACTGTAATGCCGCTTACGTCCAATAACAACCGGTCAAGTTTACGCAGCAGTGCGCCCCCACCTATGATCACCATTCCGCGCTCGGAGATGTCGGAAGCCAGTTCCCGCGGGCACTGTTCCAGTGAAACCATTACCGCGCTCTCAATACCGGTCAGCGGTTCCTGCAGTGCTTTGTGGATTTCATTTGAGTTCAGGGTAAAACCGCGTGGAACACCTTCGGCCAg                                                                                  <  2:75169/244‑1 (MQ=255)                                                                                                                                                    gaagtggggatacAGCAACAACGACTGGAATGCCGGTTTCTTCCACTAAAAAACGGTCAAGGTTCCACAGCCGTGCGCCACCACCGATGAGCACCATGCCTCGCTCGGAGATGTCGGTAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCATGCAGTGCTCCGAGGATTTCATTGTAGTCCAGGGTATAACCGCGGGGAACACCTTCTTACAGGTTa                                                                              <  2:243923/239‑1 (MQ=255)                                                                                                                                                      aGCGGGTCGTCAGCAACAACGCCTGGAACTCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCTCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGAAACCATTACCGCGCTCACAATTCCGGTCATCTGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCATGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACg                                                                            <  2:36330/251‑1 (MQ=255)                                                                                                                                                                     acaacttcTGGAATGCCTGTTTCTTCCATTAACAAACGGTCAAGGTTCCTCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCTGTCAGCGGTTCCTGCAGTGCTTCGAGTATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGAt                                                             <  1:17972/244‑1 (MQ=255)                                                                                                                                                                     acaacGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGAt                                                             >  1:146302/1‑251 (MQ=255)                                                                                                                                                                                aaTGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACtt                                                  <  1:168787/251‑1 (MQ=255)                                                                                                                                                                                 aTGCCGGTTTCTTCCATTAACAACCGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCTGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTc                                                 <  1:313548/251‑1 (MQ=255)                                                                                                                                                                                    ccGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATc                                              <  1:200543/251‑1 (MQ=255)                                                                                                                                                                                    ccGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATc                                              >  1:120703/1‑251 (MQ=255)                                                                                                                                                                                            ttcCATTAACAAACGGTCAAGGTTACGAAGCAGTGCGCCACCACCGGTGAGCACAATGCCGCGCTCGAAGATGTCGGAAGCCAGTTCCGGCGGGCACTGATCCAGTGCAACCATTACCGCGCTCACAATACCGGTAACCGGTTCCTGCAGTGCTACGAGGATTTAGTTGGAGTTCAGGGGCAAACCCCGGGGAACCCCTTCTGCCAGAGTACAGACACGAACTTCGAACTCACCGACTTCCTAGCCCggct                                      >  2:198986/1‑249 (MQ=255)                                                                                                                                                                                             tcctttaccaagccgTCAAGTTTACGCAGCAGTGCCCCCCCACCGATGAGCACCATGCCGCGCGCGGAGATGTCGGAATCCATTTCCGGCGGGCACTGTTCCAGTGCAACCAGTCCCGCGCTCACAATACCGGTCAGCGGCTCCTGCAGTGCTGCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGAt                                      <  2:302837/236‑1 (MQ=255)                                                                                                                                                                                              ccATTAACCAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATaa                                    <  1:13833/251‑1 (MQ=255)                                                                                                                                                                                                                  ttACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGGTTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCtt                >  1:84045/1‑251 (MQ=255)                                                                                                                                                                                                                      gcagcaGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATa            <  1:307442/251‑1 (MQ=255)                                                                                                                                                                                                                          cagTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCCGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTAAGGCCACGAACTTCGATTTCACGGACTTCATCACCCGGCTAAACCGAACCGATTTCGTGCTTGATACGtt        >  1:354994/1‑251 (MQ=255)                                                                                                                                                                                                                                gcCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGAGCACTGGTCCAGTGCAACCATTACCGCGCTCACAATACCGTTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGGTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCgg  >  1:306680/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 TTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGG  >  NC_000913/3399974‑3400447

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