Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,213,803 T→C S253P (TCT→CCT)  rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,213,8030TC78.7% 309.7 / 85.4 136S253P (TCT→CCT) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base T (17/12);  new base C (61/46);  total (78/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 138 total aligned reads displayed.

GCGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCT  >  NC_000913/3213755‑3213852
                                                |                                                 
gCGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTc                                                  <  1:4513783/50‑1 (MQ=255)
 cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCt                                                 <  2:4880597/50‑1 (MQ=255)
 cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                                 >  2:4394454/1‑50 (MQ=255)
  gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                                 >  2:3814508/1‑49 (MQ=255)
  gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                                 >  1:4255975/1‑49 (MQ=255)
  gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                                 >  1:2663190/1‑49 (MQ=255)
  gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                                 >  1:44207/1‑49 (MQ=255)
   aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                                >  2:2245538/1‑49 (MQ=255)
    aaGGTCGCAGTCACGCTACCGCTCAGGAAGTGATCCTGAAACTGCCTg                                                >  1:4406262/1‑48 (MQ=255)
    aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTg                                                >  2:536948/1‑48 (MQ=255)
    aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTg                                                >  2:5102247/1‑48 (MQ=255)
    aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                                >  2:4097214/1‑48 (MQ=255)
       gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                           >  1:540753/1‑50 (MQ=255)
        tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt                                          >  1:3557651/1‑50 (MQ=255)
        tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt                                          >  2:3723925/1‑50 (MQ=255)
         cGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                         >  1:973158/1‑50 (MQ=255)
          gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTAtt                                         <  1:214731/49‑1 (MQ=255)
          gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTc                                        <  1:3886739/50‑1 (MQ=255)
          gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                         >  2:2545908/1‑49 (MQ=255)
          gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                        <  1:13238/50‑1 (MQ=255)
           cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCa                                       <  1:54012/50‑1 (MQ=255)
           cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCa                                       <  2:5091436/50‑1 (MQ=255)
           cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                       <  2:4390661/50‑1 (MQ=255)
            aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaa                                      <  1:3860037/50‑1 (MQ=255)
            aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                       >  2:3815567/1‑49 (MQ=255)
             gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCaaa                                     <  1:2802762/50‑1 (MQ=255)
             gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                     >  1:1828023/1‑50 (MQ=255)
              tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                    <  1:4888755/50‑1 (MQ=255)
               cACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                     >  1:3653445/1‑48 (MQ=255)
               cACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                     >  2:4907657/1‑48 (MQ=255)
                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAg                                  <  2:4073658/50‑1 (MQ=255)
                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACa                                   <  1:5298793/49‑1 (MQ=255)
                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACa                                   <  1:2889123/49‑1 (MQ=255)
                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACa                                   <  1:1379826/49‑1 (MQ=255)
                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACa                                   <  1:6980/49‑1 (MQ=255)
                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                                  <  1:1025606/50‑1 (MQ=255)
                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                                  <  1:3151383/50‑1 (MQ=255)
                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                                  >  1:4960409/1‑50 (MQ=255)
                 cGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGt                                 >  2:1997604/1‑50 (MQ=255)
                 cGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGt                                 <  2:2474957/50‑1 (MQ=255)
                 cGCTACCGCTCAGGAAGAAATCCTGAAACTGCCTGAAGTATTCAAACAGt                                 >  2:3708590/1‑50 (MQ=39)
                  gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                                >  1:4193790/1‑50 (MQ=255)
                  gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                                <  2:1037963/50‑1 (MQ=255)
                  gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                                >  2:1153343/1‑50 (MQ=255)
                  gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                                <  1:3610402/50‑1 (MQ=255)
                   cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                               >  2:2709960/1‑50 (MQ=255)
                   cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                               <  1:193046/50‑1 (MQ=255)
                   cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                               <  2:3572938/50‑1 (MQ=255)
                   cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                               <  1:3071744/50‑1 (MQ=255)
                     aCCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                             >  1:1796515/1‑50 (MQ=255)
                      ccGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGc                            <  1:2000982/50‑1 (MQ=255)
                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                              <  1:4288560/48‑1 (MQ=255)
                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                              <  1:1563647/48‑1 (MQ=255)
                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                              <  2:2256719/48‑1 (MQ=255)
                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                             <  1:2797451/49‑1 (MQ=255)
                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGc                            <  1:2004305/50‑1 (MQ=255)
                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGc                            >  2:1813604/1‑50 (MQ=255)
                          tCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTgg                        >  1:3199575/1‑50 (MQ=255)
                          tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                        >  1:4784471/1‑50 (MQ=255)
                          tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                        >  2:1824008/1‑50 (MQ=255)
                          tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                        >  1:3817291/1‑50 (MQ=255)
                          tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                        >  1:3436620/1‑50 (MQ=255)
                           cAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGt                       >  2:3940157/1‑50 (MQ=255)
                           cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                       >  2:1771441/1‑50 (MQ=255)
                            aGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGt                       <  2:603190/49‑1 (MQ=255)
                            aGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGt                       <  1:4479782/49‑1 (MQ=255)
                            aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                       <  2:2536039/49‑1 (MQ=255)
                            aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                      <  1:4210002/50‑1 (MQ=255)
                            aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                      <  2:2205967/50‑1 (MQ=255)
                            aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                      >  2:2446741/1‑50 (MQ=255)
                             ggAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGc                     >  1:3741022/1‑50 (MQ=255)
                              gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                    >  2:1738145/1‑50 (MQ=255)
                               aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg                   <  1:1098146/50‑1 (MQ=255)
                               aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg                   <  1:5090002/50‑1 (MQ=255)
                                agagATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGa                  <  2:349693/50‑1 (MQ=255)
                                 gagaTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg                   >  1:1010234/1‑48 (MQ=255)
                                 gagaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg                   >  1:5339751/1‑48 (MQ=255)
                                 gagaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGaa                 <  1:1257510/50‑1 (MQ=255)
                                  agaTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg                >  1:4686889/1‑50 (MQ=255)
                                  agaTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg                >  1:319046/1‑50 (MQ=255)
                                  agaTCCTGAAACTGCCTGAAGTCTTCAAACAGTTCCGCCTGGTGCCGAAg                >  2:3510834/1‑50 (MQ=39)
                                   gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGc               <  2:570080/50‑1 (MQ=255)
                                   gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGc               <  2:4557593/50‑1 (MQ=255)
                                    aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg                >  1:3982975/1‑48 (MQ=255)
                                     tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg             >  1:4074530/1‑50 (MQ=255)
                                       cTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt            <  1:2053740/49‑1 (MQ=255)
                                       cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt            <  1:396327/49‑1 (MQ=255)
                                        tGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt          >  1:315212/1‑50 (MQ=255)
                                        tGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt          >  1:1016550/1‑50 (MQ=255)
                                         gAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt          >  1:2655952/1‑49 (MQ=255)
                                         gAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt          >  1:192012/1‑49 (MQ=255)
                                         gAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt          >  1:1107746/1‑49 (MQ=255)
                                          aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt          <  1:1990147/48‑1 (MQ=255)
                                          aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt          <  2:1271937/48‑1 (MQ=255)
                                          aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa        >  2:2720203/1‑50 (MQ=255)
                                          aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa        >  1:5015459/1‑50 (MQ=255)
                                          aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa        >  1:911828/1‑50 (MQ=255)
                                           aaCTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc       >  2:2546084/1‑50 (MQ=255)
                                           aaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc       <  1:1544531/50‑1 (MQ=255)
                                             cTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTa     <  2:3746046/50‑1 (MQ=255)
                                                |                                                 
GCGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCT  >  NC_000913/3213755‑3213852

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