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,8130CT97.7% 129.6 / ‑6.5 43*127* (TAG→TAAyigGPRK11371 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/21);  minor base G (0/1);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATATAT  >  NC_000913/4002695‑4002948
                                                                                                                      |                                                                                                                                       
tGGTGTTGCCCCCCAGAATCCAGGAAAACCATTTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAAt                                                                                                                         <  1:1256050/135‑1 (MQ=255)
tGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAAt                                                                                                                         >  2:1256050/1‑135 (MQ=255)
tGGTGTTGCCCCCAAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCCTTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                     >  2:29246/1‑139 (MQ=255)
        ccccccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCaa                                                                                                             >  1:208273/1‑139 (MQ=255)
                  tCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                     >  1:1395341/1‑121 (MQ=255)
                  tCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                     <  2:1395341/121‑1 (MQ=255)
                            ccATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAg                                                                                                            >  1:751653/1‑120 (MQ=255)
                            ccATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAg                                                                                                            <  2:751653/120‑1 (MQ=255)
                                                 aaTCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCaaaaaa                                                                    >  2:306208/1‑139 (MQ=255)
                                                        aTTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTAGTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATTCTGAAAGTGGAAAAATAACAACTTAATAAAATCATCg                                                             <  2:1160560/139‑1 (MQ=255)
                                                                       ttAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAAc                                              <  2:208273/139‑1 (MQ=255)
                                                                                 aTCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCa                                                                                                              >  2:1524001/1‑65 (MQ=255)
                                                                                 aTCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCa                                                                                                              <  1:1524001/65‑1 (MQ=255)
                                                                                   cGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCtt                                  <  1:306208/139‑1 (MQ=255)
                                                                                     ttCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCg                                >  2:907070/1‑139 (MQ=255)
                                                                                      tCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCgg                               <  1:907070/139‑1 (MQ=255)
                                                                                       cATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGaa                                                                                      <  2:1021112/83‑1 (MQ=255)
                                                                                       cATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGaa                                                                                      >  1:1021112/1‑83 (MQ=255)
                                                                                             ttCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt                        >  2:1382320/1‑139 (MQ=255)
                                                                                             ttCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt                        >  1:1865/1‑139 (MQ=255)
                                                                                                cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAATCAACGAACCATTGCCTGAACAGGGAAAATCTTTGGCTAACATtttg                    >  2:557704/1‑139 (MQ=255)
                                                                                                   aCTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATGACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGc                              >  1:1044431/1‑127 (MQ=255)
                                                                                                   aCTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGc                              <  2:1044431/127‑1 (MQ=255)
                                                                                                   aCTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGtgat                  >  2:303633/1‑139 (MQ=255)
                                                                                                    cTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGtgatg                 >  1:1154302/1‑139 (MQ=255)
                                                                                                        tAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATag             <  1:303633/139‑1 (MQ=255)
                                                                                                        tAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATag             >  1:244625/1‑139 (MQ=255)
                                                                                                         aaaTATTCCATAATAATTGTTTTAATAATTCAGAAGGGGAAACATTTAGAAACTTGAAGGTGGAAAAAAACAAATCCAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag           <  1:557704/139‑1 (MQ=255)
                                                                                                          aaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag           <  2:1154302/139‑1 (MQ=255)
                                                                                                          aaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag           <  1:1382320/139‑1 (MQ=255)
                                                                                                               ttcATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag           <  2:212968/132‑1 (MQ=255)
                                                                                                               tCCATAATTATTGTTTTGATGAATCAGTAGGTGCa                                                                                                              >  1:856786/1‑35 (MQ=255)
                                                                                                               tCCATAATTATTGTTTTGATGAATCAGTAGGTGCa                                                                                                              <  2:856786/35‑1 (MQ=255)
                                                                                                               tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag           >  2:265981/1‑134 (MQ=255)
                                                                                                               tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag           <  1:265981/134‑1 (MQ=255)
                                                                                                               tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATag             >  2:810571/1‑132 (MQ=255)
                                                                                                               tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATa              <  2:724091/131‑1 (MQ=255)
                                                                                                               tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTAAACAGGCAAAATCTTCGGCTATCATTGTGATGATa              >  1:724091/1‑131 (MQ=255)
                                                                                                               tCCATAATTATTGTTTTGATGAATCAGAAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATAATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag           >  1:212968/1‑134 (MQ=255)
                                                                                                               tCCATAAGTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATag             <  1:810571/132‑1 (MQ=255)
                                                                                                                  aTAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGatata   <  2:244625/139‑1 (MQ=255)
                                                                                                                  aTAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGatata   <  2:1865/139‑1 (MQ=255)
                                                                                                                   tAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGatatat  <  1:29246/139‑1 (MQ=255)
                                                                                                                      |                                                                                                                                       
TGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATATAT  >  NC_000913/4002695‑4002948

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