Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,151 A→G F134S (TTT→TCT)  rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,1510AG98.3% 211.3 / ‑4.3 60F134S (TTT→TCT) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (29/30);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCA  >  NC_000913/2867059‑2867250
                                                                                            |                                                                                                   
tAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTtctc                                                                                              >  1:380005/1‑100 (MQ=255)
  gTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTa                                                                                            <  2:1378467/100‑1 (MQ=255)
   tACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTAc                                                                                           <  1:162986/100‑1 (MQ=255)
    aCGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTAcc                                                                                          >  2:455561/1‑100 (MQ=255)
         tttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACcgcgcg                                                                                     <  2:710462/100‑1 (MQ=255)
             gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATc                                                                                 >  1:1289903/1‑100 (MQ=255)
             gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATc                                                                                 >  2:1466988/1‑100 (MQ=255)
             gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATc                                                                                 <  2:49802/100‑1 (MQ=255)
              ttCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCa                                                                                >  1:633078/1‑100 (MQ=255)
                  tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGccc                                                                            >  1:470492/1‑100 (MQ=255)
                      cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTATCGCGCGGATCAGCCCCAgg                                                                        >  2:572356/1‑100 (MQ=255)
                      cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAgg                                                                        >  1:1293192/1‑100 (MQ=255)
                             tCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCct                                                                 <  1:525694/100‑1 (MQ=255)
                              cAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCctc                                                                >  2:452336/1‑100 (MQ=255)
                               aaTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCctct                                                               >  2:586473/1‑100 (MQ=255)
                                 tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTc                                                             <  1:193511/100‑1 (MQ=255)
                                 tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTc                                                             <  2:957324/100‑1 (MQ=255)
                                  cgccTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                            >  1:1143929/4‑100 (MQ=255)
                                  cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                            <  2:889548/100‑1 (MQ=255)
                                  cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                            >  1:64951/1‑100 (MQ=255)
                                  cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                            >  1:264811/1‑100 (MQ=255)
                                   gTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                           <  1:1093377/100‑1 (MQ=255)
                                   gTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                           <  1:576208/100‑1 (MQ=255)
                                    tCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGAt                                                          <  1:899591/100‑1 (MQ=255)
                                     cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATa                                                         <  2:69210/100‑1 (MQ=255)
                                     cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATa                                                         <  2:884171/100‑1 (MQ=255)
                                     cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATa                                                         >  1:817603/1‑100 (MQ=255)
                                        gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAACCCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                      <  1:1334393/100‑1 (MQ=255)
                                        gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                      <  2:533178/100‑1 (MQ=255)
                                        gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                      >  2:84982/1‑100 (MQ=255)
                                        gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                      <  1:640486/100‑1 (MQ=255)
                                          gAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTc                                                    >  1:1440928/1‑100 (MQ=255)
                                              ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                <  1:905238/100‑1 (MQ=255)
                                               caccaGGTTGCGTATGTTGAGAATCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                               >  1:1294415/1‑100 (MQ=255)
                                                 ccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAAc                                             <  1:977860/100‑1 (MQ=255)
                                                          gTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCa                                    >  1:44205/1‑100 (MQ=255)
                                                           tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCAc                                   >  1:81024/1‑100 (MQ=255)
                                                            aTGTTGAGAAGCGGAATCCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACg                                  <  1:886471/100‑1 (MQ=255)
                                                            aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACg                                  >  2:623348/1‑100 (MQ=255)
                                                            aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACg                                  <  1:663095/100‑1 (MQ=255)
                                                             tGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                 >  2:1307410/1‑100 (MQ=255)
                                                              nttGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAt                                >  1:742044/2‑100 (MQ=255)
                                                                tGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                              >  2:117561/1‑100 (MQ=255)
                                                                tGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                              >  1:163809/1‑100 (MQ=255)
                                                                 gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                             >  2:1205833/1‑100 (MQ=255)
                                                                  agaAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGcc                            <  2:792125/100‑1 (MQ=255)
                                                                  agaAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGcc                            >  2:1452315/1‑100 (MQ=255)
                                                                    aaGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCAt                          <  1:432174/100‑1 (MQ=255)
                                                                      gCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATaa                        <  2:264810/100‑1 (MQ=255)
                                                                      gCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATaa                        >  1:811561/1‑100 (MQ=255)
                                                                        ggAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                      <  1:357195/100‑1 (MQ=255)
                                                                              cACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                <  2:440932/100‑1 (MQ=255)
                                                                              cACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                >  2:341229/1‑100 (MQ=255)
                                                                              cACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                <  1:319111/100‑1 (MQ=255)
                                                                                    ccGGGGCAGACTTATCTACCGCGCGGATCAGCCTCAGGGTGACCTGTTCGATGAGGTCCAGCAGAGCCAGAACACGATTGCTATCACGGGCGACAAtt            >  2:1193118/1‑98 (MQ=255)
                                                                                       ggTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac       <  2:550728/100‑1 (MQ=255)
                                                                                       ggTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac       >  1:33448/1‑100 (MQ=255)
                                                                                        gTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc      <  1:1332544/100‑1 (MQ=255)
                                                                                        gTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc      <  2:735967/100‑1 (MQ=255)
                                                                                        gTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc      <  1:341230/100‑1 (MQ=255)
                                                                                            aaCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAccacca  >  2:1234144/1‑100 (MQ=255)
                                                                                            |                                                                                                   
TAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCA  >  NC_000913/2867059‑2867250

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