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%	148.7 / NA	48	E162E (GAG→GAA)	yceF	m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/26);  minor base G (0/1);  total (21/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.78e-01

GTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCC  >  NC_000913/1145871‑1146344                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           gttTGAAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAAGGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCGGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGAAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGt                                                                                                                                                                                                                                 <  2:116965/251‑1 (MQ=255)   ttGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTa                                                                                                                                                                                                                                <  1:458297/250‑1 (MQ=255)    tGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGCCGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAAt                                                                                                                                                                                                                              >  2:65408/1‑251 (MQ=255)                gTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTc                                                                                                                                                                                                                  >  1:413717/1‑251 (MQ=255)                            aataaaagaaaaaaaaccagaaaaggggCGGTTGGCCTGATAAGGCGTTAAAGCAGAACCAGGCGGCCAGCACCGATGGCCGGAGGCGACGTAACCACAACCGGCGAAAAAAACTATCCCATCAGAGGGTTTGTCCCATCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAAGCCTTCACTCTTAAAGCt                                                                                                                                                                                                      <  2:351815/226‑1 (MQ=255)                                 aaGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCcg                                                                                                                                                                                                 >  1:25296/1‑251 (MQ=255)                                    cATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACCGCGTAATGCCAAAACCTTCACTCTTAAAGCTACCCgcgc                                                                                                                                                                                              >  2:264735/1‑251 (MQ=255)                                                                     aaGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCAGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCCCGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTTCAGGGGATTCTCTTTACGCCCGGAATTAt                                                                                                                                                             >  2:190674/1‑251 (MQ=255)                                                                      aGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGAAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCTCTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATc                                                                                                                                                            <  1:247971/251‑1 (MQ=255)                                                                       ggCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATc                                                                                                                                                            >  1:31514/1‑250 (MQ=255)                                                                                   gcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCg                                                                                                                                               >  1:295238/1‑251 (MQ=255)                                                                                             gccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGAc                                                                                                                                     <  1:209026/251‑1 (MQ=255)                                                                                                         tGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGCCGGAAATGGACGt                                                                                                                         >  2:2759/1‑251 (MQ=255)                                                                                                                   gACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTAAAAAGGCTc                                                                                                               >  2:263533/1‑251 (MQ=255)                                                                                                                        aaCCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACt                                                                                                           <  1:175393/250‑1 (MQ=255)                                                                                                                              agccGGCGAAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGttt                                                                                                    <  2:12275/249‑1 (MQ=255)                                                                                                                                  ggCGCAAAAAACTTTCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCTTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAAGCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAg                                                                                                <  1:454788/251‑1 (MQ=255)                                                                                                                                  ggCGCAAAAAACTATCCCATCAGCGGGTTTTCCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAg                                                                                                <  2:31514/251‑1 (MQ=255)                                                                                                                                    cgcAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCa                                                                                                                                                                                        <  2:308701/161‑1 (MQ=255)                                                                                                                                    cgcAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCa                                                                                                                                                                                        >  1:308701/1‑161 (MQ=255)                                                                                                                                    cgcAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGAt                                                                                              <  1:461848/251‑1 (MQ=255)                                                                                                                                    cgcAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGAt                                                                                              <  1:83585/251‑1 (MQ=255)                                                                                                                                        aaaaaCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGccc                                                                                          >  1:371579/1‑251 (MQ=255)                                                                                                                                         aaaaCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTACACGTTCAAACAGCGTAATGCCAAATCCTTCCCTCTTAAAGCTACCCGAGCAGTGCCGGGGATGATCTTTACGCACGTACTTATCAATCTCCGCCTCCCTCAGGTGCCGGAAATGGACGTCAAAAGGCACCAATTCTGTTTGCAAATGccc                                                                                          >  2:197744/1‑250 (MQ=255)                                                                                                                                           aaCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCCGATGCCCAtt                                                                                       >  2:115439/1‑251 (MQ=255)                                                                                                                                            aCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTc                                                                                      <  1:159712/251‑1 (MQ=255)                                                                                                                                                    aTCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAAtt                                                                              <  1:65408/251‑1 (MQ=255)                                                                                                                                                        gCGGGTGTTTCCCCTCCCGACGTAACAGCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGGTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGTTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGATGTCAAAAGGCTCCACTTCTGTTTGCAGATGCGCATTCGCCGAATTAaac                                                                          <  1:402157/251‑1 (MQ=255)                                                                                                                                                             tttttCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCAc                                                                                                                                                                     <  1:394642/155‑1 (MQ=255)                                                                                                                                                             tttttCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCAc                                                                                                                                                                     >  2:394642/1‑155 (MQ=255)                                                                                                                                                               tttCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGccc                                                                   >  1:211374/1‑250 (MQ=255)                                                                                                                                                                 tcccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGccag                                                                  <  2:71689/250‑1 (MQ=255)                                                                                                                                                                  cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGccagtc                                                                <  2:371579/251‑1 (MQ=255)                                                                                                                                                                  cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGccagt                                                                 <  1:263533/250‑1 (MQ=255)                                                                                                                                                                       cccGACGGAACCTCTGACACCGCGGGATTAGCGGTAAACCAACCAGTGTGTTAGGATCCCGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCATCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTGTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATTGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACTCCAGTCCAGta                                                           <  2:108470/251‑1 (MQ=255)                                                                                                                                                                         cGACGTAACATCTGACACAGCGCGATTAGCGGTAACCCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTAGCCGAATTAAACAACGCCAGTCCAGtat                                                          >  2:447132/1‑250 (MQ=255)                                                                                                                                                                                   tCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGt                                                                                                                         <  1:361491/177‑1 (MQ=255)                                                                                                                                                                                   tCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGt                                                                                                                         >  2:361491/1‑177 (MQ=255)                                                                                                                                                                                          caGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGCTCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTACACAACGCCAGTCCAGTATAGAAGGAGCCGATATTGc                                        >  2:288850/1‑251 (MQ=255)                                                                                                                                                                                             cgcgATTAGCGGTAAACCAACCAGCGTTTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGAGGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATGCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGc                                     <  2:384785/251‑1 (MQ=255)                                                                                                                                                                                                 aTTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTg                                                                                                   >  2:107542/1‑185 (MQ=255)                                                                                                                                                                                                 aTTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTg                                                                                                   <  1:107542/185‑1 (MQ=255)                                                                                                                                                                                                    agctgTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCtt                               <  1:299619/246‑1 (MQ=255)                                                                                                                                                                                                                   aGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACg               >  1:181312/1‑251 (MQ=255)                                                                                                                                                                                                                       tgtTAGGATCACGCCCCTCTAAACGTTCAAACAGAGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGcc           >  1:114470/1‑250 (MQ=255)                                                                                                                                                                                                                           aGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATATAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCAttt        <  1:415426/250‑1 (MQ=255)                                                                                                                                                                                                                            ggctccgcccccTCTAACCGTTAAAACAGCGTAAGGCCACATCCGTCACTCTTAAAGCTACCCGCGCAGGGCAGGGGAAGCGCTTTACGCACGTAATTATCAATCGCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCAattt       <  2:340161/242‑1 (MQ=255)                                                                                                                                                                                                                                cACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTcc  <  1:2759/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           GTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCC  >  NC_000913/1145871‑1146344

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