Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,196,981 C→A 8.7% P170P (CCG→CCT)  bacA ← undecaprenyl pyrophosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,196,9810CA 8.7% 137.7 / 3.4 46P170P (CCG→CCT) bacAundecaprenyl pyrophosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base A (2/2);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

AGGAAGCCCCAGCTTTTGTAGAGATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGCAACAAACCGCCAACGACCAGCGCATACATCACA  >  CP009273/3196835‑3197128
                                                                                                                                                  |                                                                                                                                                   
aGGAAGCCCCAGCTTTTGTAGAGATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCgg                                                                                                                                                   >  1:481329/1‑149 (MQ=255)
aGGAAGCCCCAGCTTTTGTAGAGATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCgg                                                                                                                                                   >  1:309395/1‑149 (MQ=255)
       cccAGCTTTTGTAGAGATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGGGAAACCAGGc                                                                                                                                                  >  2:435660/1‑143 (MQ=255)
       cccAGCTTTTGTAGAGATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCAGGc                                                                                                                                                  >  2:435761/1‑143 (MQ=255)
         cAGCTTTGGTAGAGATCGAGCGCCGTTGCGCCCATCATCGTCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAgcgc                                                                                                                                          <  1:370322/149‑1 (MQ=255)
            cTTTTGTAGAGATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAg                                                                                                                                             <  2:186309/143‑1 (MQ=255)
            cTTTTGTAGAGATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAg                                                                                                                                       >  1:369726/1‑149 (MQ=255)
             ttttGTAGAGATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAccgcaaaa                                                                                                                                      >  1:119346/1‑145 (MQ=255)
                   taggaTCGGGCCCCGTTGCGCCCATCAGCAACGGCAACGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGa                                                                                                                                      <  2:378670/140‑1 (MQ=255)
                   agagATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGa                                                                                                                                      >  2:152254/1‑143 (MQ=255)
                       aTCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAAc                                                                                                                            >  1:453575/1‑149 (MQ=255)
                                  ttGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCa                                                                                                                       >  2:44792/1‑143 (MQ=255)
                                  tcgcgcCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCCCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATa                                                                                                                 <  1:272714/147‑1 (MQ=255)
                                       cccATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCGTAAATGc                                                                                                            >  1:411945/1‑149 (MQ=255)
                                       cccATCATCATCCGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGCAACCCGGCCACAGCGCCAGACCCTGGACACAGCCCATCATAAATGc                                                                                                            >  1:381369/1‑149 (MQ=255)
                                          atcatcatcGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCAaaaa                                                                                                               <  2:158039/143‑1 (MQ=255)
                                          atcatcatcGGCACCACCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATaaa                                                                                                               >  2:288404/1‑143 (MQ=255)
                                           tcatcatcGGCACCCCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGa                                                                                                        <  1:190140/149‑1 (MQ=255)
                                            catcatcGGCACCGCCAGCAGGAACGAAACCTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATg                                                                                                             <  2:152191/143‑1 (MQ=255)
                                            catcatcGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATg                                                                                                             <  2:363262/143‑1 (MQ=255)
                                              tcatcGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGcc                                                                                                           >  2:389671/1‑143 (MQ=255)
                                                 tcGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAACCAGCCAGTCATAAATGCCTGACGATAg                                                                                                  >  1:63332/1‑149 (MQ=255)
                                                 tcGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGa                                                                                                        >  2:389689/1‑143 (MQ=255)
                                                    gCACCGCCAGCACGTACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGa                                                                                                     <  2:236520/143‑1 (MQ=255)
                                                       ccGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAg                                                                                                  <  2:446373/143‑1 (MQ=255)
                                                              cagGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATACCATCaa                                                                                     >  1:117813/1‑149 (MQ=255)
                                                               agGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCACCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATatc                                                                                          <  2:319202/143‑1 (MQ=255)
                                                                    cGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCg                                                                               >  1:98068/1‑149 (MQ=255)
                                                                              gAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCgg                                                                     <  1:433651/149‑1 (MQ=255)
                                                                                            gCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCg                                                             <  2:246853/143‑1 (MQ=255)
                                                                                             cctaCCCCGCTCAGCCTCCCACCTACACTGCTAACGTCGGGATGGGAGATTCCCGGCCACAGCGCCAGACACTGGAAACAGCCAAGCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCgg                                                            <  2:49264/140‑1 (MQ=255)
                                                                                                   cccATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAg                                                      <  2:106738/143‑1 (MQ=255)
                                                                                                       tCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCAt                                                  >  2:58815/1‑143 (MQ=255)
                                                                                                             tCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCCGGCATTCGGCGGCAATc                                      >  1:287892/1‑149 (MQ=255)
                                                                                                               ccACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAg                                    <  1:171626/149‑1 (MQ=255)
                                                                                                                  cccgAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGcaa                                 <  1:50574/146‑1 (MQ=255)
                                                                                                                  ccTGAAATGGTCGCCCCGGAACGGGAGAAACCAGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGcaa                                 <  1:483013/149‑1 (MQ=255)
                                                                                                                  ccTGAAATGGTCGCCCCGGAACGGGAGAAACCAGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGcaa                                 <  1:483115/149‑1 (MQ=255)
                                                                                                                    tGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCa                                     <  2:409362/143‑1 (MQ=255)
                                                                                                                    ggAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCa                                     <  2:409476/142‑1 (MQ=255)
                                                                                                                           gTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGCAACAAACCGcc                        <  1:433670/149‑1 (MQ=255)
                                                                                                                                 ccGGGACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGCAACAAACCGCCAACGAc                  <  1:168605/149‑1 (MQ=255)
                                                                                                                                 ccGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAAACAGCAACAAACCGACAACGAc                  >  1:5790/1‑149 (MQ=255)
                                                                                                                                  cGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGCAACAAACCGCCAACGAcc                 <  1:311981/149‑1 (MQ=255)
                                                                                                                                   ggAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGCAACAAACCGCCaa                      <  2:434082/143‑1 (MQ=255)
                                                                                                                                       cGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGCAACAAACCGCCAACGACCAgcgc            <  1:291787/149‑1 (MQ=255)
                                                                                                                                                 ccGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGCAACAAACCGCCAACGACCAGCGCATACATcaca  <  1:470659/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGGAAGCCCCAGCTTTTGTAGAGATCGAGCGCCGTTGCGCCCATCATCATCGGCACCGCCAGCAGGAACGAAAACTCGGAAGCAGCGTAACGGCTCACCCCCATCAGCATCCCACCTGAAATGGTCGCCCCGGAACGGGAGAAACCCGGCCACAGCGCCAGACACTGGAAACAGCCAATCATAAATGCCTGACGATAGGTCATATCATCAAGACCCGGCGCACGCGGCTCTTTCGGCTTCAGGCATTCGGCGGCAATCAGCAACAAACCGCCAACGACCAGCGCATACATCACA  >  CP009273/3196835‑3197128

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