Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,431,830	G→A	A928T (GCT→ACT)	stfR →	Rac prophage; putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,431,830	0	G	A	100.0%	121.8 / NA	37	A928T (GCT→ACT)	stfR	Rac prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/15);  total (22/15)

GGTTATGGTTTTGAGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGC  >  NC_000913/1431695‑1431896                                                                                                                                        |                                                                   ggTTATGGTTTTGAGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACttt                                                                >  2:777659/1‑140 (MQ=255)      tGGTTTTGAGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACa                                                           >  1:277577/1‑140 (MQ=255)       ggTTTTGAGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACaa                                                          <  2:223929/140‑1 (MQ=255)            tGAGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTg                                                     >  2:40424/1‑140 (MQ=255)            tGAGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTg                                                     >  2:1258386/1‑140 (MQ=255)             gAGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGc                                                    >  2:571326/1‑140 (MQ=255)             gAGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGc                                                    >  1:508346/1‑140 (MQ=255)              aGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCt                                                   <  2:426946/140‑1 (MQ=255)              aGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCt                                                   >  1:37459/1‑140 (MQ=255)                   gACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTAccc                                              >  2:340021/1‑140 (MQ=255)                   gACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCT‑CCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCg                                             >  2:724204/1‑140 (MQ=255)                      tGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGaa                                           >  2:1038885/1‑140 (MQ=255)                       gggCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGaaa                                          <  1:1038885/140‑1 (MQ=255)                        ggCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAAc                                         >  1:789408/1‑140 (MQ=255)                        ggCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAAc                                         >  2:495316/1‑140 (MQ=255)                        ggCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAAc                                         >  1:542146/1‑140 (MQ=255)                         gCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACt                                        <  1:615920/140‑1 (MQ=255)                           aGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTg                                      >  1:715528/1‑140 (MQ=255)                              aGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAg                                   <  2:789408/140‑1 (MQ=255)                                         cGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCGATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCg                        >  2:245677/1‑140 (MQ=255)                                         cGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCg                        >  1:546167/1‑140 (MQ=255)                                         cGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCg                        >  2:372118/1‑140 (MQ=255)                                         cGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCg                        >  1:870591/1‑140 (MQ=255)                                          gAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGt                       <  2:1248782/140‑1 (MQ=255)                                          gAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGt                       >  2:752590/1‑140 (MQ=255)                                          gAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGt                       >  2:769418/1‑140 (MQ=255)                                          gAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGt                       >  1:1206949/1‑140 (MQ=255)                                           aaaaaTCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTc                      <  2:1060733/140‑1 (MQ=255)                                                tCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTCAGGCg                 >  1:998829/1‑140 (MQ=255)                                                 cTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTCAGGCGt                <  2:886951/140‑1 (MQ=255)                                                 cTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTCAGGCGt                <  2:1141992/140‑1 (MQ=255)                                                 cTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTCAGGCGt                <  1:495316/140‑1 (MQ=255)                                                 cTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTCAGGCGt                <  1:40424/140‑1 (MQ=255)                                                   tCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCGATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTCAGGCGTGa              <  1:245677/140‑1 (MQ=14)                                                     caccaGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTCAGGCGTGATc            <  1:1089854/140‑1 (MQ=255)                                                      accaGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTCAGGCGTGATcc           <  2:546167/140‑1 (MQ=255)                                                               cTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATCCCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGACTTTTGACAAATCTGCTTACCCGAAACTTGCAGTTGCTTATCCGTCAGGCGTGATCCCTGATATGc  <  1:5748/140‑1 (MQ=2)                                                                                                                                        |                                                                   GGTTATGGTTTTGAGGAAGACTGGGCAGAAGTTTATACCTCGAAAAATCTTCCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGC  >  NC_000913/1431695‑1431896

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: ‑