Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,636,594 T→C 7.6% K402E (AAG→GAG)  lon_1 ← Lon protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,636,5940TC 7.6% 159.4 / 2.7 52K402E (AAG→GAG) lon_1Lon protease
Reads supporting (aligned to +/- strand):  ref base T (15/33);  new base C (2/2);  total (17/35)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCCTGCACCAGCTTGCCGGGCTGGGCGCCGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGAGGCCGGAGTGG  >  ACB122/1636465‑1636724
                                                                                                                                 |                                                                                                                                  
gggCCTGCACCAGCTTGCCGGGCTGGGCGCCGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCg                                                                                                                                <  2:224281/134‑1 (MQ=255)
    cTGCACCAGCTTGCCGGGCTGGGCGCCGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCt                                                                                                                           <  4:139702/135‑1 (MQ=255)
        aCCAGCTTGCCGGGCTGGGCGCCGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTc                                                                                                                       <  8:63230/135‑1 (MQ=255)
         ccAGCTTGCCGGGCTGGGCGCCGATATAGGGGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCt                                                                                                                      <  7:100199/135‑1 (MQ=255)
          ccgCTTGCCGGGCTGGGCTCCGATATAGGGGCTGCGGTGTCCCTTGATCTCGGCTTCATCACGCATGCCGCCGCCGCTGAAGCGATCGCACGGCCGGCCGAGCGATTCGGCGATGGACTTGCAGATGCTGGTCtt                                                                                                                     <  1:165055/133‑1 (MQ=255)
          cAGCTTGCCGGGCTGGGCGCCGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCtt                                                                                                                     <  3:263489/135‑1 (MQ=255)
          cAGCTTGCCGGGCTGGGCGCCGATATAGGTGCGGCGGGGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCtt                                                                                                                     <  7:29394/135‑1 (MQ=255)
          cAGCGTGCCGGGGTGGGCGCCGATATAGGTGCGGCGGTGGCCGTTGATCTCGGCTTCATCTCGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCtt                                                                                                                     <  3:243626/135‑1 (MQ=255)
                    ggCTGGGCGCCGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGAGGGCCTTGCCGATGCTGGTCTTGCCCGCCCcg                                                                                                           >  7:110057/1‑135 (MQ=255)
                         ggCGCCGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCcgggcgg                                                                                                      >  8:249486/1‑135 (MQ=255)
                             ccGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGccc                                                                                                  <  6:244974/135‑1 (MQ=255)
                             ccGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGccc                                                                                                  <  2:288166/135‑1 (MQ=255)
                             ccGATATAGGTGCGGCGGGGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGccc                                                                                                  <  4:6269/135‑1 (MQ=255)
                                         cggcggTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCAc                                                                                       >  2:48011/1‑134 (MQ=255)
                                            cggTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTCGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCCCCAACAGCCCGaac                                                                                   >  8:143700/1‑133 (MQ=255)
                                            cggTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTCGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATc                                                                                   >  8:143698/1‑135 (MQ=255)
                                              gTGGCGCTTGATCTCGGCTTCATCACGCATGCCGCCGCCGCTGAAGCGATAGAACGGACGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCACCAGCACGCTCGa                                                                                 >  6:203217/1‑135 (MQ=255)
                                                    cTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGc                                                                            >  6:103990/1‑134 (MQ=255)
                                                    cTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGc                                                                            >  4:213873/1‑134 (MQ=255)
                                                     ttGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTg                                                                          <  4:220585/135‑1 (MQ=255)
                                                                ttCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCcgg                                                                                                          >  6:38220/1‑92 (MQ=255)
                                                                ttCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCcgg                                                                                                          <  5:38220/92‑1 (MQ=255)
                                                                ttCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCttt                                                               <  1:185963/135‑1 (MQ=255)
                                                                ttCATCACGCATGCCGCCGACGCTGAACCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCttt                                                               <  3:245141/135‑1 (MQ=255)
                                                                          aTGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCTACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGccc                                                     >  8:207655/1‑135 (MQ=255)
                                                                          aTGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCACCAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGccc                                                     >  6:167747/1‑135 (MQ=255)
                                                                                ccgACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCCCCAACAGCACGCTCGAGCCGCTGGTCTCGCCGTTCCCGGCGCCCCCcgcg                                               >  5:184131/1‑135 (MQ=255)
                                                                                ccgACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACgg                                                 >  5:293103/1‑133 (MQ=255)
                                                                                          aGCGATTGAAGGCCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCGTTTCCAGGCGCCCACGGCGAGGAATTCg                                      <  1:274591/134‑1 (MQ=255)
                                                                                            cGATAGAACGGCCGTCCGAGCGATTCGGCGATGGACTTGCCGATGCGGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGCTCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAgg                                   <  7:45101/135‑1 (MQ=255)
                                                                                              aTAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGAt                                 <  1:48011/135‑1 (MQ=255)
                                                                                              aTAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGGCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGAt                                 <  5:248995/135‑1 (MQ=255)
                                                                                                          gccgAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGAt                     >  1:299502/1‑135 (MQ=255)
                                                                                                          gccgAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGAt                     >  1:214177/1‑135 (MQ=255)
                                                                                                            cgAGGGATTCGGCGATGGACTCGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATAGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATgt                   <  7:143700/135‑1 (MQ=255)
                                                                                                            cgAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGc                                                >  5:115223/1‑106 (MQ=255)
                                                                                                            cgAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGc                                                <  6:115223/106‑1 (MQ=255)
                                                                                                            cgAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATgt                   <  8:169820/135‑1 (MQ=255)
                                                                                                            cgAGCGATTCGGCGATGGACTCGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATgt                   <  7:143698/135‑1 (MQ=255)
                                                                                                                  agtCGGCGATGGACTTGCCGCTGCTGGTCTTGCCCCCCCCGGGCGGGCCCACCAACCGCACGAGCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGa             <  1:133532/133‑1 (MQ=255)
                                                                                                                  aTTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGa             <  6:90693/135‑1 (MQ=255)
                                                                                                                  aTTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGa             <  5:103990/135‑1 (MQ=255)
                                                                                                                  aTTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGa             <  1:2434/135‑1 (MQ=255)
                                                                                                                  aTTCGGCGATGGACGTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGa             <  3:58380/135‑1 (MQ=255)
                                                                                                                    tCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGAgg           <  3:202783/135‑1 (MQ=255)
                                                                                                                     cGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGAgg           <  2:299502/134‑1 (MQ=255)
                                                                                                                       gCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGAt                                 >  4:36123/1‑110 (MQ=255)
                                                                                                                       gCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGAt                                 <  3:36123/110‑1 (MQ=255)
                                                                                                                              aCTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGAGGCCGGAGTgg  <  6:274114/134‑1 (MQ=255)
                                                                                                                              aCTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGAGGCCGGAGTgg  <  4:187639/134‑1 (MQ=255)
                                                                                                                              aCTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGAGGCCGGAGTgg  <  7:213280/134‑1 (MQ=255)
                                                                                                                              aCTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGAGGCCGGAGTgg  <  2:214177/134‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                  
GGGCCTGCACCAGCTTGCCGGGCTGGGCGCCGATATAGGTGCGGCGGTGGCCCTTGATCTCGGCTTCATCACGCATGCCGCCGACGCTGAAGCGATAGAACGGCCGGCCGAGCGATTCGGCGATGGACTTGCCGATGCTGGTCTTGCCCACCCCGGGCGGGCCCACCAACAGCACGATCGAGCCGCTGATCTCGCCTTTCCAGGCGCCCACGGCGAGGAATTCGAGGATGCGCTCCTTGATGTCGTCGAGGCCGGAGTGG  >  ACB122/1636465‑1636724

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