Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	21,895	C→T	R163R (CGC→CGT)	ribF →	bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	21,895	0	C	T	98.5%	195.7 / NA	66	R163R (CGC→CGT)	ribF	bifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (34/31);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACA  >  NC_000913/21648‑22133                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               aCCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGCCAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGAGAAGGGGGCGTGCGTATc                                                                                                                                                                                                                                             >  1:582087/1‑251 (MQ=255)     tGCAGAGTGTAGCGTTGCTTACATGCTGTGCGTGCGTTTCGACAGGCGTTTCTCTGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTAAAGATTTTGCGCGTAAAATTTCTTGCCGTAGGTGATACTTTCCGCTTTGGCGCTGGTCGTGAAGGCGCTTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAACCTTTTCGCGAATGTGGCGGGCGTATcagc                                                                                                                                                                                                                                          >  2:586820/1‑250 (MQ=255)        agagGGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCCGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGTCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAgcagca                                                                                                                                                                                                                                      <  1:400775/251‑1 (MQ=255)                    tGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTc                                                                                                                                                                                                                          >  1:155567/1‑251 (MQ=255)                     gATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTGCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTc                                                                                                                                                                                                                          <  1:530511/250‑1 (MQ=255)                        tACGCGCTGGGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTATGGTGAAGCTTTTGCGCGTCAAATTTCTAGCCGTAGGGGATGCTTTCCGCTTTGGCGCTGGTCGTGTAGGCGATTTCTTGTTATTACAGAAAGATGGCATGGAATACGGCTTCGATATAACCAGTACGCAAACTTTTTGCGAGGGAGGCGTGCGTATCAGCAGCAcagctgtacggcaagc                                                                                                                                                                                                                      >  2:526148/1‑236 (MQ=255)                         aCGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGc                                                                                                                                                                                                                      <  1:228554/250‑1 (MQ=255)                            tGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCCCCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCCGGcgcgg                                                                                                                                                                                                                  >  2:136620/1‑247 (MQ=255)                             gCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCTCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGcgctg                                                                                                                                                                                                                  <  1:356426/250‑5 (MQ=255)                             gCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCACACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGcgctgg                                                                                                                                                                                                                 >  1:521275/1‑246 (MQ=255)                                        gtTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTTGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGTCGCTGGCGGATGACAAt                                                                                                                                                                                                      <  2:295247/251‑1 (MQ=255)                                                 ggCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAAtctggctctg                                                                                                                                                                                             <  1:472167/251‑1 (MQ=255)                                                   cGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAAtctggctctggc                                                                                                                                                                                           <  1:275416/251‑1 (MQ=255)                                                     tttCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGCTTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAAtctggctctggc                                                                                                                                                                                           >  1:257506/1‑249 (MQ=255)                                                         gcggcgTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCagaa                                                                                                                                                                                       >  1:498970/1‑248 (MQ=255)                                                             ggTTAACCGAGCAAAATTTCATCAGCGATCGTCTGGGTAAGCTGTTGCGCGTAAAATGTCTTGCCGTAGGGGATGATTTCCTCTGTGGCTCTGGTCGGGAATGCGATGTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCGTCGATATCACCAGTACGCAAGCTTTGTGCGAAGGTGGCGTGCGTATCAGCAGCACCGATGTGCGTCAGGCGCTGGCTGAGGACAATCTGGCTCTGGCAGAAAGTTTa                                                                                                                                                                                 <  2:257340/250‑1 (MQ=255)                                                                tAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGCTTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTCCAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGTAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCCGAAAGTTTACTg                                                                                                                                                                              >  1:625406/1‑251 (MQ=255)                                                                    cgcgCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGc                                                                                                                                                                          <  2:481855/251‑1 (MQ=255)                                                                     gcgcCAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCATACAGTTTACTGGGGc                                                                                                                                                                          >  2:480159/1‑250 (MQ=255)                                                                     gcgcAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACGTGTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCa                                                                                                                                                                         <  2:629756/251‑1 (MQ=255)                                                                       gcAAAAGTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGTCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGGGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAAGCTGGCTCTGGCAGAAAGTTTACTGGGGCAcc                                                                                                                                                                       <  2:558764/251‑1 (MQ=255)                                                                        cAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCAccc                                                                                                                                                                      <  2:257506/251‑1 (MQ=255)                                                                         aaaaTTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCg                                                                                                                                                                     >  2:178537/1‑251 (MQ=255)                                                                                catcaGCGATCTTCTGGTGAAGCTTTTGCGCGTAATATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTTTTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGcc                                                                                                                                                               >  2:392694/1‑250 (MQ=255)                                                                                 atcaGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGCTGACAATCTGGCTCTGGCAGAAAGTTTATTGGGGCACCCGTTTGCCAt                                                                                                                                                             >  2:400393/1‑251 (MQ=255)                                                                                    aGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCAtctc                                                                                                                                                          >  1:515490/1‑251 (MQ=255)                                                                                        aTCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGCTGATTTACGCTTTGGCGCTGGTCGCGAACGCGATTTCTCATTCTTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGCACGCAAAATTTTTGAGAAGGTGGCGCGCGCACCAGCAGCACAGCTGCGCGTCAGGCGCTGGCGGATGACACTCCGCATCTGGCAGACAGTTTACTGGGGCACCCGTTTGCCGTCTCCgtg                                                                                                                                                      >  1:287760/1‑249 (MQ=255)                                                                                         tcttctGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGACAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGc                                                                                                                                                     >  1:62639/1‑251 (MQ=255)                                                                                            tctGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCgtg                                                                                                                                                  >  2:415924/1‑251 (MQ=255)                                                                                             ctGGTGACGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGCTTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTCTTACAGACAGCTGGCATGGAATACGGCTTCGATATAGCCAGAACGCAAACTTTTTGCGCGGGTGGCGTGCGTATCCGCAGAACCGCTGTCCGTCAGGCGCTGGCGGATCACAATCGGGCTCTGGCAGCAAGTTTAATGGGGCAACCGATTGCCATCTCCGGGCgtgt                                                                                                                                                 >  1:614913/1‑251 (MQ=255)                                                                                               ggTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCgtg                                                                                                                                                  <  2:493910/248‑1 (MQ=255)                                                                                                    agcttgtgctggTAAAATTGCTTGCCGTAGGTGATGATTGCCGCGTTGGCGCTGGTCGTGAAGTCGATTGCTTGTTATTACTGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCATCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCGTGCTCGGGCAGTAAGTTCACGAGGGCACCGGTTTGCCATCACCTTGCGTGTATTCCAc                                                                                                                                          <  2:111892/240‑1 (MQ=255)                                                                                                        tttGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACgg                                                                                                                                        <  1:550006/249‑1 (MQ=255)                                                                                                              cgTAAAATTTCTTGCCGTAGGTGATGATTTACGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAAt                                                                                                                                >  1:499869/1‑251 (MQ=255)                                                                                                                   aaTTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGGGATGAATTAggg                                                                                                                           >  1:337098/1‑251 (MQ=255)                                                                                                                     tttCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGgcg                                                                                                                         >  2:175870/1‑251 (MQ=255)                                                                                                                      ttCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGgcg                                                                                                                         <  1:43668/250‑1 (MQ=255)                                                                                                                      ttCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAACCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAAGgcgc                                                                                                                        >  1:248262/1‑251 (MQ=255)                                                                                                                           gCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGAGAGCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACta                                                                                                                    <  2:215740/250‑1 (MQ=255)                                                                                                                             cGTAGGTGATGATTTCCGCTTTGGCGATGGTCGTGAATGCGATTTCTTGTTATTACAGAATGCTGGCATGGAATACGGCTTCGATATCCCCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACACTCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACtat                                                                                                                   >  2:61151/1‑249 (MQ=255)                                                                                                                                 ggTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGttt                                                                                                             <  1:117932/251‑1 (MQ=255)                                                                                                                                   tgatgaTTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTCGTCCACGGTGATGACTTAGGGCGCAATATCGGtttac                                                                                                           >  2:29118/1‑249 (MQ=255)                                                                                                                                    gatgatTTCCGCTTTGGCGATGGTCGTGATGGCGATTTCTTGTTATTACAGACAGATGGCATGGAATACGGCTTCGCTATCACCAGTACGCAAACTTTTTGCGACGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCCGGCGCTGGCGGATGACAATCTGGATCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGCCCACGGTGTTGAATTAGGGCGCCATATAGGTTTccc                                                                                                          >  2:39631/1‑251 (MQ=255)                                                                                                                                                  tGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCCTGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGtt                                                                                             <  2:15034/250‑2 (MQ=255)                                                                                                                                                     cgcTGGACGTGAAGGCGATTTCTTGTTATTACAGCAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGAGGCGTGCGTATCAGCAGCACCGATGTGCGTCAGGCGCTGGCGGATGACCATCTGGCTCTGGCAGAACGTTTAATGGGGCACAACTATGCCATCTTCGGGAGCGTAGTCAACGGTGATGACTTAGGGCGCACTATCGGATTCCAGACGGCGAATGTACCGc                                                                                         >  1:114440/1‑251 (MQ=255)                                                                                                                                                                         tCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGGACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTcccc                                                                      <  1:582260/250‑1 (MQ=255)                                                                                                                                                                            tGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTTCGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTGCCCCgg                                                                    <  2:72649/249‑1 (MQ=255)                                                                                                                                                                                          gCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGGGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGc                                                    <  2:513795/251‑1 (MQ=255)                                                                                                                                                                                           cTGGCATGGAATACGGTTTCGATATCACCAGTACGCAAACTTTTTGCGAAGATGGAGTGCGTATCAGCAGCACCGCTGTGCGTAAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAACGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAGAGGGGTTTATGCg                                                   >  2:45342/1‑251 (MQ=255)                                                                                                                                                                                             ggCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCgg                                                  <  1:501254/250‑1 (MQ=255)                                                                                                                                                                                              gCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGTGGTTTATGCGGTa                                                <  1:414446/251‑1 (MQ=255)                                                                                                                                                                                               cATGGCATACGGCTTCGCTATCACCAGTACGCAACCTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGt                                                 >  2:587125/1‑249 (MQ=255)                                                                                                                                                                                                   gAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTCGAAGt                                           >  1:215071/1‑251 (MQ=255)                                                                                                                                                                                                   gAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGAAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAACGGGGTTTATGCGGTAGAAgg                                           >  2:354637/1‑250 (MQ=255)                                                                                                                                                                                                    aaTACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTg                                          <  1:603366/251‑1 (MQ=255)                                                                                                                                                                                                       aCGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTg                                       >  1:64003/1‑251 (MQ=255)                                                                                                                                                                                                        cGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTgg                                      <  1:250817/251‑1 (MQ=255)                                                                                                                                                                                                               gATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAAGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCgg                               <  1:208971/251‑1 (MQ=255)                                                                                                                                                                                                                 tatCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTg                             >  1:18065/1‑251 (MQ=255)                                                                                                                                                                                                                      ccAGTACGCAAACTGTTTGCGAAGGTGGCGTGCGTCGCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAg                        <  2:248208/251‑1 (MQ=255)                                                                                                                                                                                                                            cgaaaaCTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGTCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCTCTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGGTTATGCGGTATAAGGGCTGGTCCTCGGTGAAAAGCCGTTa                  <  2:365422/248‑1 (MQ=255)                                                                                                                                                                                                                             gCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGTTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTAc                 <  1:390417/251‑1 (MQ=255)                                                                                                                                                                                                                              cAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTAcc                >  1:260503/1‑251 (MQ=255)                                                                                                                                                                                                                               aaaCTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTAccc               >  1:426182/1‑251 (MQ=255)                                                                                                                                                                                                                                 aCTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCATCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGCTGAATTAGGGCGCACTATCGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCgg             >  2:277139/1‑251 (MQ=255)                                                                                                                                                                                                                                            aGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTAAAAAGCCGTTACCCGGCGTGGCAAACa  <  1:29118/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               ACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACA  >  NC_000913/21648‑22133

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