Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,313,236	T→C	intergenic (‑58/+106)	rbfA ← / ← infB	30s ribosome binding factor/translation initiation factor IF‑2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,313,236	0	T	C	100.0%	147.9 / NA	40	intergenic (‑58/+106)	rbfA/infB	30s ribosome binding factor/translation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/15);  total (25/15)

GAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACGTTGGATCTCGATGATTTCGAATAC  >  NC_000913/3313088‑3313380                                                                                                                                                     |                                                                                                                                                 gAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGTATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCgg                                                                                                                                                >  2:525192/1‑151 (MQ=255)   ggATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGa                                                                                                                                               >  2:282599/1‑150 (MQ=255)         ttAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAggg                                                                                                                                           <  1:319273/148‑1 (MQ=255)         ttAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAgg                                                                                                                                            <  2:170581/147‑1 (MQ=255)                            tGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATa                                                                                                                       >  2:511923/1‑149 (MQ=255)                               gagCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGAcg                                                                                                                  >  2:194940/1‑151 (MQ=255)                                        ctTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACgcgc                                                                                                               <  2:288304/145‑1 (MQ=255)                                         ttttcGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTc                                                                                                        >  1:35580/1‑151 (MQ=255)                                                aTTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATc                                                                                                   >  1:356974/1‑149 (MQ=255)                                                 tttCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAgg                                                                                                >  2:607783/1‑151 (MQ=255)                                                       gcgcTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATcc                                                                                          >  2:382353/1‑151 (MQ=255)                                                       gcgcTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATcc                                                                                          >  2:198006/1‑151 (MQ=255)                                                                     gggCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGc                                                                            >  2:287930/1‑151 (MQ=255)                                                                      ggCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCg                                                                           >  2:556010/1‑151 (MQ=255)                                                                      ggCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCg                                                                           >  1:566365/1‑151 (MQ=255)                                                                        cGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGt                                                                          >  2:586200/1‑150 (MQ=255)                                                                         gACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTaa                                                                        <  1:197529/151‑1 (MQ=255)                                                                           ccAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAAc                                                                      >  2:558406/1‑151 (MQ=255)                                                                              aaTTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGT‑AACGCCt                                                                  >  1:308282/1‑151 (MQ=255)                                                                                 tCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCtt                                                                 >  1:439479/1‑150 (MQ=255)                                                                                         cATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCAt                                                        <  2:137393/151‑1 (MQ=255)                                                                                         cATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCAt                                                        >  2:249020/1‑151 (MQ=255)                                                                                            aataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATgg                                                      <  2:648574/150‑1 (MQ=255)                                                                                                  ttCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCt                                                >  2:6717/1‑150 (MQ=255)                                                                                                      ccTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTaaaa                                            <  1:170357/150‑1 (MQ=255)                                                                                                            aaaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTa                                      >  2:538224/1‑150 (MQ=255)                                                                                                             aaaaaGGGGCTGTTAGCCCCTTTTTTAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg                                    <  1:410846/151‑1 (MQ=255)                                                                                                             aaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg                                    <  2:298381/151‑1 (MQ=255)                                                                                                              aaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGc                                   >  1:254533/1‑151 (MQ=255)                                                                                                               aaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCa                                  <  2:28521/151‑1 (MQ=255)                                                                                                               aaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCa                                  <  2:231238/151‑1 (MQ=255)                                                                                                                           aGCCCCTTTTTAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACGTTgg                      >  1:29593/1‑151 (MQ=255)                                                                                                                                  ttttAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACGTTGGAtctc                 <  2:13508/149‑1 (MQ=255)                                                                                                                                    ttAAAATTAATTTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACGTTGGATCTCgatg             >  2:249194/1‑151 (MQ=255)                                                                                                                                        aattaattTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACGTTGGATCTCGATGAttt         <  1:519574/151‑1 (MQ=255)                                                                                                                                         attaattTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGGACGTTGGATCTCGATGATTTc        >  2:321259/1‑151 (MQ=255)                                                                                                                                           taattTCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACGTTGGATCTCGATGATTTCg       >  1:207427/1‑150 (MQ=255)                                                                                                                                                tCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAACCGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACGTTGGATCTCGATGATTTCGAATAc  <  1:140165/150‑1 (MQ=255)                                                                                                                                                tCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCt                                        <  2:112973/112‑1 (MQ=255)                                                                                                                                                tCAGGCGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCt                                        >  1:112970/1‑112 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACGTTGGATCTCGATGATTTCGAATAC  >  NC_000913/3313088‑3313380

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