Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,067,300	C→T	G8R (GGA→AGA)	insH‑6 ←	CP4‑44 prophage; IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,067,300	0	C	T	100.0%	76.9 / NA	24	G8R (GGA→AGA)	insH‑6	CP4‑44 prophage; IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/14);  total (10/14)

TGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAGCGCCA  >  NC_000913/2067165‑2067436                                                                                                                                        |                                                                                                                                         tGACTTCCCCCCTGCTTTCCCATGGCAGAATCTGCCCCCCGCCGGACAAGAAAATCTCTTTTCCGGCCCGACGGCCCCTACCCCCGAATTCGCCGCCGGCGCAGACCAACTGATGACGCATGATGAACCCTGTTCCAtcg                                                                                                                                      >  2:211141/1‑138 (MQ=2)                                                                       aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGct                                                                >  2:642078/1‑140 (MQ=255)                                                                         ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGctct                                                              <  2:728896/140‑1 (MQ=255)                                                                         ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGctct                                                              <  1:611805/140‑1 (MQ=255)                                                                           cgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGctctct                                                            <  2:695117/140‑1 (MQ=255)                                                                               tACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGt                                                        <  2:655257/140‑1 (MQ=255)                                                                               tACTGCTGAATTCACTGTCGGCGAAGGTAAGCTGATGACTCATGATGACCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGt                                                        <  2:305342/140‑1 (MQ=255)                                                                                 ctgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTaa                                                      >  1:808776/1‑140 (MQ=255)                                                                                   gctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAgg                                                    <  2:127364/140‑1 (MQ=255)                                                                                     tgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCa                                                  <  1:642078/140‑1 (MQ=255)                                                                                     tgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCa                                                  <  2:22162/140‑1 (MQ=255)                                                                                         ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCGAGAACACGACAGTTTCGCTCTGTCGGTAAGGCACTGt                                              >  1:249035/1‑140 (MQ=17)                                                                                         ttCACTGTCGGCGAAGGTAAGTCGATGACTCCTGATGAATCCTGTTCTTTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGt                                              <  2:249035/140‑1 (MQ=18)                                                                                            aCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACGATTTCGCTCTCTCGGTAAGGCACTGTca                                            <  1:200260/139‑1 (MQ=255)                                                                                            aCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACGATTTCGCTCTCTCGGTAAGGCACTGTca                                            >  2:200260/1‑139 (MQ=255)                                                                                             cTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTcaca                                          <  2:662185/140‑1 (MQ=255)                                                                                                 cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCAt                                      >  2:57329/1‑140 (MQ=255)                                                                                                 cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCAt                                      <  1:604371/140‑1 (MQ=255)                                                                                                                   gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGaa                             >  2:631723/1‑131 (MQ=255)                                                                                                                   gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGaa                             <  1:631723/131‑1 (MQ=255)                                                                                                                   gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTg                    >  1:431143/1‑140 (MQ=255)                                                                                                                   gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTg                    >  1:172065/1‑140 (MQ=255)                                                                                                                                cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAg       >  2:823498/1‑140 (MQ=255)                                                                                                                                  cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAgcg     >  2:519269/1‑140 (MQ=255)                                                                                                                                     ttCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAGCGcca  <  2:793625/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         TGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAGCGCCA  >  NC_000913/2067165‑2067436

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