Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,173,649 G→T S33S (TCC→TCA U712_10870 ← Bifunctional protein BirA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,173,6490GT98.2% 172.4 / ‑3.8 56S33S (TCC→TCAU712_10870Bifunctional protein BirA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (24/31);  total (24/32)
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.80e-01

TTCAATCCGAAGCGAATTTCACTTTCACTGAGTTTCCCGGGCTTTTTGATAAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAGCCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTGGAACAGCCGAGGGCATCACTGATTTTTTGGCCGGAAATAAATTCACTTCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGATTG  >  CP006881/2173515‑2173783
                                                                                                                                      |                                                                                                                                      
tataaTCCAAAACGAATTTCGCTTTCACTGAGTTTTCCGGGTTTTTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGaa                                                                                                                                     <  8:200186/135‑1 (MQ=255)
    aTCCAAAACGAATTTCGCTTTCACTGAGTTTTCCTGGTTTGTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTAATTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTCAACAGcc                                                                                                                                <  7:334692/139‑1 (MQ=255)
    aTCCAAAACGAATTTCGCTTTCACTGAGTTTTCCGGGTTTTTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGcc                                                                                                                                <  5:143087/139‑1 (MQ=255)
    aTCCAAAACGAATTTCGCTTTCACTGAGTTTTCCGGGTTTTTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGcc                                                                                                                                <  4:100291/139‑1 (MQ=255)
      ccaaaacGAATTTCGCTTTCACTGAGTTTTCCGGGTTTTTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTGTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGa                                                                                                                              <  3:43951/133‑1 (MQ=255)
                        tcactGAGTTTTCCGGGTTTTTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGaa                                                                                                                                     >  1:137129/1‑114 (MQ=255)
                        tcactGAGTTTTCCGGGTTTTTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGaa                                                                                                                                     <  2:137129/114‑1 (MQ=255)
                              aGTTTTCCGGGTTTTTTGATGATCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATCTGCTTCCACACAGCAGTTCTTTAACAGCCGAGAGCATAACTGATTTTTTGGCCg                                                                                                       <  4:207217/138‑1 (MQ=255)
                              aGTTTTCCGGGTTTTTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCg                                                                                                       <  3:378553/138‑1 (MQ=255)
                                   tccGGGTTTTTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGaaat                                                                                                  >  5:11534/2‑138 (MQ=255)
                                    ccGGGTTTTTTGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGACGGaaata                                                                                                 >  8:442120/1‑138 (MQ=255)
                                        gttttttGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGc                                                                                                         <  6:308341/124‑1 (MQ=255)
                                        gttttttGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGc                                                                                                         >  5:308341/3‑126 (MQ=255)
                                          tttttGATGAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCat                                                                                          >  6:390276/1‑138 (MQ=255)
                                                     ccGATATCCTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCt                                                                                 >  8:227906/1‑137 (MQ=255)
                                                             cTTTTCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAaataa                                                                       <  8:285150/139‑1 (MQ=255)
                                                               tttCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGaa                                                                                                                                     >  6:270966/1‑75 (MQ=255)
                                                               tttCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGaa                                                                                                                                     <  5:270966/75‑1 (MQ=255)
                                                               tttCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGcc                                                                                                                                >  7:257644/1‑80 (MQ=255)
                                                               tttCTTCTAACGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGcc                                                                                                                                <  8:257644/80‑1 (MQ=255)
                                                                        aCGGCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGt                                                            >  1:319228/1‑139 (MQ=255)
                                                                          ggCTTCTACTTCATAGCCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTGGAACAGCCGAGGGCATCACTGATTTTTTGGCCGGAAATAAATTCACTTCCGGCCTGAGAAAATAATTCAATAAGGTCt                                                          <  3:382274/139‑1 (MQ=255)
                                                                          ggCTTCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCt                                                          <  5:390276/139‑1 (MQ=255)
                                                                             ttCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAaataa                                                                       >  8:80501/1‑123 (MQ=255)
                                                                             ttCTACTTCATAACCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAaataa                                                                       <  7:80501/123‑1 (MQ=255)
                                                                                         acccTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGAc                                           <  7:206355/138‑1 (MQ=255)
                                                                                         accatctTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGTCCTGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGAc                                           <  8:126247/135‑1 (MQ=255)
                                                                                               tttCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCa                                       <  2:319228/137‑1 (MQ=255)
                                                                                               tttCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCa                                       <  7:3426/137‑1 (MQ=255)
                                                                                               tttCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCa                                       <  4:91896/137‑1 (MQ=255)
                                                                                               tttCCGAAGCTCTTCAATATGCTTACACACAGCAGTTATTGAACAGCAGAGAGCTTCACTGATTTTTTGGCAGGAAATAAATTAATTGCCGTCCTGAGAAAATAATTCAATAAGGTCTTTTATTAATGTTGAACGCATg                                     >  8:406845/1‑139 (MQ=255)
                                                                                                 tCCGAAGCTATTCAATATGCTTCCACACATCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGTAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCa                                       <  7:406845/135‑1 (MQ=255)
                                                                                                    gAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTc                                <  4:128616/139‑1 (MQ=255)
                                                                                                    gAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTc                                <  7:227906/139‑1 (MQ=255)
                                                                                                           ttCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCa                           <  1:222276/137‑1 (MQ=255)
                                                                                                           ttCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCa                           <  6:11534/137‑1 (MQ=255)
                                                                                                                 atGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTACTCTa                   <  4:79947/139‑1 (MQ=255)
                                                                                                                  tGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTAt                  >  2:376655/1‑139 (MQ=255)
                                                                                                                    cTTCCACCCAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTAtgt                >  6:277644/1‑139 (MQ=255)
                                                                                                                     ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTga                                                                               >  5:319742/1‑75 (MQ=255)
                                                                                                                     ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTga                                                                               <  6:319742/75‑1 (MQ=255)
                                                                                                                     ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCt                                                          <  6:392448/96‑1 (MQ=255)
                                                                                                                     ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCt                                                          <  6:398894/96‑1 (MQ=255)
                                                                                                                     ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCt                                                          >  5:398894/1‑96 (MQ=255)
                                                                                                                     ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCt                                                          >  5:392448/1‑96 (MQ=255)
                                                                                                                     ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGcc                            >  3:369216/1‑126 (MQ=255)
                                                                                                                     ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGcc                            <  4:369216/126‑1 (MQ=255)
                                                                                                                        cacacaGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGttt            >  2:86079/1‑139 (MQ=255)
                                                                                                                             agcagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATg                                                >  3:29076/1‑98 (MQ=255)
                                                                                                                             agcagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATg                                                <  4:29076/98‑1 (MQ=255)
                                                                                                                               cagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTc                                                      <  6:179963/90‑1 (MQ=255)
                                                                                                                               cagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTc                                                      >  5:179963/1‑90 (MQ=255)
                                                                                                                               cagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGa     >  6:226415/1‑139 (MQ=255)
                                                                                                                               cagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGa     >  1:115258/1‑139 (MQ=255)
                                                                                                                               cagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGa     >  4:386832/1‑139 (MQ=255)
                                                                                                                                  ttCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGATTg  >  2:429698/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                      
TTCAATCCGAAGCGAATTTCACTTTCACTGAGTTTCCCGGGCTTTTTGATAAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAGCCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTGGAACAGCCGAGGGCATCACTGATTTTTTGGCCGGAAATAAATTCACTTCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGATTG  >  CP006881/2173515‑2173783

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: