Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 130,787 T→A intergenic (+88/‑88) acnB → / → yacL bifunctional aconitate hydratase 2/2‑methylisocitrate dehydratase/UPF0231 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273130,7870TA90.6% 171.0 / 3.3 53intergenic (+88/‑88)acnB/yacLbifunctional aconitate hydratase 2/2‑methylisocitrate dehydratase/UPF0231 family protein
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base A (21/27);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGT  >  CP009273/130656‑130918
                                                                                                                                   |                                                                                                                                   
acacCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGc                                                                                                                                 <  1:1908352/136‑1 (MQ=255)
 cacCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGtt                                                                                                                             <  1:2125405/139‑1 (MQ=255)
    cGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAACCCATCGCCGTACCGCTGttttt                                                                                                                          <  1:1191801/139‑1 (MQ=255)
    cGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGtttt                                                                                                                           >  1:852102/1‑138 (MQ=255)
    cGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGtttt                                                                                                                           >  1:2513037/1‑138 (MQ=255)
    cGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGtttt                                                                                                                           >  1:2494006/1‑138 (MQ=255)
     gagaAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTtc                                                                                                                         <  1:444396/139‑1 (MQ=255)
     gagaAAGCCGAGGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGtttt                                                                                                                           >  1:2268256/1‑137 (MQ=255)
             cGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGt                                                                                                                   <  1:64673/137‑1 (MQ=255)
              gATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAg                                                                                                                >  1:244402/1‑139 (MQ=255)
                tGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAgg                                                                                                               <  1:1035557/138‑1 (MQ=255)
                    gTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGc                                                                                                           <  1:844498/138‑1 (MQ=255)
                     tGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCg                                                                                                          <  1:1668367/138‑1 (MQ=255)
                     tGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCg                                                                                                          <  1:350429/138‑1 (MQ=255)
                        ttttCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATa                                                                                                       <  1:1946160/138‑1 (MQ=255)
                          ttCCAGACTGCGGTTTAAAAAGTCAGCGCCCGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAAt                                                                                                     <  1:2338238/138‑1 (MQ=255)
                                 cTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACAt                                                                                                <  1:794756/136‑1 (MQ=255)
                                 cTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACAt                                                                                                <  1:2275854/136‑1 (MQ=255)
                                    cGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATg                                                                                           >  1:914600/1‑138 (MQ=255)
                                               gTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAg                                                                                                                <  1:1650813/106‑1 (MQ=255)
                                                 cAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGAtt                                                                              >  1:27892/1‑138 (MQ=255)
                                                         gcgcTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGACGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATaa                                                                                                      >  1:226371/1‑106 (MQ=255)
                                                                gcATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCgg                                                               <  1:2455687/138‑1 (MQ=255)
                                                                    aaaCGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCa                                                          <  1:1451136/139‑1 (MQ=255)
                                                                         acacAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGGTTTCGGTGCATTGCCGGGAGCAGAGGa                                                     >  1:2672057/1‑139 (MQ=255)
                                                                         acacAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGa                                                     >  1:101131/1‑139 (MQ=255)
                                                                         acacAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGa                                                     >  1:1689757/1‑139 (MQ=255)
                                                                         acacAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGa                                                     >  1:26450/1‑139 (MQ=255)
                                                                          cacaATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGaa                                                    <  1:1070613/139‑1 (MQ=255)
                                                                          cacaATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGa                                                     <  1:945944/138‑1 (MQ=255)
                                                                                  ccGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTa                                             <  1:1452574/138‑1 (MQ=255)
                                                                                  ccGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTa                                             <  1:1248957/138‑1 (MQ=255)
                                                                                      tGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTAcc                                                                                                                                   <  1:574832/48‑1 (MQ=255)
                                                                                      tGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCg                                                                                                          <  1:2509413/73‑1 (MQ=255)
                                                                                      tGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGa                                         <  1:895758/138‑1 (MQ=255)
                                                                                      tGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGa                                         <  1:1525720/138‑1 (MQ=255)
                                                                                                tCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTtc                                                                                                                         >  1:1998967/1‑48 (MQ=255)
                                                                                                tCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACAt                                                                                                >  1:1086044/1‑73 (MQ=255)
                                                                                                tCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAAtt                                >  1:1846482/1‑137 (MQ=255)
                                                                                                  ccGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtc                                                                                                                         <  1:1894788/46‑1 (MQ=255)
                                                                                                  ccGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCt                             >  1:2256757/1‑138 (MQ=255)
                                                                                                  ccGGGCATTTTTTTATTTCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGCAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCt                             >  1:1411265/1‑138 (MQ=255)
                                                                                                            ttttATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAgg                                                                                                               >  1:470554/1‑46 (MQ=255)
                                                                                                             tttATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCAGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTAc                 <  1:2263329/139‑1 (MQ=255)
                                                                                                                   tCTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGaa                                                    <  1:881615/98‑1 (MQ=255)
                                                                                                                    cTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTg          >  1:203132/1‑139 (MQ=255)
                                                                                                                    cTAAACCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTg          >  1:536607/1‑139 (MQ=255)
                                                                                                                      aaaCCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGt        <  1:1627607/139‑1 (MQ=255)
                                                                                                                       aaCCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTa       <  1:1960767/139‑1 (MQ=255)
                                                                                                                        aCCATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTaa      <  1:1919955/139‑1 (MQ=255)
                                                                                                                         ccATCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTgcg                          <  1:518666/118‑1 (MQ=255)
                                                                                                                            tCGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGt  >  1:495751/1‑139 (MQ=255)
                                                                                                                             cGCCGTACCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATgg                                          >  1:2305314/1‑98 (MQ=255)
                                                                                                                                   |                                                                                                                                   
ACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGT  >  CP009273/130656‑130918

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