Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,964,793 G→A A150V (GCG→GTG)  tap ← methyl‑accepting chemotaxis protein IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,964,7930GA84.0% 71.0 / 4.9 26A150V (GCG→GTG) tapmethyl‑accepting chemotaxis protein IV
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (14/7);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ATCATGCTGATAAACACCAGTGCCGAAATCTGATAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCGCAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATCCCCTGACGTGCTCGTTGCCCGCAGTCATCGCCATAA  >  NZ_CP009273/1964655‑1964931
                                                                                                                                          |                                                                                                                                          
atCATGCTGATAAACACCAGTGCCGAAATCTGATAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGgcg                                                                                                                                            >  1:123844/1‑139 (MQ=255)
               acCAGTGCCGAAATCTGATAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATcc                                                                                                                              >  1:217388/1‑138 (MQ=255)
                   gTGCCGAAATCTGATAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAg                                                                                                                         <  2:160821/139‑1 (MQ=255)
                        gAAATCTGATAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCGCAGTGAGGAAATCCGAAAGTTGAt                                                                                                                    <  2:123844/139‑1 (MQ=255)
                         aaaTCTGATAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGAtt                                                                                                                   <  1:656073/139‑1 (MQ=255)
                         aaaTCTGATAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGAt                                                                                                                    >  2:783139/1‑138 (MQ=255)
                                aTAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACttt                                                                                                              <  2:217388/137‑1 (MQ=255)
                                aTAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTc                                                                                                             >  1:17089/1‑138 (MQ=255)
                                            gCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGcc                                                                                                >  2:746718/1‑139 (MQ=255)
                                            gCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGcc                                                                                                <  2:26637/139‑1 (MQ=255)
                                                       tGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTc                                                                                        >  2:725508/1‑136 (MQ=255)
                                                               tCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTtgt                                                                             >  1:236744/1‑139 (MQ=255)
                                                                       aTGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCgg                                                                     >  2:621264/1‑139 (MQ=255)
                                                                               tctcCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACt                                                             >  1:97078/1‑139 (MQ=255)
                                                                                     gCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTtctc                                                       >  2:581957/1‑139 (MQ=255)
                                                                                         ccAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGtc                                                   >  1:699951/1‑139 (MQ=255)
                                                                                         ccAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGtc                                                   >  1:856756/1‑139 (MQ=255)
                                                                                           aGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGtctc                                                 <  2:236744/139‑1 (MQ=255)
                                                                                                   aaGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCa                                         <  1:783139/139‑1 (MQ=255)
                                                                                                         aCGTCAAACGCATTCTGTGATCCCTGCACCGGCGCAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATcccc                                   >  2:460280/1‑139 (MQ=255)
                                                                                                               aaCGCATTCTGTGATCCCTGCACCGGCGCAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATCCCCTGACGt                             <  1:460280/139‑1 (MQ=255)
                                                                                                                     ttCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCa                                         >  2:506501/1‑121 (MQ=255)
                                                                                                                     ttCTGTGATCCCTGCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCa                                         <  1:506501/121‑1 (MQ=255)
                                                                                                                            aTCCCTGCACCGGCGCAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATCCCCTGACGTGCTCGTTGCCCg                 <  1:834419/138‑1 (MQ=255)
                                                                                                                                 tGCACCGGCCCAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATCACCTGACGTGCTCGTTGCCCGCAGTCa           >  2:148937/1‑139 (MQ=255)
                                                                                                                                  gCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATCCCCTGACGTGCTCGTTGCCCGCAGTCAt          >  2:244891/1‑139 (MQ=255)
                                                                                                                                  gCACCGGCACAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATCCCCTGACGTGCTCGTTGCCCGCAGTCAt          >  2:59461/1‑139 (MQ=255)
                                                                                                                                        gcgcAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCa                                         <  2:721378/102‑1 (MQ=255)
                                                                                                                                        gcgcAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCa                                         >  1:721378/1‑102 (MQ=255)
                                                                                                                                        gcgcAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATCCCCTGACGTGCTCGTTGCCCGCAGTCATCGCCa     >  2:360761/1‑138 (MQ=255)
                                                                                                                                          gcAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATCCCCTGACGTGCTCGTTGCCCGCAGTCATCGCCATaa  <  2:430959/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
ATCATGCTGATAAACACCAGTGCCGAAATCTGATAGTTACGCTGGCTTTGCGCACTGGCGGCTTCCAGCACATGGTTGATCTCCAGCTGCCAGGCCTCAAAGTTAACGTCAAACGCATTCTGTGATCCCTGCACCGGCGCAGTGAGGAAATCCGAAAGTTGATTACTTTCAAGCCAGGTCGCCTGATGTTCGAGATCGTTGTGCCAGCGGGCAAAACTTTTCTCCGTCTCTTTTTGCAATCCCCTGACGTGCTCGTTGCCCGCAGTCATCGCCATAA  >  NZ_CP009273/1964655‑1964931

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