Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,173,661 G→A A29A (GCC→GCT)  U712_10870 ← Bifunctional protein BirA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,173,6610GA98.1% 157.8 / ‑4.5 53A29A (GCC→GCT) U712_10870Bifunctional protein BirA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (26/26);  total (26/27)
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.30e-01

AGTTTCCCGGGCTTTTTGATAAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAGCCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTGGAACAGCCGAGGGCATCACTGATTTTTTGGCCGGAAATAAATTCACTTCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGATTGGAGAGCTTTCC  >  CP006881/2173545‑2173794
                                                                                                                    |                                                                                                                                     
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)
                                          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                                                  <  7:3426/137‑1 (MQ=255)
                                                                 tttCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCa                                                  <  2:319228/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                                      <  6:11534/137‑1 (MQ=255)
                                                                             ttCAATATGCTTCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCa                                      <  1:222276/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                                                                     >  5:392448/1‑96 (MQ=255)
                                                                                       ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCt                                                                     >  5:398894/1‑96 (MQ=255)
                                                                                       ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCt                                                                     <  6:398894/96‑1 (MQ=255)
                                                                                       ttCCACACAGCAGTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCt                                                                     <  6:392448/96‑1 (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                                                                 >  5:179963/1‑90 (MQ=255)
                                                                                                 cagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTc                                                                 <  6:179963/90‑1 (MQ=255)
                                                                                                 cagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGa                >  1:115258/1‑139 (MQ=255)
                                                                                                 cagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGa                >  6:226415/1‑139 (MQ=255)
                                                                                                 cagTTCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGa                >  4:386832/1‑139 (MQ=255)
                                                                                                    ttCTTGAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGATTg             >  2:429698/1‑139 (MQ=255)
                                                                                                      cttgAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCTCATGTCTTCAGCCACTCCTCTATGTGTTTCTGTTGATTg             >  7:395223/4‑137 (MQ=255)
                                                                                                      cttgAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGATTg             >  3:259549/4‑137 (MQ=255)
                                                                                                         gAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGcc                                       <  4:84382/108‑1 (MQ=255)
                                                                                                         gAACAGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGcc                                       >  3:84382/1‑108 (MQ=255)
                                                                                                             aGCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGATTGGAGAGCttc    <  5:277644/139‑2 (MQ=255)
                                                                                                              gCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCa                                                  >  3:164186/1‑92 (MQ=255)
                                                                                                              gCCGAGAGCATCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCa                                                  <  4:164186/92‑1 (MQ=255)
                                                                                                               ccGAGAGCAGCACTGATTTTTTGGCCGGAAATAAATTCATTGCCGGCCTGAGAAAATAATTCAATAATGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGATTTGAGAGCTTccc  >  5:365691/1‑139 (MQ=255)
                                                                                                                    |                                                                                                                                     
AGTTTCCCGGGCTTTTTGATAAGCCGATATCCTTTTCTTCTAACGGCTTCTACTTCATAGCCCTCTTTCCGAAGCTCTTCAATATGCTTCCACACAGCAGTTCTGGAACAGCCGAGGGCATCACTGATTTTTTGGCCGGAAATAAATTCACTTCCGGCCTGAGAAAATAATTCAATAAGGTCTTTTCTTAATGTTGACCGCATGTCTTCAGCCACTCCTCTATGTGTTTCTTTTGATTGGAGAGCTTTCC  >  CP006881/2173545‑2173794

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