Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,405,614 A→T 9.1% I23N (ATT→AAT)  PP_1228 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,405,6140AT 9.1% 125.0 / 3.8 44I23N (ATT→AAT) PP_1228methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base T (2/2);  total (22/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 = 8.93e-01

GCGCTGGATACGCAAGGCCTGGGTTTCGGCGTGGGCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTGGCAGCGTGGGTCATTCGGCGACGT  >  NC_002947/1405469‑1405761
                                                                                                                                                     |                                                                                                                                                   
gcgcTGGATACGCAAGGCCTGGGTTTCGGCGTGGGCCTGCAGGCGCTGGCGGGCGCTGTCTTCCAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATg                                                                                                                                                   <  6:98678/148‑1 (MQ=255)
      gATACGCAAGGCCTGGGTTTCGGCGTGGGCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCcagca                                                                                                                                            <  3:217056/149‑1 (MQ=255)
                gCCTGGGTTTCGGCGTGGGCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGcc                                                                                                                                  >  2:283109/1‑149 (MQ=255)
                  cTGGGTTTCGGCGTGGGCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCgg                                                                                                                                <  3:170432/149‑1 (MQ=255)
                  cTGGGTTTCGGCGTGGGCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCgg                                                                                                                                >  4:170432/1‑149 (MQ=255)
                           ggCGTGGGCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCTTGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCATCAGGCACAGGCCGGCGAGGAgt                                                                                                                        <  7:77645/148‑1 (MQ=255)
                                gggCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGc                                                                                                                                   >  1:209655/1‑132 (MQ=255)
                                gggCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGc                                                                                                                                   <  2:209655/132‑1 (MQ=255)
                                 ggCGTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCACCGAA‑‑‑GTGGCGACCAGCAAGGCGACATTGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCt                                                                                                                 >  7:154899/1‑149 (MQ=255)
                                 ggCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACATTGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCt                                                                                                                 >  7:154897/1‑149 (MQ=255)
                                 ggCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCt                                                                                                                 >  5:25067/1‑149 (MQ=255)
                                      gcgggCGCTGGGGGGCGCTGTCTTCGATCATGGCAGTGTTGGCCTGTTTTCTCACGGGCGGCGCTGCGATGCGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCa                                                                                                             <  1:101655/146‑1 (MQ=255)
                                       cAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCa                                                                                                             <  1:283109/147‑1 (MQ=255)
                                            gcTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACtggatg                                                                                                      <  1:53925/149‑1 (MQ=255)
                                              tGGGGGGCGCTGTCTTAGATCATGGCAGTGTTGGCCTGGTTGTCCA‑TGGGGGCGGTGCGATGGGCGTGGGTGATCTTATCCCTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACtggat                                                                                                       <  1:98442/149‑1 (MQ=255)
                                                gcgggcggTGTCTTCGAGCTTTGCAGTGTTGCGCTTGTTGACCA‑GGGCGTCGCTGAGATTGGCGTGGTTGAGCGAA‑‑‑GTGGCGACCAGCTAGGGGTCATTGGCCAGCAGGCAGAGGCCGGGGTGGAGTGTGATCTTCCACtggatggtg                                                                                                   <  8:154899/148‑3 (MQ=255)
                                                gcgggcgCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACATTGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGgag                                                                                                   <  8:154897/148‑1 (MQ=255)
                                                gcgggcgCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGgag                                                                                                   <  4:28864/148‑1 (MQ=255)
                                                        tGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACtggat                                                                                                       >  5:122366/1‑136 (MQ=255)
                                                        tGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACtggat                                                                                                       <  6:122366/136‑1 (MQ=255)
                                                            ttCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCcagc                                                                                                                                             >  2:88195/1‑94 (MQ=255)
                                                            ttCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCcagc                                                                                                                                             <  1:88195/94‑1 (MQ=255)
                                                             tCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGa                                                                                                                    <  3:135694/118‑1 (MQ=255)
                                                             tCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGa                                                                                                                    >  4:135694/1‑118 (MQ=255)
                                                                            ttGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCt                                                                      >  3:213477/1‑149 (MQ=255)
                                                                                  tGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGc                                                                                       >  1:195098/1‑126 (MQ=255)
                                                                                  tGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGc                                                                                       <  2:195098/126‑1 (MQ=255)
                                                                                         ccA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAAt                                                          >  5:92683/1‑148 (MQ=255)
                                                                                          cA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCtttt                                                                   <  3:161177/138‑1 (MQ=255)
                                                                                          cA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCtttt                                                                   >  4:161177/1‑138 (MQ=255)
                                                                                                 ggcgCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGcc                                                  >  2:183943/1‑149 (MQ=255)
                                                                                                                     tCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAAt                                                          >  2:270478/1‑121 (MQ=255)
                                                                                                                       aGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCgttgt                            >  4:35714/1‑149 (MQ=255)
                                                                                                                         cGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGttgttg                          <  3:35714/149‑1 (MQ=255)
                                                                                                                           aa‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTgg                         >  7:67523/1‑148 (MQ=255)
                                                                                                                            a‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTGGCa                       <  1:204383/149‑1 (MQ=255)
                                                                                                                                 tGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGCGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTGGCAGc                     >  4:97985/1‑149 (MQ=255)
                                                                                                                                     gACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTc                                                                  <  2:94968/100‑1 (MQ=255)
                                                                                                                                       ccAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTGGCAGCGTGGGt               <  2:134365/149‑1 (MQ=255)
                                                                                                                                        cAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTc                                                                  >  1:149780/1‑97 (MQ=255)
                                                                                                                                        cAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTc                                                                  <  2:149780/97‑1 (MQ=255)
                                                                                                                                          gCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTGGCAGCGTGGGTCAt            >  5:164777/1‑149 (MQ=255)
                                                                                                                                                 gACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTGGCAGCGTGGGTCATTCGGCGa     >  3:61666/1‑149 (MQ=255)
                                                                                                                                                 gACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTGGCAGCGTGGGTCATTCGGCGa     >  6:87353/1‑149 (MQ=255)
                                                                                                                                                    aaTGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCtttt                                                                   <  1:216763/84‑1 (MQ=255)
                                                                                                                                                    aaTGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCtttt                                                                   >  2:216763/1‑84 (MQ=255)
                                                                                                                                                    aaTGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTGGCAGCGTGGGTCATTCGGCGACGt  <  5:87353/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGCTGGATACGCAAGGCCTGGGTTTCGGCGTGGGCCTGCAGGCGCTGGCGGGCGCTGTCTTCGAGCATGGCAGTGTTGGCCTGGTTGACCA‑GGGCGGCGCTGCGATGGGCCTGGGTCAGCGAA‑‑‑GTGGCGACCAGCAAGGCGACAATGGCCAGCAGGCACAGGCCGGCGAGGAGTGTGATCTTCCACTGGATGGAGAGGCGACGCAGCGGCATGGGGGCTTTCCTTTTCGAACAAATGCAACGCCCGCGATCAGGCGCGGGCGTTGTTGGCAGCGTGGGTCATTCGGCGACGT  >  NC_002947/1405469‑1405761

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