Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	694,176	G→A	S260S (AGC→AGT)	miaB ←	isopentenyl‑adenosine A37 tRNA methylthiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	694,176	0	G	A	98.5%	220.0 / NA	68	S260S (AGC→AGT)	miaB	isopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/31);  minor base C (0/1);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGATAAAGCTGTAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAGCTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCGCAAGCTGGGCAATCTCAAACAGAATATCGT  >  NC_000913/693944‑694420                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      agaTAAAGCTGTAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCgtg                                                                                                                                                                                                                                     >  1:39441/1‑251 (MQ=255)  gaTAAAGCTGTAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCgtgt                                                                                                                                                                                                                                    <  2:16533/251‑1 (MQ=255)   aTAAAGCTGTAGCTCATTTCTAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGTGAAGCCAACGATGAAATCAGAACTGATCTTAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTCCTTCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGCCAGATGCAGGAACCTCACCAGCTCCGGCGTGTc                                                                                                                                                                                                                                   <  2:143850/251‑1 (MQ=255)   aTAAAGCTGTAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTTATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTGTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTc                                                                                                                                                                                                                                   <  2:263435/251‑1 (MQ=255)   aTAAAGCTGTAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACTGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTc                                                                                                                                                                                                                                   <  1:50448/251‑1 (MQ=255)     aaaGCTGTAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGAGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTcgc                                                                                                                                                                                                                                 >  1:431756/1‑251 (MQ=255)        gCTGTAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGAt                                                                                                                                                                                                                              <  1:282668/251‑1 (MQ=255)          tGTAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCCCCGGCGTGTCGCGATac                                                                                                                                                                                                                            >  1:414808/1‑251 (MQ=255)            tAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATacac                                                                                                                                                                                                                          >  1:63940/1‑251 (MQ=255)             aGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGATGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCAAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGAAGGAAACTCCACAGCTCCGGCGTGTCGCGAGACACt                                                                                                                                                                                                                         >  2:157889/1‑251 (MQ=255)             aGCTCATGTCGAAATTAACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTACTCGGTGGATTCACCAGGGCAGCCAACGATGAACTCAGAACTGTTCTGAAAATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCGTTGTACTCCACCGCCGGATGGGGACGGCACATAAGGTTCAGAATGAGATCGGACCCGATCTGTACCGGCAGATGCAGGAAACGAAACAGCTCCGGGGTGGCGCGATaaacc                                                                                                                                                                                                                         >  1:191896/1‑247 (MQ=255)                 cATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCg                                                                                                                                                                                                                      >  2:40172/1‑250 (MQ=255)                      cGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCgatgat                                                                                                                                                                                                                 >  1:169104/1‑250 (MQ=255)                       gAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGatat                                                                                                                                                                                                               <  1:116083/251‑1 (MQ=255)                        aaaTTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACCCTTCGATGATAtc                                                                                                                                                                                                              >  2:449992/1‑251 (MQ=255)                             gACGTCGGCAAGCAGCTTCATCGTTGTCTCGAAGTCTTCGGTGGTTGCGCCAGGGAAGCCAACGATGAAGTCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTCCCGGCAGATGCAGTAATCTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATAtcgtcg                                                                                                                                                                                                          <  2:300502/250‑1 (MQ=255)                                  cGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGaa                                                                                                                                                                                                     <  1:239260/250‑1 (MQ=255)                                   ggCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCGCCAGCTCCGGCGTGTCGCGATACACTTCGATGAAATCGTCGGGGAAtt                                                                                                                                                                                                   >  2:101169/1‑251 (MQ=255)                                    gCAATCAGCTTCAGCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGATACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTc                                                                                                                                                                                                  <  1:309665/251‑1 (MQ=255)                                      aaTCAGCTTCATCGTTTTCTAGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCAACAGCTCCGGCGTGGCGCGCTACACTTCGATGATATCGTCGGTGAATTCgg                                                                                                                                                                                                >  2:51199/1‑250 (MQ=255)                                         cAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCg                                                                                                                                                                                             >  1:160971/1‑251 (MQ=255)                                           gCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGa                                                                                                                                                                                           <  1:278563/251‑1 (MQ=255)                                                aTCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCGGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCt                                                                                                                                                                                      <  2:214302/251‑1 (MQ=255)                                                  cGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCtg                                                                                                                                                                                     <  1:101169/250‑1 (MQ=255)                                                    ttttCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCtggtg                                                                                                                                                                                  >  2:399075/1‑251 (MQ=255)                                                      ttCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCtggtgg                                                                                                                                                                                 <  1:115512/250‑1 (MQ=255)                                                                   cGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATACCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATAcga                                                                                                                                                                   <  2:443969/251‑1 (MQ=255)                                                                     gtggtTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATCTAGTCGGTGCATTCGATCGGATGGCTGGTGGTAAAGCGAATAcgatc                                                                                                                                                                 >  2:129602/1‑251 (MQ=255)                                                                        ggttCGCCAGGGAAGCCAACGATGAAATCAGAACTGAGCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTGCGCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACGCTTCGATGATATCGTCCGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATAcgatcgat                                                                                                                                                              <  2:205932/249‑1 (MQ=255)                                                                           tCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATcc                                                                                                                                                            <  1:251500/250‑1 (MQ=255)                                                                            cGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCTCGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCATGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCg                                                                                                                                                          <  1:65086/251‑1 (MQ=255)                                                                                 gggAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGAt                                                                                                                                                     <  1:449992/251‑1 (MQ=255)                                                                                    aaGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCgc                                                                                                                                                  >  2:114507/1‑251 (MQ=255)                                                                                         aaCGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCgcagca                                                                                                                                              >  2:49607/1‑250 (MQ=255)                                                                                               gAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGAGCGCAGCAACCAGAc                                                                                                                                       <  1:399075/251‑1 (MQ=255)                                                                                                 aaTCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATAAGATCGATCCCGTCGATCGAAGCAACCAGACgc                                                                                                                                     >  2:390302/1‑251 (MQ=255)                                                                                                   tCAGAACTGCTCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGCTCGATCCCGTCGATCGCAGCAAcccgccgc                                                                                                                                     >  2:20818/1‑243 (MQ=255)                                                                                                                 aaTATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGAGCCGTATGGGTACGGCCCATCAGGGTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGATGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACg                                                                                                                    >  1:219/1‑251 (MQ=255)                                                                                                                    atCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCCGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCCAACGATc                                                                                                                 >  2:270938/1‑251 (MQ=255)                                                                                                                    atCTGAACCCGCCGCACGCAGTTTACGGAAGACTGCTTTCTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGCCCGGAACCCCTCTGTACCGGCAGAAGCAGGAAACTCACCAGCTCCCGCCTGTAGCGTTACACTTCGAAGTTTTCGCCGGTAAATTCGCTTGGATAGCTGGGGGTTATGCGAATACGAACGATCCCGTCTAACGGAGCAACCAGACgcagccggacggtaaacgctc                                                                                                                  >  2:354685/1‑235 (MQ=255)                                                                                                                           cgcgCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCggtggt                                                                                                          <  1:370148/251‑1 (MQ=255)                                                                                                                           cgcgCCGCACGCAGGTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACTATCCggtggt                                                                                                          <  2:197610/251‑1 (MQ=255)                                                                                                                                cgcacgcaGTTTACGGATGATAGCTTTGTACTCCAGCGCCGTATGGGTAAGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGAACCGGCAGAAGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATAACGTCGGTGAATTCGATCGGATGGCTGGCGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGACACGGTGGTGCcgt                                                                                                     >  2:294106/1‑251 (MQ=255)                                                                                                                                        gTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTAGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGGGGTGCCGTCGTAGttt                                                                                             >  2:245705/1‑251 (MQ=255)                                                                                                                                             cGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTcacc                                                                                        <  2:63940/251‑1 (MQ=255)                                                                                                                                                        tttGTCCTCCAGCGCCGTATTGTTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCTGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCTGATGGCTGGTGGTAATGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCg                                                                              <  2:295206/250‑1 (MQ=255)                                                                                                                                                                ccAGCGCCGTATGGGTACGCCCCATCAGGGTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTc                                                                     <  2:335734/251‑1 (MQ=255)                                                                                                                                                                 cAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCt                                                                    <  1:390302/251‑1 (MQ=255)                                                                                                                                                                     gcCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCCGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCCGTGGTGCCGTCGTAGTTTTCACCCCGCCAGGCGTTCACGATCTGAcc                                                                >  2:195053/1‑251 (MQ=255)                                                                                                                                                                     gcCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGAcc                                                                <  2:351055/251‑1 (MQ=255)                                                                                                                                                                        gTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTAAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCg                                                               >  1:307251/1‑249 (MQ=255)                                                                                                                                                                        atATGGGTACGGCCGATCAGGTTCATAATTCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCCCCAGCTCCGCCATGTCGCGATACCCTTCGATGATATCGTCGGTGCATTCGATAGGATGGCTGGTGGTAAATAGAATACGATCGCTACCCTCGATCGCATCAACCAGACGC‑GCAAGGTCGGCAAACTAGCCCGCAGTTCCACCGTAGTTTTCCCCACGCCATGCTTTCCCGTTCTTACCg                                                               >  2:357701/2‑249 (MQ=255)                                                                                                                                                                                     ccATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTc                                                 >  2:111798/1‑250 (MQ=255)                                                                                                                                                                                         cAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTcacgc                                             >  1:327180/1‑250 (MQ=255)                                                                                                                                                                                                 gAATTCGATCGGAACCGATCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGGCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCTATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGa                                    <  1:281868/251‑1 (MQ=255)                                                                                                                                                                                                 gAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTg                                     >  1:55416/1‑250 (MQ=255)                                                                                                                                                                                                   aTGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGc                                  <  1:457276/251‑1 (MQ=255)                                                                                                                                                                                                     gCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCg                                <  1:136647/251‑1 (MQ=255)                                                                                                                                                                                                       gATCGGAACCGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATAAGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCGcc                              >  2:369900/1‑250 (MQ=255)                                                                                                                                                                                                               ccGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCGCAAGCTGGGc                      >  2:45770/1‑251 (MQ=255)                                                                                                                                                                                                                cGCTCTGTACCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCGCAAGCTGGGCa                     <  2:169104/251‑1 (MQ=255)                                                                                                                                                                                                                        aCCGGCAGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATTCGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCGCAAGCTGGGCAATCTCaa              <  1:49607/250‑1 (MQ=255)                                                                                                                                                                                                                              aGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCGCAAGCTGGGCAATCTCAAACAGAat       >  1:113865/1‑251 (MQ=255)                                                                                                                                                                                                                              aGATGCAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGATGGTGGTAAAGCGAATACGATCGATACCGTCGATCGCAGCAAACAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGACGCAAGCTGGGCAATCTCAAACAGaa        >  1:418628/1‑250 (MQ=255)                                                                                                                                                                                                                                   cAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCGCAAGCTGGGCAATCTCAAACAGAATAtcgt  <  2:167561/251‑1 (MQ=255)                                                                                                                                                                                                                                   cAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCGCAAGCTGGGCAATCTCAAACAGAATAtcgt  >  1:438660/1‑251 (MQ=255)                                                                                                                                                                                                                                   cAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTT‑CACCACg                                                                                     <  1:99770/167‑1 (MQ=255)                                                                                                                                                                                                                                   cAGGAAACTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTT‑CACCACg                                                                                     >  2:99770/1‑167 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      AGATAAAGCTGTAGCTCATGTCGAAATTGACGTCGGCAATCAGCTTCATCGTTTTCTCGAAGTCTTCGGTGGTTTCGCCAGGGAAGCCAACGATGAAATCAGAACTGATCTGAATATCTGGACGCGCCGCACGCAGTTTACGGATGATCGCTTTGTACTCCAGCGCCGTATGGGTACGGCCCATCAGGTTCAGAATGCGATCGGAACCGCTCTGTACCGGCAGATGCAGGAAGCTCACCAGCTCCGGCGTGTCGCGATACACTTCGATGATATCGTCGGTGAATTCGATCGGATGGCTGGTGGTAAAGCGAATACGATCGATCCCGTCGATCGCAGCAACCAGACGCAGCA‑GATCGGCAAACGATCCGGTGGTGCCGTCGTAGTTTTCACCACGCCAGGCGTTCACGTTCTGACCGAGCAGGTTGACTTCACGCACGCCCTGAGCCGCAAGCTGGGCAATCTCAAACAGAATATCGT  >  NC_000913/693944‑694420

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