Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,129,332 A→G V498V (GTA→GTG)  murJ → putative lipid II flippase MurJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,129,3320AG96.8% 101.8 / ‑4.4 31V498V (GTA→GTG) murJputative lipid II flippase MurJ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (1/29);  total (1/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TGGTGATGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCC  >  NC_000913/1129196‑1129468
                                                                                                                                        |                                                                                                                                        
tGGTGATGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTg                                                                                                                                       <  1:107640/140‑1 (MQ=255)
      tGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTc                                                                                                                                 <  1:255917/140‑1 (MQ=255)
             cGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTa                                                                                                                          <  2:684321/140‑1 (MQ=255)
                             ttACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCccgccg                                                                                                          <  2:805854/140‑1 (MQ=255)
                               aCATATCGTGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGa                                                                                                        <  1:449124/140‑1 (MQ=255)
                               aCATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGa                                                                                                        <  2:110899/140‑1 (MQ=255)
                               aCATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGa                                                                                                        <  1:243606/140‑1 (MQ=255)
                               aCATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGa                                                                                                        <  2:753103/140‑1 (MQ=255)
                               aCATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGa                                                                                                        <  1:566045/140‑1 (MQ=255)
                               aCATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGa                                                                                                        <  2:796722/140‑1 (MQ=255)
                                                             cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtg                                                                          <  1:823380/140‑1 (MQ=255)
                                                               tGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgca                                                                        <  1:205537/140‑1 (MQ=255)
                                                                     ggCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCgg                                                                  <  2:100899/140‑1 (MQ=255)
                                                                                    tAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATccggcctg                                                   <  1:391777/140‑8 (MQ=255)
                                                                                        ggCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGaa                                               <  2:1258066/140‑1 (MQ=255)
                                                                                        ggCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGaa                                               <  1:639553/140‑1 (MQ=255)
                                                                                          cGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAAt                                             <  1:187732/140‑1 (MQ=255)
                                                                                                    ggCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTAccg                                   <  1:871484/140‑1 (MQ=255)
                                                                                                         ggATTGCCGCGTACTTCGCTGCCCTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCcgcg                              <  1:1145865/140‑1 (MQ=255)
                                                                                                              gCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                         <  1:607207/140‑1 (MQ=255)
                                                                                                              gCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                         <  1:482389/140‑1 (MQ=255)
                                                                                                              gCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                         <  1:789169/140‑1 (MQ=255)
                                                                                                              gCCGCGTACTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                         <  2:1251810/140‑1 (MQ=255)
                                                                                                              accGCGTCCTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGGCGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                         <  1:1285141/139‑1 (MQ=255)
                                                                                                                      cTTCGCTGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAg                 <  1:1470709/140‑1 (MQ=255)
                                                                                                                            tGCACTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGa           <  1:973318/140‑1 (MQ=255)
                                                                                                                                cTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt       <  1:1045822/140‑1 (MQ=255)
                                                                                                                                cTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt       <  1:615906/140‑1 (MQ=255)
                                                                                                                                cTGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt       <  1:307727/140‑1 (MQ=255)
                                                                                                                                 tGGCGGTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGGCGGGGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGGTGTTGAGGTCTGACCGtt      >  1:806143/1‑140 (MQ=255)
                                                                                                                                     ggTGCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGcc  <  1:69594/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TGGTGATGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCC  >  NC_000913/1129196‑1129468

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