Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,503,361	C→T	G24E (GGG→GAG)	insG ←	IS4 family putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,503,361	0	C	T	100.0%	89.3 / NA	27	G24E (GGG→GAG)	insG	IS4 family putative transposase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/15);  total (12/15)

CCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCCCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAAC  >  NC_000913/4503133‑4503588                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    cccAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACg                                                                                                                                                                                                               >  1:140891/1‑251 (MQ=255)     ggCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGc                                                                                                                                                                                                              >  1:37200/1‑248 (MQ=255)             cGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGt                                                                                                                                                                                                    >  2:60953/1‑250 (MQ=255)                       tAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATAcc                                                                                                                                                                                         >  2:76494/1‑251 (MQ=255)                           ggCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGAt                                                                                                                                                                                     <  1:98112/251‑1 (MQ=255)                                                     tGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATaaa                                                                                                                                                            <  1:76494/250‑1 (MQ=255)                                                             agcaTGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCgg                                                                                                                                                   <  1:246188/251‑1 (MQ=255)                                                                   atgTCCAGGCGATTCACAAACTGTTGAAGAGGTTCATTACGCTCAAGCGCCATGACAACAATACAACAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTACGTGGTCGAGGTAATCCTCCAGAGAATTCAGTGGGTTACGCAGAAAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACaa                                                                                                                                             <  2:227529/251‑1 (MQ=255)                                                                       ccAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGTAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGa                                                                                                                                         <  2:184572/251‑1 (MQ=255)                                                                        cAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAg                                                                                                                                        <  2:6753/251‑1 (MQ=255)                                                                                                    tCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTg                                                                                                             <  2:74961/250‑1 (MQ=255)                                                                                                             tCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCa                                                                                                    >  1:75400/1‑250 (MQ=255)                                                                                                                aGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTg                                                                                                <  1:132574/251‑1 (MQ=255)                                                                                                                    ccATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATcc                                                                                            >  2:196706/1‑251 (MQ=255)                                                                                                                                  cacCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCa                                                                              <  1:173613/251‑1 (MQ=255)                                                                                                                                                 tCGAGGGGAAGACGGCGCTTGCGTGGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATTTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCt                                                                <  2:91594/250‑1 (MQ=255)                                                                                                                                                          aGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAAc                                                                                                                                                                                                 >  1:71551/1‑112 (MQ=255)                                                                                                                                                          aGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAAc                                                                                                                                                                                                 <  2:71551/112‑1 (MQ=255)                                                                                                                                                          aGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCAGCATTAATCTTGGTTCTGGTGTTTGTa                                                      <  1:16187/251‑1 (MQ=255)                                                                                                                                                                         aGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGCAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGACACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCt                                        >  1:207066/1‑250 (MQ=255)                                                                                                                                                                         aGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCt                                        >  1:91693/1‑250 (MQ=255)                                                                                                                                                                               aCAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACaaaa                                  <  2:28368/250‑1 (MQ=255)                                                                                                                                                                                     ccTGATTCGGCAAGGCAACGAGAGATGAGTTAGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCt                           <  1:177751/251‑1 (MQ=255)                                                                                                                                                                                            cGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTtg                    >  2:114170/1‑251 (MQ=255)                                                                                                                                                                                                aaggcaaCGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGAc                >  1:90896/1‑251 (MQ=255)                                                                                                                                                                                                 aggcaaCGAGAGATGAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGAc                >  2:231613/1‑250 (MQ=255)                                                                                                                                                                                                               gAGTTCGGGGTCGAGGTAATCCTCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTTCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAAAAAACTATCAGCGACAAAAATATGCTCAATTTGTGACATCAGTAACAAAAc  <  2:193175/250‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    CCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATACACCAGACCATCATTTCGAGGGGAAGACGGCGCTTGCGTAGCGTTACAGTACCTGATTCGGCAAGGCAACGAGAGATGAGTTCGGGGTCGAGGTAATCCCCCAGAGAAGTCAGTGGGTTACGCAGAGAATCGTAACGGGATACCAGATCAAGAGCCTGTCCAATGTGCATAAAAAAATCCGGAAACAAGTGAGCGTTTCCGGATTCTTACACAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAAC  >  NC_000913/4503133‑4503588

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