Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,157,831 G→A 15.0% A165T (GCG→ACG)  holB → DNA polymerase III, delta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,157,8310GA15.0% 103.7 / 9.9 40A165T (GCG→ACG) holBDNA polymerase III, delta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (17/17);  new base A (2/4);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCA  >  W3110S.gb/1157763‑1157896
                                                                    |                                                                 
tCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTgcg                                                                 >  1:582582/1‑71 (MQ=255)
                      aCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCg                                           >  1:998028/1‑71 (MQ=255)
                      aCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCg                                           >  1:3501146/1‑71 (MQ=255)
                          ttACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGAc                                       >  1:4012151/1‑71 (MQ=255)
                           tACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTg                                         <  1:4034881/68‑1 (MQ=255)
                           tACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGAcc                                      >  1:1028860/1‑71 (MQ=255)
                            aCTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCt                                     >  1:494759/1‑71 (MQ=255)
                             cTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTACGCCGCCGCCGGAACAGTACGCCGTGACCTg                                    >  1:5012840/1‑71 (MQ=255)
                                gCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCt                                 >  1:2334165/1‑71 (MQ=255)
                                   aCATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCttt                               <  1:2601232/70‑1 (MQ=255)
                                        aCGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGaa                                                     >  1:1134105/1‑43 (MQ=255)
                                         cGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGc                                             >  1:4606661/1‑50 (MQ=255)
                                         cGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCg                                           >  1:340853/1‑52 (MQ=255)
                                         cGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCAcgcg                         <  1:2044117/70‑1 (MQ=255)
                                         cGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCAcgcg                         <  1:1239252/70‑1 (MQ=255)
                                         cGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCAcgcg                         <  1:1269899/70‑1 (MQ=255)
                                         cGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCAcgcg                         <  1:4006749/70‑1 (MQ=255)
                                                gTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGAc                  >  1:2048562/1‑70 (MQ=255)
                                                gTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGAc                  <  1:612575/70‑1 (MQ=255)
                                                gTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGAc                  >  1:1202461/1‑70 (MQ=255)
                                                  tGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAAt               >  1:1817180/1‑71 (MQ=255)
                                                    tCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGt             <  1:331944/71‑1 (MQ=255)
                                                      ggTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTca           <  1:2265357/71‑1 (MQ=255)
                                                          aCATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAg        <  1:3433021/70‑1 (MQ=255)
                                                          aCATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAg        <  1:1232730/70‑1 (MQ=255)
                                                          aCATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAg        <  1:394411/70‑1 (MQ=255)
                                                          aCATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAg        >  1:3835970/1‑70 (MQ=255)
                                                          aCATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAg        <  1:1464866/70‑1 (MQ=255)
                                                            aTTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAg        <  1:3028222/68‑1 (MQ=255)
                                                            aTTACCTTACGCCGCCGCCGGAACAGTACGCCGTGACCTGGCttt                               <  1:247274/45‑1 (MQ=255)
                                                            aTTACCTTACGCCGCCGCCGGAACAGTACGCCGTGACCTGGCttt                               <  1:3847792/45‑1 (MQ=255)
                                                            aTTACCTTACGCCGCCGCCGGAACAGTACGCCGTGACCTGGCttt                               <  1:2713061/45‑1 (MQ=255)
                                                            aTTACCTTACGCCGCCGCCGGAACAGTACGCCGTGACCTGGCttt                               <  1:2702706/45‑1 (MQ=255)
                                                             ttACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATg    >  1:2360906/1‑71 (MQ=255)
                                                             ttACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATg    >  1:4147814/1‑71 (MQ=255)
                                                              tACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAg                      >  1:1358619/1‑52 (MQ=255)
                                                               aCCTTACGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGa      >  1:5053798/1‑67 (MQ=255)
                                                                 cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCttt                               <  1:1788979/40‑1 (MQ=255)
                                                                 cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCttt                               <  1:1139219/40‑1 (MQ=255)
                                                                  ttGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATg    <  1:4801406/66‑1 (MQ=255)
                                                                   tGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTcaca         >  1:1618097/1‑60 (MQ=255)
                                                                    gcgcCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCa  >  1:1298649/1‑66 (MQ=255)
                                                                    gcgcCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCa  >  1:2600754/1‑66 (MQ=255)
                                                                    gcgcCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCa  >  1:3035841/1‑66 (MQ=255)
                                                                    |                                                                 
TCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCA  >  W3110S.gb/1157763‑1157896

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