Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,128,645 C→A 7.5% D82Y (GAT→TAT)  wcaF ← putative acyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,128,6450CA 7.5% 185.5 / 3.2 53D82Y (GAT→TAT) wcaFputative acyl transferase
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base A (2/2);  total (30/23)
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.52e-01

CAATCACAATAGGCGTGGCGTTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACT  >  NC_000913/2128507‑2128778
                                                                                                                                          |                                                                                                                                     
cAATCACAATAGGCGTGGCGTTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCg                                                                                                                                      >  2:780661/1‑140 (MQ=255)
cAATCACAATAGGCGTGGCGTTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCg                                                                                                                                      >  2:283566/1‑140 (MQ=255)
              gtggcgtTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTa                                                                                                                        <  1:721003/140‑1 (MQ=255)
              gtggcgtTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTa                                                                                                                        <  1:720688/140‑1 (MQ=255)
                 gcgtTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATc                                                                                                                     >  1:33980/1‑140 (MQ=255)
                 gcgtTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATCTGCGCCAATGGTTATTACACCGAGGGTATATAAAATGACGTCAACGCCGACCCACGCGTAATc                                                                                                                     >  1:34019/1‑140 (MQ=255)
                            gAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTa                                                                                                          >  1:332676/1‑140 (MQ=255)
                            gAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTa                                                                                                          >  1:48499/1‑140 (MQ=255)
                            gAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTa                                                                                                          >  1:474522/1‑140 (MQ=255)
                                   tGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCa                                                                                                   <  1:445356/140‑1 (MQ=255)
                                      cTTGCATGGTCGTGGCTACCTGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACgg                                                                                                <  2:132574/140‑1 (MQ=255)
                                          cATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAg                                                                                            >  1:1086743/1‑140 (MQ=255)
                                          cATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAg                                                                                            >  2:1144710/1‑140 (MQ=255)
                                          cATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAg                                                                                            >  1:35717/1‑140 (MQ=255)
                                               tCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAAt                                                                                       >  1:452308/1‑140 (MQ=255)
                                               tCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAAt                                                                                       >  1:452284/1‑140 (MQ=255)
                                                    gCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATAGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTa                                                                                  >  2:542682/1‑140 (MQ=255)
                                                    gCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATAGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTa                                                                                  >  2:628908/1‑140 (MQ=255)
                                                        ccGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTga                                                                              <  1:539174/140‑1 (MQ=255)
                                                         cGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTgac                                                                             <  2:205459/140‑1 (MQ=255)
                                                         cGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTgac                                                                             <  2:202683/140‑1 (MQ=255)
                                                         cGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTgac                                                                             <  1:824821/140‑1 (MQ=255)
                                                         cGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTgac                                                                             <  1:824727/140‑1 (MQ=255)
                                                         cGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTgac                                                                             <  1:1035939/140‑1 (MQ=255)
                                                                   ataaCTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATaa                                                                   >  1:615658/1‑140 (MQ=255)
                                                                   ataaCTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATaa                                                                   >  1:615580/1‑140 (MQ=255)
                                                                     aaCTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACt                                                                 >  2:960104/1‑140 (MQ=255)
                                                                        ttttttGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACg                                                              >  2:845531/1‑140 (MQ=255)
                                                                        ttttttGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACg                                                              >  2:674073/1‑140 (MQ=255)
                                                                                gATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTcc                                                      <  2:1018796/140‑1 (MQ=255)
                                                                                gATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTcc                                                      <  2:762400/140‑1 (MQ=255)
                                                                                          aaTGTGCGCCAATTGTTATTTCACCTAGGGAATATAAATTGACGTCATCGCCGACCCACGAGTAATCACCTAAGTTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTATTTCCTATTTTTGCt                                            <  2:596537/140‑1 (MQ=255)
                                                                                              tgCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGa                                        >  2:478734/1‑140 (MQ=255)
                                                                                                cgcCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAAt                                      >  1:889012/1‑140 (MQ=255)
                                                                                                     aTGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGACTAATCg                                 >  1:72462/1‑140 (MQ=255)
                                                                                                                  cGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGcc                    >  1:79778/1‑140 (MQ=255)
                                                                                                                  cGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGcc                    >  2:507324/1‑140 (MQ=255)
                                                                                                                  cGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGCTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGcc                    >  1:79805/1‑140 (MQ=255)
                                                                                                                        tatataAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAg              <  1:414391/140‑1 (MQ=255)
                                                                                                                        tatataAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAg              <  1:341139/140‑1 (MQ=255)
                                                                                                                        tatataAATTGACGTCATAGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAg              <  1:542682/140‑1 (MQ=255)
                                                                                                                        tatataAATTGACGTCATAGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAg              <  1:628908/140‑1 (MQ=255)
                                                                                                                           ttaATTTGACGTCTTCGCCGTCCCTCGCGTTATCACCTTTGGTTTTTTTCCTCGGATAGGTAATTTTTAGTGACGGACGTATAACTACGTTTTTTCCTATTTTTGCTCCGAATTAACGTAATAAAAAAGCCCGCCAGCGa           <  2:599790/139‑1 (MQ=255)
                                                                                                                           ataAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGa           <  1:540648/140‑1 (MQ=255)
                                                                                                                           ataAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGa           <  1:540675/140‑1 (MQ=255)
                                                                                                                           ataAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGATAATTTCCACGTATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGa           <  2:336226/140‑1 (MQ=255)
                                                                                                                              aaTTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATAc        <  2:248746/140‑1 (MQ=255)
                                                                                                                              aaTTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAAGTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATAc        <  1:453743/140‑1 (MQ=255)
                                                                                                                               aTTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACa       >  2:143406/1‑140 (MQ=255)
                                                                                                                               aTTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACa       >  2:143294/1‑140 (MQ=255)
                                                                                                                               aTTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACa       >  2:143133/1‑140 (MQ=255)
                                                                                                                                    cGTCATCGCCGACCCACTCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCATATTTTTGCTCAGAATAAACGTAATAAAAAATCCCGCCAGCGATACAATACt  >  2:593045/1‑140 (MQ=255)
                                                                                                                                    cGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGACCGCCAGCGATACAAAACt  >  2:592564/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                     
CAATCACAATAGGCGTGGCGTTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAATTTCCACGGATAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACT  >  NC_000913/2128507‑2128778

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