Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,146,110	C→T	E162E (GAG→GAA)	yceF ←	m(7)GTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,146,110	0	C	T	97.9%	143.7 / NA	47	E162E (GAG→GAA)	yceF	m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/26);  minor base G (0/1);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCC  >  NC_000913/1145868‑1146344                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           gttgttTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACac                                                                                                                                                                                                                                    >  1:385625/1‑250 (MQ=255)            aaccTGCGTACCGCGCACTAAGACACGCATATAAACCTGAAAATGGGCGGTGGGCCGGCTACGGCGTTTCCGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGAGACGTAACCATATCCGGCGCAAAAAACTATCCCTTCAGCGGGTTTTTCCCTTCCCGACGTAACTTCTGACACAGCGCTATTAGCGTTAACCCAACCAGCGTGTTAGGATCCCGCCCCTCTTAACGTTCATACAGCTTAATGctaa                                                                                                                                                                                                                         <  2:462599/248‑4 (MQ=255)                  cGTAACGCACGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAACGCCAAATCCtt                                                                                                                                                                                                                   >  2:321143/1‑251 (MQ=255)                             atAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAg                                                                                                                                                                                                        <  1:172032/250‑1 (MQ=255)                               aaTACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCt                                                                                                                                                                                                      >  1:403024/1‑251 (MQ=255)                                aTACAATCATAGAAACCTGATAATGGGCGTTTGGACTGATAAGGCGTTTGCACCGCATCAGGCCGCCAGCCCCGATTGCCGGATGCGACGTAACCACAGCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTTTGATGAGCACGCCCCTCGAAACTTTCAAACAGCGTAATGCCGAATTCGTCCCTTTTAAAGCTa                                                                                                                                                                                                     <  2:612213/251‑1 (MQ=255)                                   aaaGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTAccc                                                                                                                                                                                                  <  2:403024/250‑1 (MQ=255)                                    aaGCTTATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTGTTTCCCTTCCCGACGTAACATCTGACACAGCGCGAGTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCGAAACGGTCAAACAGCGTAATGCCAAATCCTGCACTCTTAAAGCTACCcg                                                                                                                                                                                                 <  2:409696/251‑1 (MQ=255)                                          ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCTCGCAGt                                                                                                                                                                                           <  1:27452/251‑1 (MQ=255)                                          ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAg                                                                                                                                                                                            <  2:139231/250‑1 (MQ=255)                                                  gATAATGGGCTGTTGGCCTGATACGGCGTTTCCGCCGCAGCAGGCCGCCAGCACCTATTGCCGGATGCGACGTAACCCCATCCGGCGCATAAAACTATCCCACCAGCGGGTTTTCCCCTTCCCGACGTAACATCTGACACAGCGCGAGTAGCGGTAAACCAACCAGCGTGTTAGGATCACTCCCTTCTAAACGTTCAAACAGCTTAATGCCAAATCCTGCACTCGTAAAGCTACCCTCGCAGTTCAGGGGa                                                                                                                                                                                   <  2:385295/251‑1 (MQ=255)                                                   aTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGAt                                                                                                                                                                                  <  1:139986/251‑1 (MQ=255)                                                          gCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCttt                                                                                                                                                                           <  1:396719/251‑1 (MQ=255)                                                                ggCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCAc                                                                                                                                                                     >  1:46806/1‑251 (MQ=255)                                                                 gCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCCGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGCTTAGCTGGAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACATTCAAACAGCGTAATGCCACATCCTTCACTCTTAAAGCTACCAGCGCAGTGCAGGGGATGCTCTCTACGCACg                                                                                                                                                                    >  1:12342/1‑251 (MQ=255)                                                                        aaGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTAt                                                                                                                                                             >  2:338195/1‑251 (MQ=255)                                                                        aaGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCAATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTAAGCACGTAATTa                                                                                                                                                              <  1:486104/250‑1 (MQ=255)                                                                           gCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGATTTTCCCTTCCCGACGTAACATCTGACACAGCGCGCTTACCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCCGGGGATGCTCTTTACGCACGTAATTATCa                                                                                                                                                           >  1:604781/1‑250 (MQ=255)                                                                              tttccgccgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCgcgc                                                                                                                                                                                              <  2:429470/208‑1 (MQ=255)                                                                              tttACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCgcgc                                                                                                                                                                                              >  1:429470/1‑212 (MQ=255)                                                                                tACGCCGCATCAGGCCGCCAGCACCAATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAACACGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGCTGCTCTTTACGCACGTAATTATCACTCTcc                                                                                                                                                     >  2:40119/1‑251 (MQ=255)                                                                                 aCGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTcc                                                                                                                                                     <  1:128400/250‑1 (MQ=255)                                                                                       cATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGc                                                                                                                                              <  1:303965/251‑1 (MQ=255)                                                                                            ggCCGCCAGCACCGATTGCCGGATGCGACGGAACCACATCCGGCGCAAAAAACTATACCATCAGCGGGTTGTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACATCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAgg                                                                                                                                         <  2:304854/251‑1 (MQ=255)                                                                                                         gATTGCCGGATGCGACTTAACCACATCCGGCGCAAAAAACTATCCCAGCAGCGGGTTTTTCCCGTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATGATCAATCACCGCCTCGCTCAGGTGACGGATATGGa                                                                                                                            <  2:607337/251‑1 (MQ=255)                                                                                                          aTTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTGCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCATCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTTTCAATCTCCGCCGCGTTCATGTGACGGAAATGGAc                                                                                                                           <  2:269460/251‑1 (MQ=255)                                                                                                             gCCGGATGCGACGTAACCCCATCCGGCGCAAAAAACTATCCCATCAGCGGGTGTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGAGGCTCTTTACGCACGTAATGATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTc                                                                                                                        <  2:561396/251‑1 (MQ=255)                                                                                                                      gACGTAACCACATACGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTAAGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTc                                                                                                               >  2:575116/1‑251 (MQ=255)                                                                                                                         gTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTc                                                                                                                                                >  1:124169/1‑215 (MQ=255)                                                                                                                         gTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTc                                                                                                                                                <  2:124169/215‑1 (MQ=255)                                                                                                                                       cgcAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGAt                                                                                              <  1:2746/251‑1 (MQ=255)                                                                                                                                                          aGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTaaa                                                                           <  1:178323/251‑1 (MQ=255)                                                                                                                                                               gTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGATAAACCAACCAGCGTGTTAGGACCACGCCCCTCTAAACGTTCAAACAGCGTAATCCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGAAGGGGATGCTCTTTAAGCACGTAATTATCAATATCCGCCTCGCTCAGGGGACGGAAATGGCCGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAaacac                                                                        >  1:630456/1‑248 (MQ=255)                                                                                                                                                                tttttCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGACTTAAACAACGc                                                                     >  2:245831/1‑251 (MQ=255)                                                                                                                                                                   ttcccctcccGACTTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTGTGCAGATGCCCATTCGCCGAATTAAACAACGcca                                                                   <  2:385625/250‑1 (MQ=255)                                                                                                                                                                     cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTa                                                                                                                                                                          >  1:286249/1‑145 (MQ=255)                                                                                                                                                                     cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCGCTTTa                                                                                                                                                                          <  2:286249/145‑1 (MQ=255)                                                                                                                                                                       cttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCTCTCAGGTGACGGAAATGGCCGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGCATTAAACAACGccagtcc                                                               >  2:207052/1‑250 (MQ=255)                                                                                                                                                                                     aTCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGAGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGATCCACTTCTGTTTGCGGATGCCCATTCGACGACTTAACCAACGCCAGTCCAGTATAGAAGGGGAc                                                >  2:529805/1‑251 (MQ=255)                                                                                                                                                                                          acacAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGatat                                           <  1:575116/251‑1 (MQ=255)                                                                                                                                                                                                          ggTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCg                           >  1:565591/1‑251 (MQ=255)                                                                                                                                                                                                               accaaccaGCGTGTTAGGATCCCGCCCTTCTATACGTTCAAACAGCGTCATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGCGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAACAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAAtt                      >  2:317439/1‑251 (MQ=255)                                                                                                                                                                                                                   accaGCGTGTTAGGATCCCGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTAACGATATTGCCGCTGGCTTTGCGTAATTGc                    >  1:176194/1‑249 (MQ=255)                                                                                                                                                                                                                       gCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACgg              <  1:224010/251‑1 (MQ=255)                                                                                                                                                                                                                       gCGTGGTATGATTACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACgg              <  1:352830/251‑1 (MQ=255)                                                                                                                                                                                                                                   cACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTcc  <  1:245831/251‑1 (MQ=255)                                                                                                                                                                                                                                   cACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGAAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTcc  >  1:55527/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           GTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCC  >  NC_000913/1145868‑1146344

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