Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,938,227 T→G 32.4% R533R (CGT→CGG)  pqiB → Intermembrane transport protein PqiB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,938,2270TG32.4% 33.3 / 3.9 34R533R (CGT→CGG) pqiBIntermembrane transport protein PqiB
Reads supporting (aligned to +/- strand):  ref base T (15/7);  new base G (3/9);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 2.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAACTCGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTACC  >  ACB122/2938093‑2938359
                                                                                                                                       |                                                                                                                                     
gCCAACTCGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGt                                                                                                                                       <  1:124053/135‑1 (MQ=255)
gCCAACTCGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGt                                                                                                                                       <  2:336705/135‑1 (MQ=255)
  cAACTCGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTgg                                                                                                                                     >  4:175594/1‑135 (MQ=255)
  cAACTCGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTgg                                                                                                                                     >  2:86399/1‑135 (MQ=255)
  cAACTCGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCGGCTTCGTgg                                                                                                                                     >  8:360347/1‑135 (MQ=255)
    aCTCGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTggcc                                                                                                                                   <  3:298638/135‑1 (MQ=255)
     cTCGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGTGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGAGTCGTggccc                                                                                                                                  >  2:75201/1‑134 (MQ=255)
       cGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCCCCCCGCGTCGCTGAATAGTggccgcc                                                                                                                                >  6:369173/1‑133 (MQ=255)
       cGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGAGTCGTggccggc                                                                                                                                >  2:215986/1‑135 (MQ=255)
                       ttCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTc                                                                                                                >  3:332147/1‑135 (MQ=255)
                          gccgcAGCGCGAGCTAATGGGTGAACCGCTGGCCGAGCTGGCGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATAACCTGGG‑CCGCCACCCCGAGGCGCTGATGCGTGGCCGGCCTAAAGCGGCCGGGTCCg                                                                                                              <  4:332147/134‑1 (MQ=255)
                               aGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCAGCCGGGTCCGCAGAcc                                                                                                        <  8:128570/135‑1 (MQ=255)
                                        gaATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCGGTTTCGGGGCCCGCCGAAAGCGGCCGGGTCCGCCGGCCGTGGAcccg                                                                                               >  1:263442/1‑134 (MQ=255)
                                              gTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCACCACCCCGAGTCGCTGATTCCGGGCCGGGCGCAAGCGGCCGGCGCCGCAGACCGTGGCCCCTcaggag                                                                                         >  8:341631/1‑135 (MQ=255)
                                                 aaccgcgGGACGAGGGGGAGCGCATGGCGCGTGCGCTCCGGGCCCTTGCCGATACCCTGGG‑CCGCCCCCCCGAGTCGCTGATGCGTGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACa                                                                                     <  4:241472/132‑1 (MQ=255)
                                                 aaaCGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGGTTAGTTGCCCGCCGCAAGCGGGCGGGTGCGCCGGCCGTGGACCCCGAGGAG‑ACa                                                                                     >  6:51854/1‑135 (MQ=255)
                                                           cGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAACGGCCGGGTGCGCAGACCATGAACCCGGAGGAG‑GCACCCGATGaaa                                                                           >  1:107172/1‑135 (MQ=255)
                                                            gAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCCGGTCCGCAGACCATGAACCCGGAGGAG‑CCAGCCGATGCAAc                                                                          >  1:64745/1‑135 (MQ=255)
                                                            gAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCGGCTTCGTGGCCGGCCGAAACCGGCCGGGGGCGCCGCCCGTGGACCCCGAGGGG‑ACAGCCGGGGCAAc                                                                          >  7:161839/1‑135 (MQ=255)
                                                                tGGCGCGCATGTCGCGTTCGCTACGTGCCCTTGCCGATTACCGGGC‑CCGCCCCCCCGAGGCGCTGATGCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACt                                                                      <  2:290569/135‑1 (MQ=255)
                                                                gggtgcgcAGGTCGCGTTCGCACCGTGACCTTGCGTATTCCGGGGG‑CCGCCACCCCGAGGCGCTGATCCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACt                                                                      <  3:269177/131‑1 (MQ=255)
                                                                    gcgcATGTCGCGTTCGCTACGTGACCTTTCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAATCGGCCCGGTCCGCAGACCTTGCAACCAGAGGAG‑CCACCCGCTTCAACCACTgcg                                                                   >  8:201910/1‑134 (MQ=255)
                                                                           tCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGCTTCGTGGCCGGCCGCAAGCGGCCGGGTCCGCAGCCCTTGAACCCGGAGGAG‑CCAGCCGCTGCAACGCCTGCGCAATGCAt                                                           >  4:295995/1‑135 (MQ=255)
                                                                               gTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCCCCCCGCGTCGCGGCTTCGGGGCCGGCCGCAAGCCGCCGGGGGCGCCGCCCGTGCACCCGGAGGCGCCCCGCCG‑GGCCACGCCGGCGCCATGCAtgtgt                                                       >  5:136785/1‑135 (MQ=255)
                                                                                            ccTTGCCGATTCCCGGGC‑CCGCCACCCCGCGGCGCTGATGCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCctgctg                                          <  2:18737/135‑1 (MQ=255)
                                                                                                      tACCCGGGCCCGCCACCCCGCGGCGCTGATGGGGGGCCGGCCGCACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTc                                 <  8:52898/135‑1 (MQ=255)
                                                                                                      aaCCTGGC‑CCCCCACCCCGATGCCCTGATGCGTGCCCGCCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTc                                 <  6:257960/133‑1 (MQ=255)
                                                                                                          gggg‑CCGCCCCCCCGCGGCGCTGATGCGGGCCCGCCCTAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCggggg                            <  6:176613/134‑1 (MQ=255)
                                                                                                          gggg‑CCCCCCCCCCGCGTCCCGGATAGGGGGCCGCCCGACCGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCggggg                            <  4:234125/134‑1 (MQ=255)
                                                                                                          gggg‑CCCCCCCCCCGCGGCGCTGATGCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCggggg                            <  5:193340/134‑1 (MQ=255)
                                                                                                                   ccccccGCGCCGCTGAGGCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACCGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTgcagca                    <  4:33537/133‑1 (MQ=255)
                                                                                                                     ccccGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTACAGCAGt                  >  1:201325/1‑135 (MQ=255)
                                                                                                                           ggcggtGATGCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCg            <  7:360347/130‑1 (MQ=255)
                                                                                                                           gTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCACACCTTTAACCCGGACGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCCGCAGTACGcc             >  4:290535/1‑134 (MQ=255)
                                                                                                                                   ttCGTGGCCGGCCTAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTa    <  3:175594/135‑1 (MQ=255)
                                                                                                                                      gTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTAcc  <  8:202247/134‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GCCAACTCGGCCATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGG‑CCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAG‑ACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTACC  >  ACB122/2938093‑2938359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑