Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,185,776	T→C	F121L (TTC→CTC)	yqiG →	putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,185,776	0	T	C	100.0%	220.8 / NA	80	F121L (TTC→CTC)	yqiG	putative membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (52/28);  total (52/28)

TATTGAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATTTCGATAAAGGGAAAATGCAATTATCCATATTTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGG  >  NC_000913/3185703‑3185847                                                                          |                                                                        tATTGAAAAAACTCAGCAAACCCTGATCGATGGTTGTTACCCTATCGAAACAGAAAAACCAATTACAACTTAtctc                                                                       >  1:96118/1‑76 (MQ=255)  aTTGAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCg                                                                      <  1:2219226/76‑1 (MQ=255)    tGAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGAt                                                                    <  1:2679572/76‑1 (MQ=255)      aaaaaaCTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATaa                                                                  >  1:3707023/1‑76 (MQ=255)         aaaCTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGACAAACAAATTACAACTTATCTCGATAAAgg                                                               >  2:1648546/1‑76 (MQ=255)         aaaCTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAgg                                                               >  2:1607862/1‑76 (MQ=255)           aCTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGa                                                             >  2:1423237/1‑76 (MQ=255)             tCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGaaa                                                           <  1:4850715/76‑1 (MQ=255)             tCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGaaa                                                           >  1:3331801/1‑76 (MQ=255)              cAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGaaaa                                                          >  1:1205460/1‑76 (MQ=255)              cAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGaaaa                                                          <  2:3290317/76‑1 (MQ=255)              cAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGaaaa                                                          >  1:1541027/1‑76 (MQ=255)               aGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAAt                                                         <  1:3319103/76‑1 (MQ=255)                  aaaaCCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCa                                                      >  2:1123288/1‑76 (MQ=255)                   aaaCCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAACCAAATTACAACTTATCTCGATAAAGGGAAAATGCaa                                                     >  1:787369/1‑76 (MQ=255)                   aaaCCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCaa                                                     >  1:2580087/1‑76 (MQ=255)                   aaaCCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCaa                                                     <  2:4024151/76‑1 (MQ=255)                   aaaCCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCaa                                                     >  2:377407/1‑76 (MQ=255)                    aaCCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAAt                                                    >  1:2623843/1‑76 (MQ=255)                        tGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATc                                                <  2:1764859/76‑1 (MQ=255)                        tGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATc                                                <  2:2273203/76‑1 (MQ=255)                         gATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATcc                                               <  1:2809211/76‑1 (MQ=255)                         gATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATcc                                               >  2:1182133/1‑76 (MQ=255)                         gATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATcc                                               <  1:99874/76‑1 (MQ=255)                         gATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATcc                                               <  1:4678276/76‑1 (MQ=255)                         gATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATcc                                               >  2:4222690/1‑76 (MQ=255)                         gATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGAAATTATcc                                               >  2:87969/1‑76 (MQ=255)                         gATCGATGGTTGGTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATcc                                               <  1:2854083/76‑1 (MQ=255)                          aTCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATccc                                              <  2:2356699/76‑2 (MQ=255)                          aTCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCa                                              >  1:3161070/1‑76 (MQ=255)                             gATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCatat                                           >  1:2419237/1‑76 (MQ=255)                              aTGGTTGTTACCATACCGAAAAAGAAAAACAAATTACAAATTATCGCGATAAAGGGTAAATGCAATTATCCatatc                                          >  2:983654/1‑75 (MQ=255)                               tGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCatatct                                         <  1:2091240/76‑3 (MQ=255)                               tGGGTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCatatct                                         <  2:4096233/76‑3 (MQ=255)                                ggTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCatatctg                                        <  2:2180301/76‑4 (MQ=255)                                 gttgttACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGc                                       >  2:1726835/1‑76 (MQ=255)                                  ttgttACCCTATCGAAAAAGAAAAACAAATTCCAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGcc                                      >  2:3367502/1‑75 (MQ=255)                                    gttACCCTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCAcc                                    <  1:3011657/76‑1 (MQ=255)                                        cccTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAg                                >  1:2407976/1‑76 (MQ=255)                                        cccTATCGAAAAAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAg                                >  1:137154/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:799505/1‑76 (MQ=255)                                               gaaaaagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGt                         <  2:1502910/76‑1 (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:4682729/1‑76 (MQ=255)                                                   aagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaaa                     >  1:2841158/1‑76 (MQ=255)                                                   aagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTaaa                     >  2:2899786/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                   >  2:4516070/1‑76 (MQ=255)                                                     gaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAAt                   >  1:1461207/1‑76 (MQ=255)                                                       aaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATAc                 >  1:2517592/1‑76 (MQ=255)                                                        aaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACa                >  1:1715744/1‑76 (MQ=255)                                                        aaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACa                >  2:4742981/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:2585127/76‑1 (MQ=255)                                                              aTTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATg          <  2:507892/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:4286502/1‑76 (MQ=255)                                                                 aCAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCaa       >  1:2764896/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)                                                                          |                                                                        TATTGAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATTTCGATAAAGGGAAAATGCAATTATCCATATTTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGG  >  NC_000913/3185703‑3185847

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