Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,942,694 A→G 8.5% L162P (CTG→CCG)  surA_2 ← Chaperone SurA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,942,6940AG 8.5% 135.8 / 4.3 47L162P (CTG→CCG) surA_2Chaperone SurA
Reads supporting (aligned to +/- strand):  ref base A (17/26);  new base G (2/2);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

CGCACGGATTGGCGCGCCTGCGCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCT  >  ACB122/3942580‑3942821
                                                                                                                  |                                                                                                                                
cgcACGGATTGGCGCGCCTGCGCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCGGTCTGTCCTGGTTGGCATc                                                                                                              >  6:65167/1‑134 (MQ=255)
cgcACGGATTGGCGCGCCTGCGCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCGGGCTGTCCTGGTTGGCATc                                                                                                              >  6:65169/1‑134 (MQ=255)
         tGGCGCGCCTGCGCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGcc                                                                                                    >  5:293991/1‑135 (MQ=255)
         tGGCGCGCCTGCGCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGcc                                                                                                    >  3:93278/1‑135 (MQ=255)
             gcgcctgcgcctCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCa                                                                                                <  1:171399/135‑1 (MQ=255)
                 ctgcgcctCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTc                                                                                            <  7:66011/135‑1 (MQ=255)
                 ctgcgcctCCCGCTGCCCGATGAGCGGTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTc                                                                                            <  8:58413/135‑1 (MQ=255)
                   gcgcctCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGc                                                                                          >  5:233611/1‑135 (MQ=255)
                        tCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCctcg                                                                                     <  3:23970/135‑1 (MQ=255)
                        tCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTGCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCctcg                                                                                     <  5:127171/135‑1 (MQ=255)
                             cTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCgg                                                                                <  3:30617/135‑1 (MQ=255)
                             cTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCgg                                                                                <  5:240292/135‑1 (MQ=255)
                                  cGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCGGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGt                                                                           <  5:65167/135‑1 (MQ=255)
                                  cGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGt                                                                           <  6:293991/135‑1 (MQ=255)
                                  cGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGt                                                                           <  4:93278/135‑1 (MQ=255)
                                  cGAGGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCGGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGt                                                                           <  5:65169/135‑1 (MQ=255)
                                   gATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTc                                                                          >  5:67278/1‑135 (MQ=255)
                                   gATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTc                                                                          >  5:60522/1‑135 (MQ=255)
                                    aTGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTGCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCGCGGCGCGTCg                                                                         <  6:128665/135‑1 (MQ=255)
                                       agggTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGTCCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCTTCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACt                                                                      <  1:285777/132‑1 (MQ=255)
                                        gCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACt                                                                      <  3:288828/134‑1 (MQ=255)
                                        gCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACt                                                                      <  4:250004/134‑1 (MQ=255)
                                         cGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTgc                                                                    <  8:20227/135‑1 (MQ=255)
                                              cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTg                                                                <  7:67265/134‑1 (MQ=255)
                                              cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTg                                                                <  6:233611/134‑1 (MQ=255)
                                                     cAAGGCCAGGCCGTGTTCGCGTATGACCAATGTATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTGGGCATCAAAAAATGCCTCCACTTCGCCCTCGCGCGGCGCGTCGACTGCGCTGAGCAGGGc                                                        <  1:13792/135‑1 (MQ=255)
                                                      aaGGCCAGGCCGTGTTCGCGTATGTCCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGc                                                        >  4:49965/1‑134 (MQ=255)
                                                           caggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTacac                                                  >  4:65655/1‑135 (MQ=255)
                                                                cGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGCGCAGGGCGTACACTTGCt                                             >  4:219977/1‑135 (MQ=255)
                                                                cGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCt                                             >  8:229857/1‑135 (MQ=255)
                                                                cGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCt                                             >  7:65665/1‑135 (MQ=255)
                                                                     tCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgcc                                        >  8:283806/1‑135 (MQ=255)
                                                                     tCGCGTATGACCAATGGATTAGCCCTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCAGCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGGCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgc                                         <  4:37248/134‑1 (MQ=255)
                                                                               ccATTGGATTATCACTTTGTTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGAATAAAAAAATGCCTCCACTTCCCCCTCGCTCGGCGCGTCGACTGCGCTGAGCTGGGCGTACACTTGCTGGGCCGCCATCTCa                               <  2:83073/134‑1 (MQ=255)
                                                                               ccACTGGATTATCCCTTTGTTGCTGCGCTCCTTGCAGTGGGTCGGGTGGGCCAAAAAAAAATCCCGCCCCTTCCCCCTGGCGCGGCGCGGCGCCGGGGCTGCGCCGGGCGTACACTTGCTGGGCCGCCATCTCa                               <  6:164734/134‑1 (MQ=255)
                                                                               ccAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCAt                              >  8:132608/1‑135 (MQ=255)
                                                                                 aaTGTATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTGGGAATCAAAAAATGCCTCCACTCCGCCCGCGCGCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGc                            <  3:49965/135‑1 (MQ=255)
                                                                                         aTCAATTTGGTTCTGCGGGCGTTGAAGTCTGGCGTGGGTGGCATCAAAAACTGCCTCCAGTTCGCCCTCGCGGGGCGGGGCTACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgtgt                     <  6:60522/134‑1 (MQ=255)
                                                                                          tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCa                               >  3:212359/1‑123 (MQ=255)
                                                                                          tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGAATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCa                               <  4:212359/123‑1 (MQ=255)
                                                                                           cACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTGGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGGCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTAt                   <  5:144560/134‑1 (MQ=255)
                                                                                                tGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCctcg                                                                                     >  7:301348/1‑63 (MQ=255)
                                                                                                tGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCctcg                                                                                     <  8:301348/63‑1 (MQ=255)
                                                                                                   ttCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCa                               <  3:242482/114‑1 (MQ=255)
                                                                                                   ttCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCa                               >  4:242482/1‑114 (MQ=255)
                                                                                                   ttCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTAc          >  3:124279/1‑135 (MQ=255)
                                                                                                           tCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCt  >  4:123204/1‑135 (MQ=255)
                                                                                                                  |                                                                                                                                
CGCACGGATTGGCGCGCCTGCGCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCT  >  ACB122/3942580‑3942821

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: