Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,601	C→T	A56A (GCG→GCA)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,601	0	C	T	96.0%	140.4 / NA	51	A56A (GCG→GCA)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/36);  minor base A (2/0);  total (14/37)
Fisher's exact test for biased strand distribution p-value = 7.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCC  >  NC_000913/3815357‑3815843                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   ccATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGAGTGGGCGCCATCGGGATCCCCTTTTAAACAGGGCCAAACAGCAGAACGAACTCAACTCCAGAATTCACCAAAGCTTc                                                                                                                                                                                                                                               >  1:377798/1‑250 (MQ=255)         ccgcgAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCGTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGAGGCGCCGCGAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTGTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAATGCTTCAGCGTaaaa                                                                                                                                                                                                                                      <  2:439511/250‑1 (MQ=255)                 cGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCt                                                                                                                                                                                                                               <  1:258214/250‑1 (MQ=255)                       tCCGTCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGt                                                                                                                                                                                                                         <  1:74729/250‑1 (MQ=255)                           gcgttgttcatctcagctcccAGCATGACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCTTGGTCTTTTGCCTCGTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGAGGCTCCGCCAGCGCCACGGCTGTAGTGTTCGCAATCGGGATCCCTTTATAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCGTCAGCGTAAATACTGCCTAACAGGGCCAg                                                                                                                                                                                                                    <  2:402945/232‑1 (MQ=255)                              gTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAggtc                                                                                                                                                                                                                 <  1:60538/251‑4 (MQ=255)                                    ccatcatCTACCAGCATTACGAGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTTAGGATGCTCCGCCAGCGCCACGGCAGTAGTGGGCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGccc                                                                                                                                                                                                           <  2:356393/250‑1 (MQ=255)                                           cTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGGTCCCTTTGTAAGCAGGGCCAAACAGCAGATCCAACTCAATGCCAGAATCCACCAATGCTTCAACGTAAAAAAGGCCTAACAAGGCCAGGGCGCGCCCGGGAtt                                                                                                                                                                                                     >  1:262865/1‑250 (MQ=255)                                             tcccgCATTACGCGTCCTTGTAACGCGCTACCAACCAGACTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTa                                                                                                                                                                                                    <  2:283666/245‑1 (MQ=255)                                             aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTAtt                                                                                                                                                                                                     <  1:456851/248‑1 (MQ=255)                                             aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTaa                                                                                                                                                                                                   <  1:274004/250‑1 (MQ=255)                                                               tttagcgcgCTCCCAACCAGAGTTCCGCCTTAAACGTGGTCTGTTGCCTATTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGTTGCGCCGGCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCGTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGCAATCCACAAATGCTACAGCGTAAAAACGGCCTAACAGTCCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTgaag                                                                                                                                                                                  <  2:358261/244‑1 (MQ=255)                                                               tGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGGAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCGTCAGCTTAAAAACGGCCTAACAGTGCCAGGTCGCGCCGGGTATTAAACAGCCCGGCGTTgaag                                                                                                                                                                                  <  2:435376/249‑1 (MQ=255)                                                                    cgcgCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGGAGGGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCCGAAACCACCCATGCTTCCGCGTAAAAACGGCCTAACAGTGCCAGGGCGCGCCCGGGATTAAACAGCCCGGCGTTGGAGAAATAgg                                                                                                                                                                           >  1:183821/1‑251 (MQ=255)                                                                        cTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTACAGCAGTACGGCAGGTAACGGTCAGGATGCTCCGCCAGCGCCCCGGCGGTAGTTGGCGCCATCCGGATCCCTTTTTAAGCAGGGGCAAAAAACAAATAGAAATCAAATCCAAAAACCACCAAAGCTTCAACGTCAAAAAGACCCAAAAATGCCAGGTCGCGCCCCGGATTAAACAACCCCGCCGTGAAGAAAtcagggct                                                                                                                                                                       >  1:248260/1‑244 (MQ=255)                                                                           ccaaccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCCGCGCCACGGCGGGAGTGGTCGCAATTGGGATCCCCTTGTAAGCAGGGCCAAACAGCCAAACGGAATCAAAGCCCGAATCCACCAATGCTTTAGCGGAAAAAAGGCCCAACAAGGCCAGGTCGCCCCCCGTATTAAACCACCCCGCGGTGAAGAAAAAGGGGGtttt                                                                                                                                                                    >  2:407299/1‑251 (MQ=255)                                                                              accaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGTCCTGATGCTCCGCCAGCGCCACGGCGGGAGGGGGCGCAATCCGGATCCCTTTGTAAGCAGGGGCAAAAAGCCGAACGGACTCAAAGCCAGAATCCCCCAATGGTTCAACGTAAAAACGGCCTAACAATGCCCGGTCGCCCCCCGGGTTAAACAGCCCCGCGtttgagaaaaaaggggtttttcc                                                                                                                                                                 >  1:74925/1‑230 (MQ=255)                                                                                caGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGcc                                                                                                                                                               <  1:192732/251‑1 (MQ=255)                                                                                      gccgccgcctccGTGGTCTTTTCCCTCTTTGCGGTTAAAGCAGTACTGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGTTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACACAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTgaa                                                                                                                                                                                   <  2:197342/215‑1 (MQ=255)                                                                                                             cctcttTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATTATGCACCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACGCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGGCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACttt                                                                                                                                  <  2:278356/249‑1 (MQ=255)                                                                                                                     cGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtg                                                                                                                          <  1:89547/251‑1 (MQ=255)                                                                                                                     cGGTTAAAGCAGTACGGCAGGTCAAGGTCAGGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtg                                                                                                                          <  2:345057/251‑1 (MQ=255)                                                                                                                         tAAAGGAGTACGTCAGGTCAAGGTGATGATGCTCCGCCAGCGCCACGGCGGTATTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCATATCGAACTCAATGCCAGAATCCACCAAGGCTTCAGCGTAAAAACGGCCTAACAGTTCCAGGTCGCGCCCGGTATTAAACAGCCCGGCTTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgctg                                                                                                                       <  1:276761/250‑2 (MQ=255)                                                                                                                                      caggtcaaggtcATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAtt                                                                                                         <  2:183821/251‑1 (MQ=255)                                                                                                                                              ggtcATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACGTTAACACCTGCTTGCTAAGCGCAAATTCAATAAAc                                                                                                 <  1:355585/251‑1 (MQ=255)                                                                                                                                               gtcatGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACt                                                                                                <  2:168690/247‑1 (MQ=255)                                                                                                                                                 catGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTgg                                                                                              <  1:396520/249‑1 (MQ=255)                                                                                                                                                       gCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTg                                                                                         <  1:83703/250‑1 (MQ=255)                                                                                                                                                                     acggCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCCCCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTGTCATGcc                                                                          <  2:274031/248‑1 (MQ=255)                                                                                                                                                                          ggtagaggtcgCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGCACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATACACTGGCGCTGATATGGTTTCATGCCTTCGc                                                                     >  2:141577/11‑251 (MQ=255)                                                                                                                                                                                    gCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCCATGCTTCAGCGTAAAAACGGCCTAACAGTGCAAGGTCGCGCCCGGTATTAAACAGCACGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTCAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTAAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCtt                                                           >  2:260652/1‑251 (MQ=255)                                                                                                                                                                                     cAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCtt                                                           >  1:481161/1‑250 (MQ=255)                                                                                                                                                                                           ggATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTc                                                    <  1:323898/251‑1 (MQ=255)                                                                                                                                                                                            gATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCt                                                   <  1:313753/251‑1 (MQ=255)                                                                                                                                                                                             aTCCCTTTGTAAGCAGGGCCAAACAGCAGATCTAACTCAATGCCAGAAGCCACCAATGTCTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCTCGCCCGGTATTAAACAGCCCGGCTTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTTAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCCTGCCTTCGCTCCTCGTCTGACTTTTCTa                                                  <  2:390435/251‑1 (MQ=255)                                                                                                                                                                                                   ttGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAc                                             <  1:290162/250‑1 (MQ=255)                                                                                                                                                                                                       aagcatgtccAAACAGCAGATCGAACTCAATGCAAGAATCCTCCAATGCTTCAGCGTAAAAACGGCCGAACAGTGCCAGGTCGCGCCCTGTATTAAACAGCCCGGCTTTGAAGAAATAGGGGTTTTTGCGCCCGGATTTCAGCGTAAACTCGACAAACTTTAACACCTGCGTGCTAAGCGCATATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCATCATCTTACTTTTCTACAGACaaaa                                         <  1:292605/242‑1 (MQ=255)                                                                                                                                                                                                             gtccAAACATCAGATCGCACTCAATGCCAGCATCACCCAATGCTTCATCGTAAAAACGGCCTAACAGTGCCATGTCGCCCCCGGTCTTTAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAATCTTTAACTCCTGCTTGCTAAGCGCCAATTCAATAATCGGGCGCTGATATGGTTTCATGACTTCGCTCCTCATCTTTCGTTTCTACAGACaaaaaaaaagc                                  <  2:140440/249‑4 (MQ=255)                                                                                                                                                                                                             ggCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGc                                  <  1:407299/251‑1 (MQ=255)                                                                                                                                                                                                              gCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCg                                 <  1:352363/251‑1 (MQ=255)                                                                                                                                                                                                                cAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGAc                               >  1:333823/1‑251 (MQ=255)                                                                                                                                                                                                                      gcagATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatca                         <  2:349876/251‑1 (MQ=255)                                                                                                                                                                                                                         gATCGAACTCAATGCCAGAATCCACCAATGCNTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc                      <  2:228067/251‑1 (MQ=255)                                                                                                                                                                                                                              aaCTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCCCCtt                 >  2:3097/1‑251 (MQ=255)                                                                                                                                                                                                                               aCTCAATGCCAGAAGCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACGTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTa                <  2:171302/251‑1 (MQ=255)                                                                                                                                                                                                                                      gccagAATCCACCAATGCTTCAGCGTAAGAACGACCTAACATTGCCAAGTCGCGCCCGGTATTAAACAGCCCGGCGTGGAAAAAAGAGGGGCTTTGGCGCCCGGATTTCAGCGTAAACTAGCCAAACTTTAACACCGGCTGACTAAGCGCAAATTCAATAAACTGACGCTGATATGGTGTCATGCCTTCGCTACTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCATTAAAGATc          >  1:250196/5‑250 (MQ=255)                                                                                                                                                                                                                                      gccagAATCCACCAATGCTTCAGCGTAAAAANGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCTTAGAAGAAATAGGGGCTGTTGCGCCCGGATTTCAGTGTAAACGCGCCAAACGTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACAGGCGCTGATATGGGTTCATGCCTTCGCTCCGCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa         <  2:248260/247‑1 (MQ=255)                                                                                                                                                                                                                                         agAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGtt      <  1:399951/250‑1 (MQ=255)                                                                                                                                                                                                                                              ccaccaATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc  <  1:149456/250‑1 (MQ=255)                                                                                                                                                                                                                                               caccaATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc  >  1:463088/1‑249 (MQ=255)                                                                                                                                                                                                                                               caccaATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc  >  1:2336/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   CCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCC  >  NC_000913/3815357‑3815843

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