Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,252,012 G→A A250V (GCT→GTT)  nusA ← transcription elongation factor NusA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,252,0120GA97.7% 160.8 / ‑3.1 44A250V (GCT→GTT) nusAtranscription elongation factor NusA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (22/21);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

GACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCG  >  REL606/3251939‑3252084
                                                                         |                                                                        
gACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACa                                                                        <  1:1360110/76‑1 (MQ=255)
    atatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTa                                                                    <  2:221203/76‑1 (MQ=255)
    atatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTa                                                                    <  2:523791/76‑1 (MQ=255)
        cGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCgg                                                                >  1:1756851/1‑76 (MQ=255)
          aTACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGAt                                                              <  2:1776841/76‑1 (MQ=255)
           tACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATc                                                             <  1:1401309/76‑1 (MQ=255)
            aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg                                                            <  1:1160488/76‑1 (MQ=255)
            aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg                                                            <  2:1832063/76‑1 (MQ=255)
            aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg                                                            >  2:1441360/1‑76 (MQ=255)
            aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg                                                            >  2:515734/1‑76 (MQ=255)
            aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg                                                            <  2:613514/76‑1 (MQ=255)
            aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg                                                            >  1:426005/1‑76 (MQ=255)
               cTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATa                                                         >  2:1356376/1‑76 (MQ=255)
                 cGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATCCCTACGCAAACACCTACCGGATCGATCCg                                                       >  2:969126/1‑76 (MQ=255)
                 cGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACg                                                       <  2:1570988/76‑1 (MQ=255)
                     aCCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTa                                                   >  2:899767/1‑76 (MQ=255)
                     aCCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTa                                                   >  2:354775/1‑76 (MQ=255)
                      cccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTAt                                                  <  2:987226/76‑1 (MQ=255)
                       ccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATc                                                 <  2:434010/76‑1 (MQ=255)
                        cAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCg                                                >  1:700564/1‑76 (MQ=255)
                         aGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGt                                               >  2:1646868/1‑76 (MQ=255)
                            tCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTgg                                            <  2:998990/76‑1 (MQ=255)
                            tCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTgg                                            >  2:1693729/1‑76 (MQ=255)
                            tCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTgg                                            >  2:1567998/1‑76 (MQ=255)
                                tAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGtttt                                        <  2:281308/76‑1 (MQ=255)
                                 aGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTc                                       <  2:904732/76‑1 (MQ=255)
                                  gACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCa                                      <  1:937833/76‑1 (MQ=255)
                                      ccGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcgc                                  >  2:885637/1‑76 (MQ=255)
                                       cGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcgcg                                 >  2:819288/1‑76 (MQ=255)
                                          cTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                              >  2:146847/1‑76 (MQ=255)
                                          cTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                              >  1:1985969/1‑76 (MQ=255)
                                           tGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAttt                             >  2:194518/1‑76 (MQ=255)
                                             aaCACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc                           <  1:1640524/76‑1 (MQ=255)
                                             aaCACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc                           >  2:265278/1‑76 (MQ=255)
                                             aaCACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc                           <  2:827364/76‑1 (MQ=255)
                                                      gcCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGaa                  <  1:48093/76‑1 (MQ=255)
                                                        cACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAcc                >  1:1303223/1‑76 (MQ=255)
                                                          cGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCg              >  2:706872/1‑76 (MQ=255)
                                                          cGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCg              <  2:1442598/76‑1 (MQ=255)
                                                             aTACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGAt           <  2:1195101/76‑1 (MQ=255)
                                                             aTACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGAt           <  2:1262384/76‑1 (MQ=255)
                                                                ccTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgc        <  1:885637/76‑1 (MQ=255)
                                                                  tACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGa      >  2:1218250/1‑76 (MQ=255)
                                                                    cGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGc    >  1:940922/1‑76 (MQ=255)
                                                                      cAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCg  <  2:779839/76‑1 (MQ=255)
                                                                         |                                                                        
GACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCG  >  REL606/3251939‑3252084

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