Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,251,818 C→T A315T (GCG→ACG)  nusA ← transcription elongation factor NusA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,251,8180CT91.2% 109.0 / ‑0.5 34A315T (GCG→ACG) nusAtranscription elongation factor NusA
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (11/20);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTG  >  REL606/3251750‑3251889
                                                                    |                                                                       
gTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTa                                                                  >  1:940669/1‑76 (MQ=255)
gTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTa                                                                  <  1:295215/76‑1 (MQ=255)
   gAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCg                                                               >  1:1576537/1‑76 (MQ=255)
       tCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCtt                                                           <  1:1426896/76‑1 (MQ=255)
        cccAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCCTc                                                          >  1:180802/1‑76 (MQ=255)
         ccAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCa                                                         <  2:1273813/76‑1 (MQ=255)
         ccAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCa                                                         <  1:246326/76‑1 (MQ=255)
         ccAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCa                                                         <  1:526906/76‑1 (MQ=255)
               gCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACNGGCTTCAACggag                                                   <  1:1078472/76‑3 (MQ=255)
               gCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCg                                                   <  1:813284/76‑1 (MQ=255)
                cTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCCTCAACGGCGa                                                  >  2:833483/1‑76 (MQ=255)
                 tCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGat                                                 <  1:1046418/76‑1 (MQ=255)
                   aGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGGGCCAGGCTACCGGCTTCAACGGGGatat                                               >  1:1517224/1‑76 (MQ=255)
                     cTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATcc                                             <  1:422724/76‑1 (MQ=255)
                          aaGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTTCACGGCGATATCCATggg                                        >  2:261229/1‑75 (MQ=255)
                          aaGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGgt                                        <  2:986768/76‑1 (MQ=255)
                                  cGCACGTCCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTAt                                <  2:421680/76‑1 (MQ=255)
                                   gCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATc                               >  2:916356/1‑76 (MQ=255)
                                    cACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTTAACGGCGAGATCCATGGGGGGGTTATCt                              >  2:544220/1‑76 (MQ=255)
                                     aCGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCtt                             <  1:1121940/76‑1 (MQ=255)
                                     aCGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCtt                             >  1:630190/1‑76 (MQ=255)
                                       gTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCa                           <  1:1426620/76‑1 (MQ=255)
                                       gTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCCACGGCGGTATCCCTGGTGTGTTTATCttta                           >  2:360106/1‑74 (MQ=255)
                                          cTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATcc                        <  2:1465413/76‑1 (MQ=255)
                                            gACCGTTACGGCCAGTCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATccac                      >  1:863436/1‑76 (MQ=255)
                                            gACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATccac                      >  2:752563/1‑76 (MQ=255)
                                                  tACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGAt                <  1:1015642/76‑1 (MQ=255)
                                                     ggCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGa             <  1:57217/76‑1 (MQ=255)
                                                      gCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGaa            <  1:1590854/76‑1 (MQ=255)
                                                       ccAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAg           <  2:764649/76‑1 (MQ=255)
                                                         aaTCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCa         >  1:749179/1‑76 (MQ=255)
                                                         aaTCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCa         >  1:596307/1‑76 (MQ=255)
                                                              ccTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTc    <  1:341910/76‑1 (MQ=255)
                                                               cTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCt   <  1:255291/76‑1 (MQ=255)
                                                                tGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTg  <  1:1555623/76‑1 (MQ=255)
                                                                cgcgtCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTg  <  2:416072/75‑1 (MQ=255)
                                                                    |                                                                       
GTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTG  >  REL606/3251750‑3251889

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