Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,185,807	T→C	F131S (TTT→TCT)	yqiG →	putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,185,807	0	T	C	100.0%	252.9 / NA	88	F131S (TTT→TCT)	yqiG	putative membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (56/32);  total (56/32)

GTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATTTCGATAAAGGGAAAATGCAATTATCCATATTTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTG  >  NC_000913/3185735‑3185878                                                                         |                                                                        gttgttACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGc                                                                      <  1:598075/76‑1 (MQ=255)    gttACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCAcc                                                                   >  2:3621598/1‑76 (MQ=255)    gttACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCAcc                                                                   <  1:3312430/76‑1 (MQ=255)     ttACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCt                                                                  <  2:5166657/76‑1 (MQ=255)      tACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTc                                                                 <  2:262986/76‑1 (MQ=255)       aCCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCa                                                                >  2:2733906/1‑76 (MQ=255)       aCCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCa                                                                >  1:3954777/1‑76 (MQ=255)        cccTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAg                                                               <  1:4859933/76‑1 (MQ=255)         ccTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAgg                                                              >  2:4380912/1‑76 (MQ=255)           tATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCa                                                            >  2:4560965/1‑76 (MQ=255)           tATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCa                                                            <  2:5486635/76‑1 (MQ=255)            aTCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCAt                                                           >  2:4327443/1‑76 (MQ=255)            aTCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCAt                                                           <  2:1917465/76‑1 (MQ=255)            aTCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCAt                                                           >  2:4284159/1‑76 (MQ=255)              cGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATgg                                                         <  2:5604140/76‑1 (MQ=255)               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCGGGCATGGt                                                        >  2:5696082/1‑76 (MQ=255)               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCGGGCATGGt                                                        <  1:5696084/76‑1 (MQ=255)               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGt                                                        >  2:4220517/1‑76 (MQ=255)               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGt                                                        >  1:769474/1‑76 (MQ=255)               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGt                                                        <  1:5827783/76‑1 (MQ=255)               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGt                                                        >  1:5225605/1‑76 (MQ=255)               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGt                                                        >  1:4415167/1‑76 (MQ=255)                aaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGtt                                                       <  2:1861199/76‑1 (MQ=255)                aaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGtt                                                       <  1:2094247/76‑1 (MQ=255)                aaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGtt                                                       <  1:4889714/76‑1 (MQ=255)                aaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGtt                                                       >  1:1392305/1‑76 (MQ=255)                 aaaagaaacaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTa                                                      <  1:5371815/76‑1 (MQ=255)                 aaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTa                                                      >  1:4418359/1‑76 (MQ=255)                 aaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTa                                                      <  2:3039852/76‑1 (MQ=255)                 aaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTa                                                      <  2:5599252/76‑1 (MQ=255)                  aaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaa                                                     >  2:1983810/1‑76 (MQ=255)                  aaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaa                                                     <  2:4267542/76‑1 (MQ=255)                  aaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCAGCACCTCAGGCATGGTTaa                                                     >  2:4393372/1‑76 (MQ=255)                     gaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAAt                                                  >  2:3848529/1‑76 (MQ=255)                     gaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAAt                                                  <  2:3206337/76‑1 (MQ=255)                      aaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATa                                                 >  1:2219229/1‑76 (MQ=255)                         aaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACaa                                              <  2:723817/76‑1 (MQ=255)                           cAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAg                                            >  2:5467682/1‑76 (MQ=255)                             aaTTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGAt                                          >  1:4607348/1‑76 (MQ=255)                              aTTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATg                                         <  2:4266097/76‑1 (MQ=255)                              aTTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATg                                         <  2:3067563/76‑1 (MQ=255)                              aTTACAACTTATCTCGATAAAAGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATg                                         >  2:3475714/1‑76 (MQ=255)                               ttACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGc                                        >  1:3678830/1‑76 (MQ=255)                               ttACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGc                                        <  1:5227694/76‑1 (MQ=255)                                tACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCa                                       >  1:1333531/1‑76 (MQ=255)                                tACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCa                                       >  1:5771669/1‑76 (MQ=255)                                 aCAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCaa                                      >  1:1447897/1‑76 (MQ=255)                                 aCAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCaa                                      >  2:5105577/1‑76 (MQ=255)                                  cAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCaaa                                     <  2:2219229/76‑1 (MQ=255)                                   aaCTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAAc                                    >  2:2855332/1‑76 (MQ=255)                                   aaCTGATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAAc                                    <  2:223269/76‑1 (MQ=255)                                    aCTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACt                                   <  1:5541115/76‑1 (MQ=255)                                    aCTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACt                                   >  1:3450244/1‑76 (MQ=255)                                     cTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTg                                  >  2:817237/1‑76 (MQ=255)                                      ttATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTgg                                 >  2:4425265/1‑76 (MQ=255)                                      ttATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTgg                                 >  1:5828333/1‑76 (MQ=255)                                        atctCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGAc                               >  2:5464828/4‑76 (MQ=255)                                        atctCGAGAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGAc                               >  1:2410054/4‑76 (MQ=255)                                          ctCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGc                             >  2:95372/2‑76 (MQ=255)                                           tCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcc                            >  2:3257410/1‑76 (MQ=255)                                           tCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcc                            >  1:1462064/1‑76 (MQ=255)                                            cGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcct                           >  2:3042968/1‑76 (MQ=255)                                            cGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcct                           >  1:3008570/1‑76 (MQ=255)                                            cGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcct                           <  2:4418358/76‑1 (MQ=255)                                            cGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGccc                           >  2:4879235/1‑75 (MQ=255)                                             gATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGATGcctc                          >  1:5799964/1‑76 (MQ=255)                                              aTAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcctcc                         >  1:1327303/1‑76 (MQ=255)                                               tAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcctcct                        >  2:3539161/1‑76 (MQ=255)                                               tAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcctcct                        >  2:3504739/1‑76 (MQ=255)                                                 aaGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGa                      <  2:5141834/76‑1 (MQ=255)                                                  aGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTgca                     >  1:3864636/1‑74 (MQ=255)                                                  aGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGaa                     >  2:5087939/1‑76 (MQ=255)                                                  aGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGaa                     <  2:592301/76‑1 (MQ=255)                                                    ggAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACt                   >  1:4723036/1‑76 (MQ=255)                                                     gAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACtt                  >  2:5359923/1‑76 (MQ=255)                                                     gAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACtt                  >  1:2318260/1‑76 (MQ=255)                                                         aTGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGa              >  1:4321281/1‑76 (MQ=255)                                                         aTGAAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACGTTGGa              >  1:1200870/1‑76 (MQ=255)                                                             aaTTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCa          <  2:3415599/76‑1 (MQ=255)                                                              aTTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCAt         >  1:1261182/1‑76 (MQ=255)                                                                tATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAGCTGGTCGCCTCCTGAACTTTGGAATCATgg       <  2:3224257/76‑1 (MQ=255)                                                                tATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATgg       <  1:695496/76‑1 (MQ=255)                                                                tATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATgg       >  2:3396704/1‑76 (MQ=255)                                                                 aTCCATATCTGCACCTCAGGCATGGGTAAAATACAAAGATGCAAACTGGACCACTCCTCAACTGTGGAATCATGGt      >  1:975032/1‑76 (MQ=255)                                                                  tCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTa     >  2:5354396/1‑76 (MQ=255)                                                                    cATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTAtt   >  1:1052398/1‑76 (MQ=255)                                                                    cATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTAtt   <  2:3954776/76‑1 (MQ=255)                                                                     atatCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTg  >  1:1267265/1‑76 (MQ=255)                                                                         |                                                                        GTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATTTCGATAAAGGGAAAATGCAATTATCCATATTTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTG  >  NC_000913/3185735‑3185878

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