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,8130CT98.3% 226.4 / ‑3.4 60Q260* (CAG→TAG) hnrresponse regulator of RpoS
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/28);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

CAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATG  >  REL606/1290741‑1290885
                                                                        |                                                                        
cAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGt                                                                       >  2:1863301/1‑76 (MQ=255)
cAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGt                                                                       <  2:547493/76‑1 (MQ=255)
cAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGt                                                                       <  2:1570367/76‑1 (MQ=255)
 aaCGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGtt                                                                      <  1:1674293/76‑1 (MQ=255)
  aCGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGttt                                                                     >  1:897930/1‑76 (MQ=255)
   cGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTc                                                                    >  1:1154067/1‑76 (MQ=255)
          cAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTAtt                                                             <  1:457610/76‑1 (MQ=255)
            gagTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAg                                                           <  2:755269/76‑1 (MQ=255)
               ttGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTg                                                        <  1:580379/76‑1 (MQ=255)
                   gcgcGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCta                                                    <  2:819493/76‑1 (MQ=255)
                   gcgcGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCta                                                    <  1:1580050/76‑1 (MQ=255)
                     gcgTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCtatt                                                  <  1:522619/76‑1 (MQ=255)
                      cgTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCtatta                                                 >  1:1779146/1‑76 (MQ=255)
                      cgTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCtatta                                                 >  1:1331340/1‑76 (MQ=255)
                      cgTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCtatta                                                 >  2:839873/1‑76 (MQ=255)
                      cgTTAATGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCtatta                                                 <  2:1974779/76‑1 (MQ=255)
                       gTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCtattat                                                >  2:419119/1‑76 (MQ=255)
                        ttATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTatc                                               >  2:1793671/1‑76 (MQ=255)
                          aTTGAAGCAGCTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTatcat                                             <  2:1388406/76‑1 (MQ=255)
                            tGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATcg                                           <  1:1428880/76‑1 (MQ=255)
                            tGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATcg                                           >  2:658936/1‑76 (MQ=255)
                            tGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATcg                                           <  2:1544885/76‑1 (MQ=255)
                              aaGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATcgcg                                         <  2:1186810/76‑1 (MQ=255)
                              aaGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATcgcg                                         >  1:940663/1‑76 (MQ=255)
                                gCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGaa                                       >  1:1488760/1‑76 (MQ=255)
                                gCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGaa                                       <  1:1279685/76‑1 (MQ=255)
                                 cAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAAc                                      <  2:1194488/76‑1 (MQ=255)
                                  aGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACt                                     >  1:1761214/1‑76 (MQ=255)
                                   ggTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTg                                    <  1:221222/76‑1 (MQ=255)
                                      aaaCCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGaaa                                 <  2:1507421/76‑1 (MQ=255)
                                      aaaCCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGaaa                                 <  1:1167654/76‑1 (MQ=255)
                                       aaCCATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGaaaa                                >  2:276692/1‑76 (MQ=255)
                                          cATTTACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAAtc                             >  2:1478603/1‑76 (MQ=255)
                                            tttACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAAtctc                           <  1:1693068/76‑1 (MQ=255)
                                             ttACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCa                          >  2:481329/1‑76 (MQ=255)
                                             ttACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCa                          >  2:130767/1‑76 (MQ=255)
                                             ttACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCa                          >  2:669470/1‑76 (MQ=255)
                                              tACTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCAt                         >  2:1245197/1‑76 (MQ=255)
                                               aCTCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCAtt                        >  1:485069/1‑76 (MQ=255)
                                                 tCCGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCt                      <  1:324540/76‑1 (MQ=255)
                                                   cGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTgg                    <  1:297270/76‑1 (MQ=255)
                                                   cGTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTgg                    >  1:1177876/1‑76 (MQ=255)
                                                    gTCAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGt                   >  2:961001/1‑76 (MQ=255)
                                                      cAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGttt                 >  2:1825022/1‑76 (MQ=255)
                                                      cAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGttt                 >  2:955876/1‑76 (MQ=255)
                                                      cAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGttt                 <  1:301486/76‑1 (MQ=255)
                                                      cAGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGttt                 >  2:538336/1‑76 (MQ=255)
                                                       aGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTc                >  2:232468/1‑76 (MQ=255)
                                                       aGGCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTc                <  2:1075251/76‑1 (MQ=255)
                                                        ggCCAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCt               >  1:1481497/1‑76 (MQ=255)
                                                           naaTCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCg            >  1:1451172/2‑76 (MQ=255)
                                                           cAATCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCg            >  1:1262/1‑76 (MQ=255)
                                                             aTCTGCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCggg          <  2:220668/76‑1 (MQ=255)
                                                                 gCCGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTg      >  2:39673/1‑76 (MQ=255)
                                                                  ccGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGa     <  2:1171875/76‑1 (MQ=255)
                                                                   cGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGaa    >  1:548355/1‑76 (MQ=255)
                                                                   cGGGGTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGaa    >  1:825693/1‑76 (MQ=255)
                                                                    ggggTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAAt   <  1:290713/76‑1 (MQ=255)
                                                                     gggTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATg  <  1:288655/76‑1 (MQ=255)
                                                                     gggTAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATg  <  1:1745393/76‑1 (MQ=255)
                                                                        |                                                                        
CAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTCCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATG  >  REL606/1290741‑1290885

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