Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,313,070	T→C	T37A (ACC→GCC)	rbfA ←	30S ribosome binding factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,313,070	0	T	C	100.0%	132.8 / NA	38	T37A (ACC→GCC)	rbfA	30S ribosome binding factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/19);  total (19/19)

TGCGGATGAAACCAGAAGCTTCTTGCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGG  >  NC_000913/3312938‑3313202                                                                                                                                     |                                                                                                                                     tGCGGATGAAACCAGAAGCTTCTTGCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGtcatc                                                                                                                               >  2:1027425/1‑140 (MQ=255)     gATGAAACCAGAAGCTTCTTGCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGc                                                                                                                           <  2:1141803/140‑1 (MQ=255)      aTGAAACCAGAAGCTTCTTGCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGcc                                                                                                                          <  2:655971/140‑1 (MQ=255)         aaaCCAGAAGCTTCTTGCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAg                                                                                                                       >  2:51313/1‑140 (MQ=255)               gAAGCTTCTTGCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAgg                                                                                                                 <  1:291931/140‑1 (MQ=255)               gAAGCTTCTTGCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAgg                                                                                                                 <  1:51313/140‑1 (MQ=255)                      cttGCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTa                                                                                                          >  2:708299/1‑140 (MQ=255)                         gCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAAtt                                                                                                       <  2:7609/140‑1 (MQ=255)                             cGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCAc                                                                                                   <  2:586017/140‑1 (MQ=255)                                                aCCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGAt                                                                                <  1:756427/140‑1 (MQ=255)                                                 ccGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGAtc                                                                               <  1:607790/140‑1 (MQ=255)                                                    cgTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATctct                                                                            >  1:503909/1‑140 (MQ=255)                                                        ttCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCCGGCGAGGAGCTTTTATTTTACGCCGCAGGGTGAGAGCGGTCTCttttt                                                                        >  2:607790/1‑140 (MQ=255)                                                                tGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGACACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTcc                                                                >  1:1164906/1‑140 (MQ=255)                                                                 gTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCt                                                               <  2:319126/140‑1 (MQ=255)                                                                 gTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCt                                                               >  2:700792/1‑140 (MQ=255)                                                                                cacaTATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCg                                                <  1:1159771/140‑1 (MQ=255)                                                                                  caTATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCggg                                              >  1:166226/1‑140 (MQ=255)                                                                                       aCTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCggg                                              >  1:665651/1‑135 (MQ=255)                                                                                       aCTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCggg                                              <  2:665651/135‑1 (MQ=255)                                                                                          ttGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCaaa                                      >  2:294224/1‑140 (MQ=255)                                                                                                 aCGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTc                               >  1:610564/1‑140 (MQ=255)                                                                                                  cGCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCg                              >  2:200914/1‑140 (MQ=255)                                                                                                   gCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGc                             >  1:554627/1‑140 (MQ=255)                                                                                                   gCCAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGc                             <  1:16928/140‑1 (MQ=255)                                                                                                    ccAGGTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGcc                            <  2:503909/140‑1 (MQ=255)                                                                                                       ggTCGCGAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCata                         >  1:997967/1‑140 (MQ=255)                                                                                                            cgAGGCATTTCGACACCGGAAACGGCGGCCCTCATGGCCACGTGAAGATCTTTAAATTCACCGTGCTGGATTGATGGGGTTTCTTTTTACATTTACTGCGCTAGCCGCTCCGGGGGAGCCAATTATTTACCCCTAATaaa                    >  2:479023/1‑140 (MQ=255)                                                                                                            cgAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCttt                                <  1:723270/128‑1 (MQ=255)                                                                                                            cgAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCttt                                >  2:723270/1‑128 (MQ=255)                                                                                                            cgAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATaaa                    >  2:998693/1‑140 (MQ=255)                                                                                                            cgAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATaaa                    >  1:246032/1‑140 (MQ=255)                                                                                                            cgAGACATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCA‑TTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAAt                   >  1:368905/1‑140 (MQ=255)                                                                                                                 cATTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATtctc               <  1:700792/140‑1 (MQ=255)                                                                                                                  aTTTCGACACCGGAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTcc              >  1:687459/1‑140 (MQ=255)                                                                                                                    ttCGACACCGGAAACGGCGGTCATCATGCCAAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGGTCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTg            <  2:785514/140‑1 (MQ=255)                                                                                                                              gAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATtctc               >  2:1107047/1‑127 (MQ=255)                                                                                                                              gAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATtctc               <  1:1107047/127‑1 (MQ=255)                                                                                                                              gAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAgg  <  1:294224/140‑1 (MQ=255)                                                                                                                              gAAACGGCGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAgg  <  2:898537/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                     TGCGGATGAAACCAGAAGCTTCTTGCAACGCTTTGATGCCCGCTTTAACCGCGTCTTCATCTTTGTCGTTGAGGAACGTCACATATACTTTGGCATACGCCAGGTCGCGAGACATTTCGACACCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGG  >  NC_000913/3312938‑3313202

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