Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,310,843 A→G 11.4% I172V (ATT→GTT)  gabP‑III → gamma‑aminobutyrate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,310,8430AG11.4% 89.2 / 3.4 35I172V (ATT→GTT) gabP‑IIIgamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base A (14/17);  new base G (2/2);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TGGTGATCACGTTGCTGCTGACGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGTTCTCGTTCCTCGGGGCCGAGGTGGTGACCATTG  >  NC_002947/3310715‑3310990
                                                                                                                                |                                                                                                                                                    
tGGTGATCACGTTGCTGCTGACGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCggggg                                                                                                                                 <  2:208197/150‑1 (MQ=255)
  gTGATCACGTTGCTGCTGACGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCGTTTTCGGCCTGTTGCCGGGGGcc                                                                                                                               >  2:363754/1‑150 (MQ=255)
  gTGATCACGTTGCTGCTGACGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCGTTTTCGGCCTGTTGCCGGGGGcc                                                                                                                               >  2:192243/1‑150 (MQ=255)
         cGTTGCTGCTGACGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCa                                                                                                                       <  1:147442/150‑1 (MQ=255)
               tgctgACGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAgcggcg                                                                                                                 >  1:759485/1‑150 (MQ=255)
                 ctgACGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGtc                                                                                                               >  1:798890/1‑150 (MQ=255)
                   gACGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGtctc                                                                                                             <  1:7228/150‑1 (MQ=255)
                    aCGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTcg                                                                                                            <  2:745846/150‑1 (MQ=255)
                          accaACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGTCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCtg                                                                                                      >  2:803210/1‑150 (MQ=255)
                          accaACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCtg                                                                                                      >  2:92282/1‑150 (MQ=255)
                             aaCCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTCTTCTCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTgg                                                                                                   <  1:225568/150‑1 (MQ=255)
                                  gTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTcgc                                                                                                           <  1:319096/137‑1 (MQ=255)
                                  gTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTcgc                                                                                                           >  2:319096/1‑137 (MQ=255)
                                           gaagaaTTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCtt                                                                                     >  1:247590/1‑150 (MQ=255)
                                                  tACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTAATCAACAACGTCACGTTCATCGCCCTGGGTGTTCGTGTCGTTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGccc                                                                              <  1:192243/150‑1 (MQ=255)
                                                  tACGGTGAGTTCGAGTTTTGGCTCGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCGTTTTCGGCCTGTTGCCTCAACCGACTGGCGTCGGTGTCGCGCGCCTGTGGGATACGGGGGGCTTCATGccc                                                                              <  1:363754/150‑1 (MQ=255)
                                                           ttCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACg                                                                                             >  1:146113/1‑126 (MQ=255)
                                                           ttCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACg                                                                                             <  2:146113/126‑1 (MQ=255)
                                                           ttCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCttt                                                                     <  2:798890/150‑1 (MQ=255)
                                                                            aCTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCGGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCttt                                                                     <  2:331598/133‑1 (MQ=255)
                                                                            aCTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCttt                                                                     >  1:331598/1‑133 (MQ=255)
                                                                            aCTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAgcgc                                                    >  2:375601/1‑150 (MQ=255)
                                                                                gTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGACATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGa                                                  <  1:66462/148‑1 (MQ=255)
                                                                                gTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGACATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGa                                                  >  2:66462/1‑148 (MQ=255)
                                                                                 tCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTggtgcggtg                                                            >  1:804137/1‑137 (MQ=255)
                                                                                 tCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTggtgcggtg                                                            <  2:804137/137‑1 (MQ=255)
                                                                                        ggtggCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCAc                                        <  1:767984/150‑1 (MQ=255)
                                                                                              cATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGtt                                  >  2:239280/1‑150 (MQ=255)
                                                                                                   tCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGTTCTCGt                             <  1:803210/150‑1 (MQ=255)
                                                                                                          cATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGTTCTCGTTCCTCgg                      <  2:537233/150‑1 (MQ=255)
                                                                                                           aTCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGTTCTCGTTCCTCggg                     >  2:302188/1‑150 (MQ=255)
                                                                                                                  tGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGTTCTCGTTCCTCGGGGCCGAgg              >  2:447156/1‑150 (MQ=255)
                                                                                                                    ggTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGTTCTCGTTCCTCGGGGCCGAggtg            <  1:239280/150‑1 (MQ=255)
                                                                                                                         ttGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGTTCTCGTTCCTCGGGGCCGAGGTGGTGAc       <  1:204972/150‑1 (MQ=255)
                                                                                                                              ccATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGTTCTCGTTCCTCGGGGCCGAGGTGGTGACCATtg  <  1:447156/150‑1 (MQ=255)
                                                                                                                                |                                                                                                                                                    
TGGTGATCACGTTGCTGCTGACGGCCACCAACCTGTTCAGCGTGAAGAATTACGGTGAGTTCGAGTTTTGGCTGGCACTGGTCAAGGTGGTGGCCATCGTCAGCTTCATCGCCCTGGGTGTTTGCGCCATTTTCGGCCTGTTGCCGGGGGCCGGCGTCAGCGGCGTCTCGCGCCTGTGGGATACGGGGGGCTTCATGCCCAACGGCTTTGGTGCGGTGCTGAGCGCGATGCTGATCACCATGTTCTCGTTCCTCGGGGCCGAGGTGGTGACCATTG  >  NC_002947/3310715‑3310990

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