Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,145,984 A→T I41F (ATC→TTC)  rpoD → RNA polymerase sigma factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,145,9840AT96.3% 95.4 / NA 27I41F (ATC→TTC) rpoDRNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (14/12);  minor base G (0/1);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGGAGCAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATC  >  REL606/3145911‑3146056
                                                                         |                                                                        
tAAGGAGCAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCagatcagttc                                                                        <  1:1149492/76‑1 (MQ=255)
 aaGGAACAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCg                                                                       >  1:1139204/1‑76 (MQ=255)
     aGCAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGa                                                                   <  1:981656/76‑1 (MQ=255)
      gCAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGAc                                                                  <  1:660304/76‑1 (MQ=255)
       cAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGAca                                                                 <  1:1154238/76‑1 (MQ=255)
            cTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGACatcatc                                                            >  2:897824/1‑76 (MQ=255)
              aTCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTCCGAAGACATCATCCa                                                          <  1:110340/76‑1 (MQ=255)
                  gACCTATGCCGAGGTCAATGACCTTATGCCGGAAGATATNGTCGAATCAGATCAGNTCGAAGTCATCATCCaaacg                                                      >  2:1197837/1‑74 (MQ=255)
                  gACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAGGACATCATCCAAATg                                                      >  1:7881/1‑76 (MQ=255)
                  gACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATg                                                      >  1:589437/1‑76 (MQ=255)
                         gCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACg                                               <  1:574708/76‑1 (MQ=255)
                            gAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACa                                            >  1:613221/1‑76 (MQ=255)
                              ggTCAATGACCATCTGCGGTAAGATATCGTCGATTCAGATCAGTTCGAAGATATCATCCAAATGATCAACGACATg                                          <  1:55055/76‑1 (MQ=255)
                                    ttaCCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCAtt                                    <  2:1093761/74‑1 (MQ=255)
                                    tGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCAtt                                    >  2:60980/1‑76 (MQ=255)
                                     gACCATCTGCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTc                                   <  2:294036/76‑1 (MQ=255)
                                      aCCATCTGNCGGGAGATATCGTCGTTTGAGATCTGGTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCa                                  <  2:202583/76‑1 (MQ=255)
                                          tCTGCCGGAAGATATCGACGATGCAGATCAGTTCGAAGACATCATCCAAATGATCCACGACATGGGCATTCAGGTg                              <  1:55966/76‑1 (MQ=255)
                                             gCCGGAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATg                           >  2:639356/1‑76 (MQ=255)
                                              ccGGAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATgg                          >  1:782006/1‑76 (MQ=255)
                                               cGGAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGga                         >  2:1003205/1‑76 (MQ=255)
                                                ggAAGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGgaa                        >  2:1132719/1‑76 (MQ=255)
                                                  aaGATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGgaaga                      <  2:756585/76‑1 (MQ=255)
                                                    gATATCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGTTGATGGaagaag                    <  2:48148/76‑1 (MQ=255)
                                                     atatCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATNGAAGAAGc                   <  1:1167822/76‑1 (MQ=255)
                                                      tatCGTCGATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCa                  >  1:976058/1‑76 (MQ=255)
                                                         cgtcgATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCg               >  1:575462/1‑76 (MQ=255)
                                                         cgtcgATTCAGAGCACTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGNTGGAAGAAGCACCg               >  1:1038328/1‑76 (MQ=255)
                                                            cgATTCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGAt            >  2:165716/1‑76 (MQ=255)
                                                               ttCAGATCAGTTCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGcc         >  2:379855/1‑76 (MQ=255)
                                                                    atcagttcGAAGACATCATCCAAATGATCAACGACATNGGCATTCAGGTGATGGAAGAAGCACCGGATGCCgatga    >  1:1185766/1‑76 (MQ=255)
                                                                    atcagttcGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCgatga    <  2:1112352/76‑1 (MQ=255)
                                                                     tcagttcGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACAGGATGCCgatgat   >  2:555448/1‑76 (MQ=255)
                                                                      cagttcGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATc  >  1:735573/1‑76 (MQ=255)
                                                                         |                                                                        
TAAGGAGCAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATC  >  REL606/3145911‑3146056

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