Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,290,440 Δ1 bp coding (405/1014 nt) hnr → response regulator of RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,290,4390A.97.4% 182.7 / ‑3.4 39coding (404/1014 nt)hnrresponse regulator of RpoS
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (17/21);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGG  >  REL606/1290365‑1290510
                                                                          |                                                                       
tAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTgaggag                                                                        <  1:552570/76‑2 (MQ=255)
 aaaGATCTGAATCGCTTGCGCGAGATGGTTTTTGCTTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTgaggaga                                                                       >  2:1662466/1‑74 (MQ=255)
 aaaGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTgaggaga                                                                       >  1:1139056/1‑74 (MQ=255)
 aaaGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTgaggaga                                                                       <  2:1104756/76‑3 (MQ=255)
  aaGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTgaggagag                                                                      >  1:1349834/1‑73 (MQ=255)
    gATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTtgaggagagga                                                                   >  2:1142729/1‑76 (MQ=255)
    gATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTtgaggagagga                                                                   >  1:1109561/1‑76 (MQ=255)
     aTCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTtgaggagaggaa                                                                  >  2:240022/1‑76 (MQ=255)
     aTCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTtgaggagaggaa                                                                  <  1:361851/76‑1 (MQ=255)
            tCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCttt                                                           <  2:1019594/76‑1 (MQ=255)
             cGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGATAGGAAAGGCtttt                                                          <  1:1211872/76‑1 (MQ=255)
             cGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCtttt                                                          <  2:221285/76‑1 (MQ=255)
                    gcgAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGAc                                                   >  1:238148/1‑76 (MQ=255)
                    gcgAGATGGTTTTNGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGAc                                                   >  1:1242124/1‑76 (MQ=255)
                     cgAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACt                                                  <  1:726407/76‑1 (MQ=255)
                      gagaTGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTg                                                 >  1:55682/1‑76 (MQ=255)
                      gagaTGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTg                                                 >  2:1181683/1‑76 (MQ=255)
                         aTGGTTTTTGCCTGNCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGa                                              >  1:993241/1‑76 (MQ=255)
                          tGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTCTCTCGCGACTGGGAt                                             >  2:1739750/1‑76 (MQ=255)
                            gTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGc                                           >  1:526632/1‑76 (MQ=255)
                               tttGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAAt                                        <  2:313615/76‑1 (MQ=255)
                                ttGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCTACTGGGATGCAATg                                       <  1:978514/76‑1 (MQ=255)
                                  gCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGt                                     >  2:1115079/1‑76 (MQ=255)
                                  gCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATACAATGGt                                     <  1:1191351/76‑1 (MQ=255)
                                   ccTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGACGCAATGGtt                                    >  1:665399/1‑76 (MQ=255)
                                      gTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGAt                                 >  2:1663922/1‑76 (MQ=255)
                                       tctcTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATa                                <  2:174330/76‑1 (MQ=255)
                                       tctcTATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATa                                >  1:521949/1‑76 (MQ=255)
                                         tctATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAAc                              <  1:959967/76‑1 (MQ=255)
                                         tctATCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAAc                              >  1:1674839/1‑76 (MQ=255)
                                            aTCCCAGCTGGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCt                           <  2:888345/76‑1 (MQ=255)
                                            aTCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCt                           >  1:1752026/1‑76 (MQ=255)
                                            aTCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCt                           <  2:221406/76‑1 (MQ=255)
                                             tCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTg                          <  1:938328/76‑1 (MQ=255)
                                             tCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTg                          <  1:953709/76‑1 (MQ=255)
                                             tCCCAGCATGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTg                          <  1:953693/76‑1 (MQ=255)
                                                    aTGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAgc                   >  1:1085259/1‑76 (MQ=255)
                                                    aTGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAgc                   >  2:459897/1‑76 (MQ=255)
                                                     tGTTTAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAgcg                  <  1:1724840/76‑1 (MQ=255)
                                                       tttAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAgcggc                <  1:1405500/76‑1 (MQ=255)
                                                         tAATTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGa              <  2:7792/76‑1 (MQ=255)
                                                          aaTTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGaa             >  1:820293/1‑76 (MQ=255)
                                                          aaTTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGaa             <  1:758478/76‑1 (MQ=255)
                                                           aTTCGCGCGTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGaaa            <  2:1111170/76‑1 (MQ=255)
                                                              ggggcgTTGAGGGGAGGAAGGGGTTTTGGGCGGCTGGGATGCAATGGTTGATAACCCTGCCGCGGCGGCGAAatta         <  2:1249923/73‑1 (MQ=255)
                                                              cgcgcgTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAatta         <  1:1649276/76‑1 (MQ=255)
                                                                 gcgTTGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAAttatta      >  1:1781218/1‑76 (MQ=255)
                                                                    ttGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAg   >  2:163504/1‑76 (MQ=255)
                                                                     tGAGGAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAgg  >  1:1089058/1‑76 (MQ=255)
                                                                          |                                                                       
TAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGG  >  REL606/1290365‑1290510

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