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,8030TC94.0% 607.2 / 15.9 168S253P (TCT→CCT) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base T (2/8);  new base C (80/78);  total (82/86)
Fisher's exact test for biased strand distribution p-value = 9.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 168 total aligned reads displayed.

CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTA  >  NC_000913/3213756‑3213849
                                               |                                              
cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              <  1:4281265/50‑1 (MQ=255)
 gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              >  1:484037/1‑49 (MQ=255)
  aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                             >  2:4951406/1‑49 (MQ=255)
  aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGa                                            <  2:5568086/50‑1 (MQ=255)
   aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGaa                                           <  1:3529524/50‑1 (MQ=255)
    aGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                          <  1:1156426/50‑1 (MQ=255)
     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGt                                         <  2:4568074/50‑1 (MQ=255)
     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGt                                         <  1:754691/50‑1 (MQ=255)
     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGt                                         <  1:99295/50‑1 (MQ=255)
      gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                        <  1:1015138/50‑1 (MQ=255)
      gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                        <  1:5627251/50‑1 (MQ=255)
      gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                        >  1:813555/1‑50 (MQ=255)
       tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt                                       <  1:4420894/50‑1 (MQ=255)
       tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt                                       >  2:4974126/1‑50 (MQ=255)
        cGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  1:3634930/1‑50 (MQ=255)
        cGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  2:4588278/1‑50 (MQ=255)
         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAttg                                     <  1:5469653/50‑2 (MQ=255)
         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      <  2:699862/49‑1 (MQ=255)
         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      <  1:4648843/49‑1 (MQ=255)
         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  2:4761452/1‑49 (MQ=255)
         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                     <  1:3555911/50‑1 (MQ=255)
         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                     <  2:3227156/50‑1 (MQ=255)
         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                     <  2:1920343/50‑1 (MQ=255)
          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCa                                    <  1:605908/50‑1 (MQ=255)
          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                     >  1:594138/1‑49 (MQ=255)
          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    <  1:1784868/50‑1 (MQ=255)
          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    <  1:3360224/50‑1 (MQ=255)
           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaa                                   <  2:4222778/50‑1 (MQ=255)
           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    >  1:3896818/1‑49 (MQ=255)
           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    >  2:901963/1‑49 (MQ=255)
           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    >  2:3106116/1‑49 (MQ=255)
           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    >  1:345818/1‑49 (MQ=255)
            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  >  1:1871455/1‑50 (MQ=255)
             tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                 <  1:1338794/50‑1 (MQ=255)
             tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                 <  1:4143183/50‑1 (MQ=255)
             tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                 <  1:5521906/50‑1 (MQ=255)
               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACa                                <  2:2548945/49‑1 (MQ=255)
               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACa                                <  1:845286/49‑1 (MQ=255)
               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                               >  1:5354237/1‑50 (MQ=255)
               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                               <  2:57151/50‑1 (MQ=255)
                 gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                             <  2:3967165/50‑1 (MQ=255)
                 gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                             <  1:2839814/50‑1 (MQ=255)
                 gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                             >  1:141129/1‑50 (MQ=255)
                  cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                            <  2:3737080/50‑1 (MQ=255)
                   tACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                           >  2:2208865/1‑50 (MQ=255)
                   tACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                           >  1:4023005/1‑50 (MQ=255)
                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                           <  1:2892351/48‑1 (MQ=255)
                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                          <  1:1448681/49‑1 (MQ=255)
                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                          <  1:1404850/49‑1 (MQ=255)
                      cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                        >  2:2845964/1‑50 (MQ=255)
                      cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                        <  2:4915246/50‑1 (MQ=255)
                      cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                        >  2:3572646/1‑50 (MQ=255)
                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                       <  1:5478945/50‑1 (MQ=255)
                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                       <  1:5307361/50‑1 (MQ=255)
                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTACCGCCt                       <  2:3833061/50‑1 (MQ=39)
                        cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTg                      <  2:4503522/50‑1 (MQ=255)
                        cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTg                      >  2:713112/1‑50 (MQ=255)
                         tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                     >  2:5386601/1‑50 (MQ=255)
                         tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                     >  2:1729178/1‑50 (MQ=255)
                          cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                    >  1:5706528/1‑50 (MQ=255)
                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                    <  2:5266368/49‑1 (MQ=255)
                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                    <  1:3322550/49‑1 (MQ=255)
                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   <  1:5046238/50‑1 (MQ=255)
                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   >  1:3482419/1‑50 (MQ=255)
                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   >  2:2700992/1‑50 (MQ=255)
                            ggAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGc                  <  2:5470482/50‑1 (MQ=255)
                            ggAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGc                  >  1:4471505/1‑50 (MQ=255)
                             gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                 >  1:2790230/1‑50 (MQ=255)
                             gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                 >  2:3134411/1‑50 (MQ=255)
                              aaGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg                <  2:4258689/50‑1 (MQ=255)
                              aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                 >  1:2883201/1‑49 (MQ=255)
                                 agaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  1:1542568/1‑50 (MQ=255)
                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  1:2523032/1‑49 (MQ=255)
                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  1:2783253/1‑49 (MQ=255)
                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  1:2246525/1‑49 (MQ=255)
                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  2:3403893/1‑48 (MQ=255)
                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  2:1538352/1‑48 (MQ=255)
                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCa           <  2:1216742/50‑1 (MQ=255)
                                    tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg          <  2:5264716/50‑1 (MQ=255)
                                     ccTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt         >  2:3019932/1‑50 (MQ=255)
                                      cTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt        >  2:195758/1‑50 (MQ=255)
                                      cTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt         <  1:3298825/49‑1 (MQ=255)
                                      cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt        >  1:2536547/1‑50 (MQ=255)
                                      cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt        <  1:1188291/50‑1 (MQ=255)
                                      cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt        >  1:2494303/1‑50 (MQ=255)
                                      cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt        >  2:1480885/1‑50 (MQ=255)
                                       tGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt       >  1:4109388/1‑50 (MQ=255)
                                         aaaCTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTg      <  1:1504848/49‑1 (MQ=255)
                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCTAAGCAGTTTGa     >  1:1100520/1‑50 (MQ=39)
                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt       <  1:2988718/48‑1 (MQ=255)
                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTg      <  1:1867631/49‑1 (MQ=255)
                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa     >  1:1016655/1‑50 (MQ=255)
                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa     >  2:3467286/1‑50 (MQ=255)
                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa     >  2:1901890/1‑50 (MQ=255)
                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa     >  1:4724884/1‑50 (MQ=255)
                                          aaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc    >  2:16686/1‑50 (MQ=255)
                                           aCTGCCTGAAGTATTCAATCAGTTCCGTCTGGTGCCGAAGCAGTTTGTCt   <  2:3860482/50‑1 (MQ=37)
                                           aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc    >  2:2262458/1‑49 (MQ=255)
                                           aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACt   <  2:805056/50‑1 (MQ=255)
                                            cTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTa  >  2:3046832/1‑50 (MQ=255)
                                               |                                              
CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTA  >  NC_000913/3213756‑3213849

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