Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,328,359 T→C 20.8% N71S (AAC→AGC)  greA ← transcription elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,328,3590TC20.8% 103.0 / 25.7 53N71S (AAC→AGC) greAtranscription elongation factor
Reads supporting (aligned to +/- strand):  ref base T (23/19);  new base C (6/5);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

CGTTACGGTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCACGAGCTGCG  >  W3110S.gb/3328292‑3328424
                                                                   |                                                                 
cGTTACGGTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTc                                                                 <  1:3299325/70‑1 (MQ=255)
cGTTACGGTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTc                                                                 <  1:2749629/70‑1 (MQ=255)
 gTTACGGTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGa                                                               >  1:1734321/1‑71 (MQ=255)
    aCGGTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAg                                                            >  1:664656/1‑71 (MQ=255)
       gTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCtt                                                         <  1:1517250/71‑1 (MQ=255)
         aGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGa                                                               >  1:4578708/1‑63 (MQ=255)
         aGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGa                                                               >  1:3199649/1‑63 (MQ=255)
          gCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGc                                                      >  1:4297463/1‑71 (MQ=255)
             ccAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTc                                                   <  1:1669837/71‑1 (MQ=255)
               aaaaaTAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGa                                                 >  1:1579385/1‑71 (MQ=255)
                aaaaTAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGAt                                                <  1:2668055/71‑1 (MQ=255)
                  aaTAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCTCGTTCGACAGCTTGGCTTCGATg                                               >  1:4001338/1‑70 (MQ=255)
                   aTAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTc                                             <  1:4331016/71‑1 (MQ=255)
                       cgcgCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTc                                                                 >  1:3539101/1‑47 (MQ=255)
                       cgcgCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTa                                         <  1:2856202/71‑1 (MQ=255)
                        gcgcCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTaa                                        <  1:3216632/71‑1 (MQ=255)
                        gcgcCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTaa                                        <  1:3393572/71‑1 (MQ=255)
                         cgcCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTa                                         <  1:2051134/69‑1 (MQ=255)
                             cATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCtt                                                         <  1:1093155/49‑1 (MQ=255)
                             cATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCt                                            >  1:2825553/1‑62 (MQ=255)
                              aTTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGCTCGa                                                               <  1:1253069/42‑1 (MQ=255)
                               ttgttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTgg                                                       >  1:4311312/1‑49 (MQ=255)
                               ttgttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTc                                                   >  1:4338137/1‑53 (MQ=255)
                               ttgttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTc                                                   >  1:3273744/1‑53 (MQ=255)
                               ttgttgGGCATTTTGGTGACATCAATCACCTGCGCGCTCGACAGCtt                                                         <  1:2264425/47‑1 (MQ=255)
                                tgttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCt                                >  1:3256930/1‑71 (MQ=255)
                                 gttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGc                                                           >  1:2440248/1‑43 (MQ=255)
                                 gttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTgg                                                       >  1:3927719/1‑47 (MQ=255)
                                 gttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTgg                                                       >  1:2289339/1‑47 (MQ=255)
                                 gttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATg                                               >  1:3630914/1‑55 (MQ=255)
                                    gggCATTTTGGTGACATCAATCACCTGCGCGCTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCg                             <  1:2693327/70‑1 (MQ=255)
                                         ttttGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTc                                                   <  1:2669606/43‑1 (MQ=255)
                                         ttttGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGaaa                       <  1:997233/71‑1 (MQ=255)
                                           ttGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGa                                                 >  1:4009064/1‑43 (MQ=255)
                                           ttGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAAcc                     >  1:1636236/1‑71 (MQ=255)
                                           ttGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAAc                      <  1:2642275/70‑1 (MQ=255)
                                              gTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAAACctg                  <  1:2799834/71‑1 (MQ=255)
                                                  cATCAATCACCTGCGCGCTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGt              >  1:2667231/1‑71 (MQ=255)
                                                   atcaatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGtt             >  1:4428749/1‑71 (MQ=255)
                                                   atcaatcaCCTGCGCGCTCGACAGCTTGGCTTCGATGTCTTTAATACGGGCTTCGCAGAAACCCTGCTGtt             >  1:1906662/1‑71 (MQ=255)
                                                   atcaatcaCCTGCGCGCTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGtt             >  1:4174132/1‑71 (MQ=255)
                                                   atcaatcaCCTGCGCGCTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGtt             >  1:4273094/1‑71 (MQ=255)
                                                   atcaatcaCCTGCGCGCTCGACAGCTTGGCTGCGATGTCttt                                          <  1:4026544/42‑1 (MQ=37)
                                                     caatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCa           >  1:702435/1‑71 (MQ=255)
                                                      aatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATAc                                     >  1:773308/1‑44 (MQ=255)
                                                      aatcaCCTGCGCGCTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCAc          >  1:1675646/1‑71 (MQ=255)
                                                       atcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCttt                                          <  1:3514150/38‑1 (MQ=255)
                                                       atcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCACg         <  1:2903715/71‑1 (MQ=255)
                                                        tcaCCTGCGCGCTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCACg         <  1:2365722/70‑1 (MQ=255)
                                                            cTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCACGAGCTg    <  1:2990003/71‑1 (MQ=255)
                                                            cTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCACGAGCTg    <  1:1440013/71‑1 (MQ=255)
                                                            cTGCGCGCTCGACAGCTTGGCTTCGATGTCTTTAATa                                      >  1:1260680/1‑37 (MQ=255)
                                                              gcgcGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCACGAGCTGCg  >  1:2425807/1‑71 (MQ=255)
                                                                   |                                                                 
CGTTACGGTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCACGAGCTGCG  >  W3110S.gb/3328292‑3328424

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