Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,002,813 C→T *127* (TAG→TAA)  yigG ← PRK11371 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,8130CT100.0% 134.1 / NA 42*127* (TAG→TAA) yigGPRK11371 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/22);  total (20/22)

AAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAG  >  NC_000913/4002677‑4002939
                                                                                                                                        |                                                                                                                              
aaaTGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTa                                                                                                                              >  1:517704/1‑139 (MQ=255)
                  tGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGTTATTTCCATTTTTCACTTAATGAACCATCCTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                            >  2:612219/1‑139 (MQ=255)
                  tGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                            >  1:336365/1‑139 (MQ=255)
                       ttGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                       <  2:517704/139‑1 (MQ=255)
                         gCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCa                                                                                                     <  2:4518/139‑1 (MQ=255)
                             cccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGGACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCa                                                                                                 <  2:454537/139‑1 (MQ=255)
                                  aaTCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGc                                                                                            <  1:715359/139‑1 (MQ=255)
                                                      aTCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATa                                                                        >  1:881596/1‑139 (MQ=255)
                                                                  aaaTCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTtgatga                                                                                                                  <  2:724553/85‑1 (MQ=255)
                                                                  aaaTCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTtgatga                                                                                                                  >  1:724553/1‑85 (MQ=255)
                                                                         tATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCa                                                                                                 >  2:649707/1‑95 (MQ=255)
                                                                         tATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCa                                                                                                 <  1:649707/95‑1 (MQ=255)
                                                                         tATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAatcatc                                                     <  2:336365/139‑1 (MQ=255)
                                                                                aTTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCAt                                              >  1:239163/1‑139 (MQ=255)
                                                                                                          atatatTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAg                                                                                                   <  2:432231/60‑1 (MQ=255)
                                                                                                          atatatTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAg                                                                                                   >  1:432231/1‑60 (MQ=255)
                                                                                                           tatatTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCa                                                                                                     <  1:416096/57‑1 (MQ=255)
                                                                                                           tatatTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCa                                                                                                     >  2:416096/1‑57 (MQ=255)
                                                                                                              atTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCa                >  1:884797/1‑139 (MQ=255)
                                                                                                               ttCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt               <  2:239163/139‑1 (MQ=255)
                                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                       >  1:905234/1‑48 (MQ=255)
                                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                       <  2:905234/48‑1 (MQ=255)
                                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgt            >  1:363834/1‑139 (MQ=255)
                                                                                                                    tACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt               >  2:325250/1‑134 (MQ=255)
                                                                                                                    tACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt               <  1:325250/134‑1 (MQ=255)
                                                                                                                     aCTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                            <  1:898130/40‑1 (MQ=255)
                                                                                                                     aCTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                            >  2:898130/1‑40 (MQ=255)
                                                                                                                     aCTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGtgat         >  2:749552/1‑139 (MQ=255)
                                                                                                                        cATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTgatgat      <  1:216129/139‑1 (MQ=255)
                                                                                                                         aTAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATa     >  2:222469/1‑139 (MQ=255)
                                                                                                                         aTAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATa     <  1:749552/139‑1 (MQ=255)
                                                                                                                          tAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATag    <  2:881596/139‑1 (MQ=255)
                                                                                                                          tAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGACTATCATTGTGATGATag    <  1:612219/139‑1 (MQ=255)
                                                                                                                            aaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag  <  2:884797/139‑1 (MQ=255)
                                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag  <  1:883103/134‑1 (MQ=255)
                                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag  >  2:883103/1‑134 (MQ=255)
                                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCt                          >  2:415781/1‑108 (MQ=255)
                                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCt                          <  1:415781/108‑1 (MQ=255)
                                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTgatgat      >  2:61919/1‑128 (MQ=255)
                                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTgatgat      <  1:61919/128‑1 (MQ=255)
                                                                                                                                     tAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtg             >  2:189413/1‑119 (MQ=255)
                                                                                                                                     tAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtg             <  1:189413/119‑1 (MQ=255)
                                                                                                                                        |                                                                                                                              
AAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAG  >  NC_000913/4002677‑4002939

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