Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	440,662	C→A	9.9%	S147I (AGC→ATC)	ispA ←	geranyltranstransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	440,662	0	C	A	9.9%	148.8 / 3.2	56	S147I (AGC→ATC)	ispA	geranyltranstransferase
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base A (3/3);  total (29/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.82e-01

GTTGCTGGTCGGCACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCA  >  NC_000913/440369‑440962                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                             gTTGCTGGTCGGCTCCCTGGCGTTTTCCCATCGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTGGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAg                                                                                                                                                                                                                                                                                                       <  1:210431/301‑1 (MQ=255) gTTGCTGGTCGGCACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCAAGCGCCAg                                                                                                                                                                                                                                                                                                       >  2:41051/1‑301 (MQ=255)    tctcttccgatctcTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtgg                                                                                                                                                                                                                                                                                                    <  2:174521/288‑3 (MQ=255)     ttGGTCGGCACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggt                                                                                                                                                                                                                                                                                                   <  2:597387/300‑4 (MQ=255)       gttCGGCACCCTGGCGTTTTCCCAACGTTGCAGTAGTTCCCACCACATCCAGGATGGCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggtat                                                                                                                                                                                                                                                                                                 <  1:40998/299‑6 (MQ=255)       gtgcGGCACCCTGGCGTTTTCCCAACGTTTCATTATCTCCCACCCCATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGa                                                                                                                                                                                                                                                                                                 <  1:445156/298‑1 (MQ=255)       ggTCGGCACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGa                                                                                                                                                                                                                                                                                                 <  2:568112/301‑1 (MQ=255)         tCGGCACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGAGATCggaaaagcacaca                                                                                                                                                                                                                                                                                               >  2:131124/1‑289 (MQ=255)               cctttgCGTTTTCCCACCGTTGCAGTTTCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCGGATGCTGTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCg                                                                                                                                                                                                                                                                                         <  1:154718/296‑1 (MQ=255)                 cTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAACACCAGTGGAGaacggaagagg                                                                                                                                                                                                                                                                                       >  1:174324/1‑292 (MQ=255)                 cTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTACGCGTCTAAATCTAATGCCTGACAACCGCACAATCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGaa                                                                                                                                                                                                                                                                                       >  2:220826/1‑301 (MQ=255)                           cccTACGTTGCAGTATCTCCCACCACATCCAGGATGTTATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTcgc                                                                                                                                                                                                                                                                             <  2:292813/301‑1 (MQ=255)                            ccAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGACGGTTTTATGACGCTGCATACGCTCAAGCGCCGACAGAGGTACGGGTATGCCTTCCGCGTCGACATCTAATGCCTGACCACCGCACAATCCGGCAATACCAATGGAGCGCGACAGATCAAAACACAACGaaacaacgtaccg                                                                                                                                                                                                                                                                            >  1:403297/1‑291 (MQ=255)                                 gTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCg                                                                                                                                                                                                                                                                       <  2:400598/301‑1 (MQ=255)                                      aGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAGGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACt                                                                                                                                                                                                                                                                  >  2:301424/1‑301 (MQ=255)                                                  cacaTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCTCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCg                                                                                                                                                                                                                                                      >  2:511977/1‑301 (MQ=255)                                                    caTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATATCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGACTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAACTCTAATGCCTGACCACCGCACATTCCGGCACTACCACTGGCGCTCGCCAGTTCACAAATCATCGAAATTCTGTCGCGGTCCGACACATCCGGCataacagc                                                                                                                                                                                                                                                    >  1:526481/1‑296 (MQ=255)                                                                             gAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACCTTCCGGCAATACCACTGGCGCCCGCCAGTTCAGACATCATCGAAATTCTCTCGCGGTCAGCCACTTCAGGCATATCGGACACCCTTAAAATCAAAAACGCAAGc                                                                                                                                                                                                                           >  1:65672/1‑301 (MQ=255)                                                                                ggcaaggcCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCAAGCGtt                                                                                                                                                                                                                        >  2:225514/1‑301 (MQ=255)                                                                                              tctctGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGTGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACAATTCCGGCATAACGGCATCGCTTAAAATCGAAAACGCAAGCGATTGGAAAGCGTcccc                                                                                                                                                                                                          >  1:311122/1‑301 (MQ=255)                                                                                                   gCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCCCCaccag                                                                                                                                                                                                     >  1:77763/1‑297 (MQ=255)                                                                                                      tACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGaa                                                                                                                                                                                                  >  2:155572/1‑301 (MQ=255)                                                                                                        ggagttcagacgtgtgctcttccgatctCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATc                                                                                                                                                                                                <  2:5862/274‑1 (MQ=255)                                                                                                                         gagaACGACTTCCTTTATCTTCGGCTTTTAATGCACCAAGGCGAACGGCGGCGCGAAGTAATGCGCCGGTGTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCAGCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCaaa                                                                                                                                                                               <  1:511305/301‑1 (MQ=255)                                                                                                                           cgacttggTCCTTTATCTCCGGCGCGTAAGGCACCAAGGAGAACGGCGGCGGGAATCAATGCGACGGTTTTATGACGATGAATACGCTCAAGGGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCGTCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCTAAACt                                                                                                                                                                             <  1:302008/294‑1 (MQ=255)                                                                                                                              cgacgTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTca                                                                                                                                                                          <  2:77853/301‑1 (MQ=255)                                                                                                                                 cgTCCTTTCTCTCCGTCGGTTAATGAACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGTTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACAt                                                                                                                                                                       <  1:47130/301‑1 (MQ=255)                                                                                                                                    ccTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCCGCATCGCTTAACATCGAGAACGCCAGCGTTTGTAAAGCGTCGCAAGCGAGaaacacaccggacaagaaccgcggagggaaaa                                                                                                                                                                    >  1:5862/1‑271 (MQ=255)                                                                                                                                                   cTTAATGCACCATGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGACTCGCGTTTGCTTCCCCAAACTTCCACTTGCAGGTTGGCACACcgc                                                                                                                                                     >  1:194765/1‑301 (MQ=255)                                                                                                                                                               gtgcgtgtggcggcGCGTATTAGTGCGCGTGGTTTATTGCGATGAATTCGTGTAATCGCGTCCAGAGTGACTTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACCTTCCGGCAATACCACGGGCGTTCGTCAGTTCAGAAATCAGCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCAGCGCTTAGAATCGAGAACGCCACCGTTTGTAAAGCGTCGCCCGCGAAAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGt                                                                                                                                         <  1:220583/293‑1 (MQ=255)                                                                                                                                                               aGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACAGCGACGCAGATCGt                                                                                                                                         >  2:201651/1‑301 (MQ=255)                                                                                                                                                                  cGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCCCGTTTGCTTCGCCAAACTTCACATGGCAGATTGGCAAACCGCCACCCACAACGtcac                                                                                                                                      >  1:295265/1‑300 (MQ=255)                                                                                                                                                                     aCGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGCGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTTGCAAACCGCGACACCAAACGtcatcaa                                                                                                                                   >  2:336775/1‑300 (MQ=255)                                                                                                                                                                                attCATTGCGCTGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGt                                                                                                                        <  2:578357/299‑1 (MQ=255)                                                                                                                                                                                aaTCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTTGAAAACCCCGACGCCAACAGCAATCACCCATTGCCGGt                                                                                                                        >  1:582861/1‑301 (MQ=255)                                                                                                                                                                                aaTCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGt                                                                                                                        <  2:195008/301‑1 (MQ=255)                                                                                                                                                                                   cAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTaaa                                                                                                                     <  2:583585/301‑1 (MQ=255)                                                                                                                                                                                   cAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTaaa                                                                                                                     <  2:527153/301‑1 (MQ=255)                                                                                                                                                                                   aattGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGGTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGACCGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTaaa                                                                                                                     <  1:225271/298‑1 (MQ=255)                                                                                                                                                                                                                 cttaaGCGCGTTCAGGGTTACGTGTTTGCCGTCCGCGTCTAAATCTAATGCCTGACCACCGGACATTCCGGCAATACCACTGGCGCTGGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGTGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATa                                                                                       <  1:281557/298‑1 (MQ=255)                                                                                                                                                                                                                 cTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGGCGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGAGTGGATa                                                                                       >  2:602963/1‑301 (MQ=255)                                                                                                                                                                                                                      gcgcGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAAATCGTCATCATCCATTGCCGGTAAAACATCATGAATTAATGAGTAAGCGAAGATACACta                                                                                  >  2:495201/1‑300 (MQ=255)                                                                                                                                                                                                                       cgcgTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGCCCTCCTACATTGCCGGTAAACCAACATCAATTAATAAGTAACCGTGAATACACTCa                                                                                 >  1:93894/1‑301 (MQ=255)                                                                                                                                                                                                                              agtggtgtgtttttGCCTTCCGCGTCTAAAGCTTATGCCTGACCACCGCACTTTCCGGCAATACCACTGGCGCTCGACAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTACAGATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGAGACACTCAGcggcgg                                                                          <  1:155421/290‑1 (MQ=255)                                                                                                                                                                                                                              agatGTACGTTTTTGCGTTCCGTGTTTAAATTTAATTCTGGTCCACCTCACATATCTGCTATACCACTGTCGCGCGCCAGGTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTGGCCAGCGAAAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAAcggcgg                                                                          <  1:494554/301‑1 (MQ=255)                                                                                                                                                                                                                                  gTATGTGTTTGCCGTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGc                                                                      <  1:201408/301‑1 (MQ=255)                                                                                                                                                                                                                                            gCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTcc                                                            >  2:13954/1‑301 (MQ=255)                                                                                                                                                                                                                                                       cTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTtgt                                                 >  2:553171/1‑301 (MQ=255)                                                                                                                                                                                                                                                           aTCTAATTGCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGGTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTa                                             <  1:301020/301‑1 (MQ=255)                                                                                                                                                                                                                                                               aaTGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGc                                         <  2:479374/301‑1 (MQ=255)                                                                                                                                                                                                                                                                  gCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGATTGTGCTAACGCCGa                                      >  2:341140/1‑301 (MQ=255)                                                                                                                                                                                                                                                                            cGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTCAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCACACGTCACATGGCAGGTTGGCAACCCGCCACGCAACTCGTCATCATCCATTGCCCGTCAAGCCTCATGAATACATGAGTACACGTCGCAACAACCAACGCAGGCAGACGGACCGTACCACGAGTTTGTGCTAACCCCCAACAAATCAcc                            >  1:501615/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                cgttCTGGCAATACCACTGGCGCTTGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGACCGGTg                        <  1:126781/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                caTTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAACGTGGATACCCTCAACGGCGGCAGAGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGAACGATg                        >  2:72775/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                    ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACCCTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGACCGGTGGCAt                    >  1:111124/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                      ggCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGACCGGTgagatcg                  >  2:126905/1‑295 (MQ=255)                                                                                                                                                                                                                                                                                                      cTCGCCGGTTTAGAATTGTTCGAAAGTCTGTCGCGGTCCGCCACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGGCGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCa  <  1:602202/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                             GTTGCTGGTCGGCACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAACCGCGACGCAGATCGTCATCATCCATTGCCGGTAAATCATCATGAATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCA  >  NC_000913/440369‑440962

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