Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,335,010	T→C	S34S (AGT→AGC)	sfsB →	putative transcriptional regulator SfsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,335,010	0	T	C	100.0%	154.8 / NA	44	S34S (AGT→AGC)	sfsB	putative transcriptional regulator SfsB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/23);  total (21/23)

CTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGTTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGCTACCATGATCCGCAGACCCA  >  NC_000913/3334873‑3335132                                                                                                                                          |                                                                                                                           cTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTc                                                                                                                          >  2:77450/1‑140 (MQ=255)    aatgaaAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGctcctc                                                                                                                       >  2:187169/1‑140 (MQ=255)         aaGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGc                                                                                                                  >  1:358718/1‑140 (MQ=255)             gACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGc                                                                                                              <  1:498583/140‑1 (MQ=255)             gACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGc                                                                                                              <  2:581327/140‑1 (MQ=255)              aCTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCg                                                                                                             >  1:393638/1‑140 (MQ=255)              aCTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCg                                                                                                             >  1:342378/1‑140 (MQ=255)                            tGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATcgcg                                                                                               >  1:173293/1‑140 (MQ=255)                            tGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATcgcg                                                                                               >  2:129365/1‑140 (MQ=255)                            tGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATcgcg                                                                                               >  2:592675/1‑140 (MQ=255)                             gAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCgcgc                                                                                              >  2:554998/1‑140 (MQ=255)                              aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGcc                                                                                             <  1:353836/140‑1 (MQ=255)                              aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGcc                                                                                             >  1:95777/1‑140 (MQ=255)                              aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGcc                                                                                             >  2:438360/1‑140 (MQ=255)                               ggAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGccc                                                                                            <  1:567905/140‑1 (MQ=255)                                      tGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCg                                                                                     >  1:263286/1‑140 (MQ=255)                                       ggAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGa                                                                                    <  1:77450/140‑1 (MQ=255)                                               aTTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGgagag                                                                            >  2:606129/1‑140 (MQ=255)                                                 ttCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAgat                                                                          <  2:426780/140‑1 (MQ=255)                                                     ttGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTTGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGatt                                                                      >  2:673339/1‑140 (MQ=255)                                                       gACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGattat                                                                    >  2:675880/1‑140 (MQ=255)                                                              aTCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGaaa                                                             <  2:393638/140‑1 (MQ=255)                                                              aTCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGaaa                                                             <  1:321526/140‑1 (MQ=255)                                                                   gCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCt                                                        <  2:467559/140‑1 (MQ=255)                                                                       aCATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGa                                                    >  1:12653/1‑140 (MQ=255)                                                                         atcatTGCGGGTTTGCGCAAGAAGGGAACGTTAATGGCGGCTGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAAc                                                  <  1:652487/140‑1 (MQ=255)                                                                          tcatTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACt                                                 <  1:673339/140‑1 (MQ=255)                                                                          tcatTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACt                                                 <  1:545871/140‑1 (MQ=255)                                                                              tGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGAcc                                             <  1:462281/140‑1 (MQ=255)                                                                                        cgcAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTAt                                   >  1:399477/1‑140 (MQ=255)                                                                                               agGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCa                            <  1:187169/140‑1 (MQ=255)                                                                                                  gAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCa                         <  1:129365/140‑1 (MQ=255)                                                                                                  gAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCa                         <  2:263286/140‑1 (MQ=255)                                                                                                    aCGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACg                       <  2:358718/140‑1 (MQ=255)                                                                                                    aCGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACg                       <  2:173293/140‑1 (MQ=255)                                                                                                     gtttccaggcggcGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGc                      >  2:137722/8‑140 (MQ=255)                                                                                                     cGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGc                      <  1:58048/140‑1 (MQ=255)                                                                                                       tCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCa                            >  2:191629/1‑132 (MQ=255)                                                                                                       tCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCa                            <  1:191629/132‑1 (MQ=255)                                                                                                       tCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGCTa                    <  1:606129/140‑1 (MQ=255)                                                                                                       tCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGCTa                    >  1:248475/1‑140 (MQ=255)                                                                                                       tCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGCTa                    >  2:669066/1‑140 (MQ=255)                                                                                                                         cGCAGAAATGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGCTACCATGATCCGCAGACCCa  <  1:592675/140‑1 (MQ=255)                                                                                                                                 tGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGCTACCATGATCCGCa        >  2:216193/1‑126 (MQ=255)                                                                                                                                 tGGTTTGAGCTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGCTACCATGATCCGCa        <  1:216193/126‑1 (MQ=255)                                                                                                                                          |                                                                                                                           CTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGTTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGGCCATCACGCTACCATGATCCGCAGACCCA  >  NC_000913/3334873‑3335132

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