Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,064,259 T→G V93G (GTG→GGG)  BW25113_RS10440 → GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,064,2590TG95.8% 69.8 / ‑3.3 24V93G (GTG→GGG) BW25113_RS10440GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (11/12);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCAC  >  NZ_CP009273/2064172‑2064345
                                                                                       |                                                                                      
tACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGc                                                                                      <  1:223004/90‑1 (MQ=255)
    gtgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGtt                                                                                  >  1:244394/1‑90 (MQ=255)
    gtgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGtt                                                                                  >  1:208945/1‑90 (MQ=255)
                    gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                                  >  2:137920/1‑90 (MQ=255)
                     aTGACCGTGCCCTGTATGTGGATGAGAATTTCTGGCGACAAATCCTTCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcca                                                                 >  1:482150/1‑89 (MQ=255)
                              cccTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACgg                                                        <  2:459247/90‑1 (MQ=255)
                                 tgtctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGc                                                     <  2:208945/90‑1 (MQ=255)
                                        ggATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCa                                              <  2:244394/90‑1 (MQ=255)
                                            gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                          >  2:460913/1‑90 (MQ=255)
                                            gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                          >  1:47268/1‑90 (MQ=255)
                                              gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATg                                        >  1:262355/1‑90 (MQ=255)
                                                   tCTGGCGACACATCCTGCAGTGCGGACCTCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGata                                   <  2:137362/90‑1 (MQ=255)
                                                        cGACACATCCTGCANTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGcc                              <  2:361299/90‑1 (MQ=255)
                                                        cGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGcc                              <  2:403737/90‑1 (MQ=255)
                                                               tCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTc                       <  1:365872/90‑1 (MQ=255)
                                                                    cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt                  <  2:61575/90‑1 (MQ=255)
                                                                      gTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGAc                                                               >  2:258023/1‑43 (MQ=255)
                                                                      gTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGAc                                                               <  1:258023/43‑1 (MQ=255)
                                                                        gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                                  <  1:538691/38‑1 (MQ=255)
                                                                        gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                                  >  2:538691/1‑38 (MQ=255)
                                                                        gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTtctc              >  1:543675/1‑90 (MQ=255)
                                                                        gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTtctc              >  1:434236/1‑90 (MQ=255)
                                                                               tCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcaca       <  2:371802/90‑1 (MQ=255)
                                                                                    agGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcacaggcac  <  2:169290/90‑1 (MQ=255)
                                                                                       |                                                                                      
TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCAC  >  NZ_CP009273/2064172‑2064345

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