Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,636,042 A→G 12.0% L438S (TTG→TCG)  der ← GTPase; multicopy suppressor of ftsJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,636,0420AG12.0% 80.7 / 2.4 33L438S (TTG→TCG) derGTPase; multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base A (17/12);  new base G (2/2);  total (19/14)
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.28e-01

GTGCTTCATCAGACGCTTACGTTTACGCATCTGGGTTGGCGTCAGGGTGTTACGCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATC  >  NC_000913/2635900‑2636188
                                                                                                                                              |                                                                                                                                                  
gtgCTTCATCAGACGCTTACGTTTACGCATCTGGGTTGGCGTCAGGGTGTTACGCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCGAGTAgc                                                                                                                                              >  2:347533/1‑149 (MQ=255)
gtgCTTCATCAGACGCTTACGTTTACGCAAACGGACTATCGTCAGGGTGTTACGCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCGAGTAgc                                                                                                                                              >  2:347450/1‑149 (MQ=255)
   cTTCATCAGACGCTTACGTTTACGCATCTGGGTTGGCGTCAGGGTGTTACGCTTATTCGCATACGGGTTTTACCCTTCTTTGAACTGAATAAGAATCGGAGATACAATTACGTCCAGCGATTTGCGGAAGTAGTTGATCAGGAAGCGCt                                                                                                                                           >  2:602888/1‑149 (MQ=255)
         caGACGCTTACGTTTACGCATCTGGGTTGGCGTCAGGGTGATACGCTTATACGCATACGGGCTCTCCCCTCCTTTGAACTGAATACGAACCGGCGACCCCATTACGTCCCGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTCGTAgc                                                                                                                                     >  2:26712/1‑148 (MQ=255)
           gtcGCTTACGTTTACGCATCTGGGTTGGCGTCAGGGTGTTACGCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGaa                                                                                                                                   <  2:114001/147‑1 (MQ=255)
                         cGCATCTGGGTTGGCGTCAGGGTGTTACGCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCttt                                                                                                                     <  2:520747/149‑1 (MQ=255)
                          gCATCTGGGTTGGCGTCAGGGTGATACGCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTc                                                                                                                    >  1:493505/1‑149 (MQ=255)
                              cggggTTGCCAACAGGCTATCACGCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCt                                                                                                                <  2:303790/147‑1 (MQ=255)
                               tGGGTTGGCGACATGTTGTGACGCTGATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCTATTTGCTGAAGTAGTTCTTCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTg                                                                                                               <  1:445166/149‑1 (MQ=255)
                                                   aCGCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAAGCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACaa                                                                                           >  1:741986/1‑149 (MQ=255)
                                                     gCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATc                                                                                         <  1:183036/149‑1 (MQ=255)
                                                         aTTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACGACAATagggc                                                                                     >  2:371253/1‑147 (MQ=255)
                                                          ttCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATcggcgg                                                                                    <  2:12613/149‑1 (MQ=255)
                                                                    ggTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAacca                                                                          >  1:142710/1‑149 (MQ=255)
                                                                     gTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGGAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGCCTTTCACCCGATTACCGTGAATCACCACAATCGCCGGGTTATAaccgc                                                                         >  2:312484/1‑147 (MQ=255)
                                                                            ccTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATCACCACCGGCGTg                                                                  >  2:245419/1‑149 (MQ=255)
                                                                            ccTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGta                                                                  >  2:431117/1‑148 (MQ=255)
                                                                             cttcttTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTGATAACCACCGGCGTgg                                                                 >  2:666618/1‑149 (MQ=255)
                                                                                     aaCTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGCGGGCATAttt                                                         >  2:279326/1‑149 (MQ=255)
                                                                                       cTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCGAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCa                                                       <  1:347450/149‑1 (MQ=255)
                                                                                       cTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCGAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCa                                                       <  1:347533/149‑1 (MQ=255)
                                                                                         gAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGc                                                     <  1:305626/149‑1 (MQ=255)
                                                                                         gAATACGAATCGGCCATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGCAATGAGGCAGGTCTTTCACCTGATTACAGTGAATCACCACGCAAGGCGGGTTATAACCACCGACCTCGGTATATTTCAGc                                                     >  2:255367/1‑149 (MQ=255)
                                                                                          aaTACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGCTATAACCACCGGCGTGGGCATATTTCAGca                                                    >  2:397305/1‑148 (MQ=255)
                                                                                                  tCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACacga                                            <  1:612851/149‑1 (MQ=255)
                                                                                                      cGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACgacgac                                        <  2:501909/149‑1 (MQ=255)
                                                                                                          cccATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACGACGACcgcg                                    <  1:392511/149‑1 (MQ=255)
                                                                                                                   tCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACGACGACCGCGTACCAgcgg                           >  1:611970/1‑149 (MQ=255)
                                                                                                                     cAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCACATTTCAGCTTCACACCACGACCGCGAACCCgcggcg                         >  1:411378/1‑149 (MQ=255)
                                                                                                                           tttGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATAGTTCAGCTTCACACGACGACCGCGTACCAGCGGCGGTTGGt                   <  2:142710/149‑1 (MQ=255)
                                                                                                                            ttGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACGACGACCGCGTACCCGCGGCggtgggtg                  >  2:502626/1‑149 (MQ=255)
                                                                                                                                ggAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCt              >  1:118010/1‑149 (MQ=255)
                                                                                                                                  aaGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGGTTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTc            <  2:660656/149‑1 (MQ=255)
                                                                                                                                         tcatcaAGTAGCGCTTGTAGGAAACAGGAAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACGACAACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGcc     >  1:491533/1‑149 (MQ=255)
                                                                                                                                            tgaaGTAGCGCTTGTAGGAATCGGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATc  <  2:215641/147‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTGCTTCATCAGACGCTTACGTTTACGCATCTGGGTTGGCGTCAGGGTGTTACGCTTATTCGCATACGGGTTTTCCCCTTCTTTGAACTGAATACGAATCGGCGATCCCATTACGTCCAGCGATTTGCGGAAGTAGTTCATCAAGTAGCGCTTGTAGGAATCAGGCAGGTCTTTCACCTGATTACCGTGAATCACCACAATCGGCGGGTTATAACCACCGGCGTGGGCATATTTCAGCTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATC  >  NC_000913/2635900‑2636188

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑