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%	218.6 / NA	79	F131S (TTT→TCT)	yqiG	putative membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (57/22);  total (57/22)

GTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATTTCGATAAAGGGAAAATGCAATTATCCATATTTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTAT  >  NC_000913/3185735‑3185876                                                                         |                                                                      gttgttACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGc                                                                    >  2:1726835/1‑76 (MQ=255)  ttgttACCCTATCGAAAAAGAAAAACAAATTCCAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGcc                                                                   >  2:3367502/1‑75 (MQ=255)    gttACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCAcc                                                                 <  1:3011657/76‑1 (MQ=255)        cccTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAg                                                             >  1:137154/1‑76 (MQ=255)        cccTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAg                                                             >  1:2407976/1‑76 (MQ=255)          cTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGc                                                           <  2:1156920/76‑1 (MQ=255)           tATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCa                                                          >  2:925290/1‑76 (MQ=255)           tATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCa                                                          >  1:4300040/1‑76 (MQ=255)            aTCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCAt                                                         >  2:3849424/1‑76 (MQ=255)              cGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATgg                                                       >  2:1453917/1‑76 (MQ=255)               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGt                                                      <  2:1502910/76‑1 (MQ=255)               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGt                                                      >  2:799505/1‑76 (MQ=255)                aaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGtt                                                     >  2:4580433/1‑76 (MQ=255)                 aaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTa                                                    >  2:2532199/1‑76 (MQ=255)                  aaagaaaaaCAAATTCCAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCCCCTCAGGCATGTTTaa                                                   >  2:335612/1‑76 (MQ=255)                   aagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaaa                                                  >  1:2841158/1‑76 (MQ=255)                   aagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaaa                                                  >  2:2899786/1‑76 (MQ=255)                   aagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaaa                                                  >  1:4682729/1‑76 (MQ=255)                   aagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaaa                                                  >  2:3061797/1‑76 (MQ=255)                   aagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCCGCACCTCAGGCATGGTTaaa                                                  <  2:743358/76‑1 (MQ=255)                   aagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAAATATCCATATCTGCACCTCAGGCATGGTTaaa                                                  <  2:4209382/76‑1 (MQ=255)                    agaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaaaa                                                 <  2:1091643/76‑1 (MQ=255)                    agaaaaaCAAATTACAACGTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaaaa                                                 <  2:1120506/76‑1 (MQ=255)                     gaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAAt                                                >  1:1461207/1‑76 (MQ=255)                     gaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAAt                                                >  2:4516070/1‑76 (MQ=255)                       aaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATAc                                              >  1:2517592/1‑76 (MQ=255)                        aaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACa                                             >  2:4742981/1‑76 (MQ=255)                        aaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACa                                             >  1:1715744/1‑76 (MQ=255)                           cAAATTACAACTTATCTCGATGAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAg                                          <  1:84635/76‑1 (MQ=255)                           cAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAg                                          >  2:4387139/1‑76 (MQ=255)                            aaaTTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGa                                         >  2:4355727/1‑76 (MQ=255)                             aaTTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGAt                                        <  1:2629773/76‑1 (MQ=255)                             aaTTACAACTTACCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGAt                                        >  2:2349212/1‑76 (MQ=255)                              aTTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATg                                       <  2:507892/76‑1 (MQ=255)                              aTTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATg                                       <  2:2585127/76‑1 (MQ=255)                              aTTACAACTTATCGCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATg                                       <  2:2608122/76‑1 (MQ=255)                               ttACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCCGCCATGGTTAAAATAAAAAGATga                                      >  2:43343/1‑75 (MQ=255)                               ttACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGc                                      >  1:1285029/1‑76 (MQ=255)                                tACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCa                                     >  1:2638878/1‑76 (MQ=255)                                tACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGCAATGGTTAAAATACAAAGATGCa                                     >  2:993245/1‑76 (MQ=255)                                 aCAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCaa                                    >  1:2764896/1‑76 (MQ=255)                                 aCAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCaa                                    >  1:4286502/1‑76 (MQ=255)                                  cAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCaaa                                   >  1:4713145/1‑76 (MQ=255)                                   aaCTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAAc                                  >  2:4387952/1‑76 (MQ=255)                                   aaCTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCCCAGGCATGGTTAAAATACAAAGATGCAAAc                                  >  2:4740883/1‑76 (MQ=255)                                      ttATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTgg                               >  2:348202/1‑76 (MQ=255)                                         tctCGATAAAGGGAAAATGCAATTATACATAACTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACg                            >  1:3524294/3‑76 (MQ=255)                                            cGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcct                         >  2:2370889/1‑76 (MQ=255)                                            cGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcct                         >  2:3040766/1‑76 (MQ=255)                                              aTAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcctcc                       >  2:2386839/1‑76 (MQ=255)                                              aTAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcctcc                       >  1:4355877/1‑76 (MQ=255)                                              aTAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcctcc                       <  2:2764896/76‑1 (MQ=255)                                               tAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGcctcct                      <  2:4470570/76‑1 (MQ=255)                                                aaaGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTg                     <  1:920662/76‑1 (MQ=255)                                                 aaGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGa                    >  2:1144265/1‑76 (MQ=255)                                                   gggAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAAc                  >  2:3384489/1‑76 (MQ=255)                                                   gggAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAAc                  >  1:1465125/1‑76 (MQ=255)                                                   gggAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAAc                  >  1:3111131/1‑76 (MQ=255)                                                    ggAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACt                 >  1:301260/1‑76 (MQ=255)                                                    ggAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACt                 <  2:1384595/76‑1 (MQ=255)                                                     gAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACtt                >  1:2042342/1‑76 (MQ=255)                                                      aaaaTGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACttt               >  1:856215/1‑76 (MQ=255)                                                        aaTGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTgg             >  1:4534282/1‑76 (MQ=255)                                                          tGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGaa           >  2:50428/1‑76 (MQ=255)                                                          tGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGaa           >  1:4562306/1‑76 (MQ=255)                                                            cAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATc         >  2:1945947/1‑76 (MQ=255)                                                            cAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATc         >  1:2401854/1‑76 (MQ=255)                                                             aaTTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCa        <  2:4682729/76‑1 (MQ=255)                                                             aaTTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCa        >  2:3487491/1‑76 (MQ=255)                                                             aaTTAGCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTCGACGCCTAAAGAAAATTGGAATCa        >  2:4612648/1‑76 (MQ=255)                                                              aTTATCCCTATCTGCACCTCAGGCATGGTTAAAATCCAAAGATGCACACTGGACGCTTCCTGAACTTTGGAATCAt       <  2:4713486/76‑1 (MQ=255)                                                               ttATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATg      <  1:4516070/76‑1 (MQ=255)                                                                tATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCAttg     >  2:272325/1‑74 (MQ=255)                                                                tATCCATAACTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATgg     <  2:3524294/76‑1 (MQ=255)                                                                 aTCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGt    >  2:4462788/1‑76 (MQ=255)                                                                 aTCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTCTGGAATCATGGt    >  2:903080/1‑76 (MQ=255)                                                                   ccATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTAt  >  1:1402307/1‑76 (MQ=255)                                                                   ccATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTAt  <  2:4812352/76‑1 (MQ=255)                                                                   ccATATCGGCACCGCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTAt  <  2:4822615/76‑1 (MQ=255)                                                                         |                                                                      GTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATTTCGATAAAGGGAAAATGCAATTATCCATATTTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTAT  >  NC_000913/3185735‑3185876

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