Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,940,742 G→C 8.2% L195V (CTG→GTG)  fliR ← Flagellar biosynthetic protein FliR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,940,7420GC 8.2% 146.1 / 7.0 49L195V (CTG→GTG) fliRFlagellar biosynthetic protein FliR
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base C (2/2);  total (24/25)
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.61e-02

CAAGACATCGGCCAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTAAAACTTTCGGT  >  ACB122/4940614‑4940871
                                                                                                                                |                                                                                                                                 
cAAGACATCGGCCAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCagcagc                                                                                                                               >  1:151053/1‑133 (MQ=255)
cAAGACATCGGCCAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCagcagc                                                                                                                               >  1:151054/1‑133 (MQ=255)
cAAGACATCGGCCAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCagcagc                                                                                                                               >  5:100925/1‑133 (MQ=255)
cAAGACATCGGCCAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCagcagc                                                                                                                               >  3:133237/1‑133 (MQ=255)
     cATCGGCCAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGa                                                                                                                        <  5:49373/135‑1 (MQ=255)
           ccAGGCCGACCCAGAAGATGACCATGCCCATGACCAAGGGCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCg                                                                                                                  <  1:178161/135‑1 (MQ=255)
           ccAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCg                                                                                                                  >  5:51857/1‑135 (MQ=255)
           ccAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCg                                                                                                                  >  6:309898/1‑135 (MQ=255)
           ccAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCACCAGCGCAGCGATCACCgg                                                                                                                 >  7:64686/1‑136 (MQ=255)
           ccAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCACCAGCGCAGCGATCACCgg                                                                                                                 >  7:58311/1‑136 (MQ=255)
            cAGGCCGACCCAGAAGATGAACATGCCCATGACAAAGGTCAGCGCGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCg                                                                                                                  >  4:248889/1‑134 (MQ=255)
             aGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCg                                                                                                                  >  2:77067/1‑133 (MQ=255)
                  gACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGACAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGat                                                                                                            >  1:134566/1‑133 (MQ=255)
                     ccAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCaa                                                                                                        <  5:230002/135‑1 (MQ=255)
                      cAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGCTCACCGGCAGAATCAAca                                                                                                      >  7:159792/1‑136 (MQ=255)
                          agaTGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAAcagca                                                                                                   <  3:351020/135‑1 (MQ=255)
                             ggAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGcc                                                                                               <  8:281051/135‑1 (MQ=255)
                              gAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGcc                                                                                               <  3:235889/135‑1 (MQ=255)
                                 cATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGc                                                                                            >  6:367782/1‑135 (MQ=255)
                                  aTGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCa                                                                                           >  4:36370/1‑135 (MQ=255)
                                       cATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGaa                                                                                      >  5:16124/1‑135 (MQ=255)
                                       cATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGaa                                                                                      >  2:149894/1‑135 (MQ=255)
                                        aTGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAaca                                                                                    >  8:174016/1‑136 (MQ=255)
                                              aaGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAg                                                                               >  3:150637/1‑135 (MQ=255)
                                                  tCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACCCCCAGCCCAt                                                                          >  8:199803/1‑136 (MQ=255)
                                                  tCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCa                                                                           >  4:209062/1‑135 (MQ=255)
                                                   cAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCAt                                                                          >  3:80838/1‑135 (MQ=255)
                                                   cAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCAt                                                                          >  4:10829/1‑135 (MQ=255)
                                                   cAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCAt                                                                          >  1:397441/1‑135 (MQ=255)
                                                       gggAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGcc                                                                     >  8:165844/1‑136 (MQ=255)
                                                         gAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGccc                                                                    <  5:101380/135‑1 (MQ=255)
                                                          aaGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCACCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCgg                                                                  <  8:64686/136‑1 (MQ=255)
                                                          aaGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCACCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCgg                                                                  <  8:58311/136‑1 (MQ=255)
                                                           aGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGACAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCgg                                                                  <  3:190087/135‑1 (MQ=255)
                                                           aGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGACAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCgg                                                                  <  3:196753/135‑1 (MQ=255)
                                                            gCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGc                                                                 <  4:150637/135‑1 (MQ=255)
                                                                aTGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGc                                                             <  6:16124/135‑1 (MQ=255)
                                                                   gagaAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTc                                                           >  4:200778/1‑134 (MQ=255)
                                                                              cAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGtt                                               <  7:291877/135‑1 (MQ=255)
                                                                              cAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAAGCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGtt                                               <  7:51804/135‑1 (MQ=255)
                                                                                 cTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGac                                            >  8:1287/1‑135 (MQ=255)
                                                                                     ggcgcggcgcGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGAccacca                                       <  7:165844/136‑1 (MQ=255)
                                                                                       cgcggcgcggGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCgg                                                                  >  8:73469/1‑107 (MQ=255)
                                                                                       cgcggcgcggGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCgg                                                                  <  7:73469/107‑1 (MQ=255)
                                                                                       cgcggcgcggGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAgc                                     <  7:174016/136‑1 (MQ=255)
                                                                                       cgcggcgcggGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAg                                      >  5:91866/1‑135 (MQ=255)
                                                                                                      tACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCgg                       <  1:160687/135‑1 (MQ=255)
                                                                                                                   aTGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTaaaa          <  6:91866/135‑1 (MQ=255)
                                                                                                                   aTGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTaaaa          <  6:100925/135‑1 (MQ=255)
                                                                                                                        caccacCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTAAAACTTTc     <  8:159792/135‑1 (MQ=255)
                                                                                                                         accaccAGCAGCGCAGCGATCACCGGCAGAATCAACCGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTAAAACTTTc     <  2:397441/134‑1 (MQ=255)
                                                                                                                          ccaccaGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTAAAACTTTCgg   <  3:209062/135‑1 (MQ=255)
                                                                                                                           caccaGCAGCGCAGAGATCACCGGCAGAATCAACAGCAGGCCGGCACCGACCACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTAAAACTTTCGGt  <  3:248889/135‑1 (MQ=255)
                                                                                                                              cagcagCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTAAAACTTTc     >  3:320804/1‑129 (MQ=255)
                                                                                                                              cagcagCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTAAAACTTTc     <  4:320804/129‑1 (MQ=255)
                                                                                                                                |                                                                                                                                 
CAAGACATCGGCCAGGCCGACCCAGAAGATGAACATGCCCATGACCAAGGTCAGCGGGAAGCCGATGGAGAAAATGTTCAGCTGCGGCGCGGCGCGGGTCATTACGCCAAAGGCAATGTTCACCACCAGCAGCGCAGCGATCACCGGCAGAATCAACAGCAGGCCGGCACCGAACACCCAGCCCATGCGCCCGGCCAGCTCCCAGAAGTGGTTGACCACCAGCGCACTCCCCACCGGCAAGGTGGTAAAACTTTCGGT  >  ACB122/4940614‑4940871

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: ‑