Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,290,813 C→T Q260* (CAG→TAG)  hnr → response regulator of RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,290,8130CT97.4% 133.7 / ‑3.1 38Q260* (CAG→TAG) hnrresponse regulator of RpoS
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/18);  total (20/18)
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.55e-01

ACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAA  >  REL606/1290743‑1290883
                                                                      |                                                                      
aCGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGttt                                                                   <  2:426386/76‑1 (MQ=255)
   gTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCg                                                                >  2:1234565/1‑76 (MQ=255)
   gTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCg                                                                <  2:169988/76‑1 (MQ=255)
    ttGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGc                                                               >  1:77002/1‑76 (MQ=255)
    ttGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGc                                                               >  1:1458822/1‑76 (MQ=255)
        cAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTAtt                                                           >  1:15064/1‑76 (MQ=255)
        cAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTAtt                                                           <  2:1526976/76‑1 (MQ=255)
        cAGAGTGGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTAtt                                                           <  2:928435/76‑1 (MQ=255)
          gagTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAg                                                         <  1:511761/76‑1 (MQ=255)
             ttGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTg                                                      >  1:1569741/1‑76 (MQ=255)
               gggCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGc                                                    >  2:1391534/1‑76 (MQ=255)
                ggCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGca                                                   <  1:154399/76‑2 (MQ=255)
                ggCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCt                                                   <  1:1577675/76‑1 (MQ=255)
                ggCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCt                                                   <  2:738285/76‑1 (MQ=255)
                      ttATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTatc                                             >  2:1430406/1‑76 (MQ=255)
                         ttGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTatcatc                                          <  2:1592333/76‑1 (MQ=255)
                          tGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATcg                                         <  2:1176225/76‑1 (MQ=255)
                            aaGCAGGAAAACCATTTACCCCGTCAGTCAAACCGCCCGGGCTATTTTCCGCTCTTAGTTGGCTATTATCATcgcg                                       >  2:1048998/1‑76 (MQ=255)
                                    aaaCCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGaaa                               >  1:1218451/1‑76 (MQ=255)
                                    aaaCCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGaaa                               >  2:838320/1‑76 (MQ=255)
                                      aCCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGaaaaa                             <  2:1636552/76‑1 (MQ=255)
                                      aCCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGaaaaa                             <  2:153157/76‑1 (MQ=255)
                                        cATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAAtc                           >  2:1330932/1‑76 (MQ=255)
                                           ttACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCa                        >  1:821679/1‑76 (MQ=255)
                                           ttACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCa                        >  1:15228/1‑76 (MQ=255)
                                            tACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCAt                       >  1:392221/1‑76 (MQ=255)
                                             aCTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCAtt                      <  2:497804/76‑1 (MQ=255)
                                               tCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCt                    <  1:1653674/76‑1 (MQ=255)
                                               tCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCt                    >  1:233285/1‑76 (MQ=255)
                                                ccGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTg                   >  1:746762/1‑76 (MQ=255)
                                                ccGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTg                   >  1:1114995/1‑76 (MQ=255)
                                                  gTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTGGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTggn                 <  1:1279494/76‑2 (MQ=255)
                                                  gTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGt                 <  2:464006/76‑1 (MQ=255)
                                                   tCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGtt                >  1:1248383/1‑76 (MQ=255)
                                                        ccAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGc           <  1:448452/76‑1 (MQ=255)
                                                             cTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTc      >  1:1227136/1‑76 (MQ=255)
                                                               gCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTg    <  2:1381145/76‑1 (MQ=255)
                                                                 cGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGaa  >  1:917272/1‑76 (MQ=255)
                                                                 cGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGaa  >  2:440204/1‑76 (MQ=255)
                                                                      |                                                                      
ACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAA  >  REL606/1290743‑1290883

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