Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,589	G→A	S60S (TCC→TCT)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,589	0	G	A	95.6%	146.1 / NA	47	S60S (TCC→TCT)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  major base A (9/34);  minor base C (2/0);  total (11/36)
Fisher's exact test for biased strand distribution p-value = 5.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTA  >  NC_000913/3815348‑3815829                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 gAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCGCTTTGCGGTTAAAGCAGTACGTCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCTGTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATTCCAGAATcctcc                                                                                                                                                                                                                                         <  2:257230/251‑4 (MQ=255)          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)                                                                                   accaaccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCaatg                                                                                                                                                                                                                                     <  1:472985/173‑3 (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)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 GAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTA  >  NC_000913/3815348‑3815829

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