Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,440,194 C→A G279V (GGT→GTT)  rpoA ← RNA polymerase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,440,1940CA83.5% 300.3 / 75.2 127G279V (GGT→GTT) rpoARNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (12/9);  new base A (53/53);  total (65/62)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

Only 100 of 128 total aligned reads displayed.

AGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACA  >  NC_000913/3440146‑3440241
                                                |                                               
aGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAcc                                                <  1:6654747/50‑1 (MQ=255)
  gTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGa                                              >  1:168785/1‑50 (MQ=255)
  gTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGa                                              >  1:5666486/1‑50 (MQ=255)
   ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGa                                              >  1:8125663/1‑49 (MQ=255)
   ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGa                                              >  1:5308328/1‑49 (MQ=255)
     aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                            <  1:107552/49‑1 (MQ=255)
     aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                            <  1:4211885/49‑1 (MQ=255)
     aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                            <  1:183170/49‑1 (MQ=255)
     aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                            >  1:5180286/1‑49 (MQ=255)
     aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                            <  1:2659789/49‑1 (MQ=255)
      ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGAtata                                          <  1:406700/50‑1 (MQ=255)
      ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGAtata                                          <  1:4484558/50‑1 (MQ=255)
       gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAg                                         >  1:2825614/1‑50 (MQ=255)
        cGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGt                                        <  1:8133986/50‑1 (MQ=255)
        cGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGt                                        >  1:3564662/1‑50 (MQ=255)
         gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTg                                       >  1:5906747/1‑50 (MQ=255)
         gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTg                                       >  1:1673752/1‑50 (MQ=255)
         gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                       >  1:6543073/1‑50 (MQ=255)
         gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                       >  1:8046154/1‑50 (MQ=255)
          ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTgg                                      <  1:3015682/50‑1 (MQ=255)
          ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTgg                                      <  1:30993/50‑1 (MQ=255)
          ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                       >  1:3448018/1‑49 (MQ=255)
            ttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                     <  1:6748264/49‑1 (MQ=255)
            ttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                     <  1:7222711/49‑1 (MQ=255)
            ttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                     <  1:3208637/49‑1 (MQ=255)
             tAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGa                                     >  1:7538948/1‑48 (MQ=255)
             tAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                     >  1:3869692/1‑48 (MQ=255)
             tAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                     >  1:5436270/1‑48 (MQ=255)
             tAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGAt                                    <  1:3855468/49‑1 (MQ=255)
              aaGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATa                                   <  1:5846893/49‑1 (MQ=255)
              aaGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAg                                  >  1:830696/1‑50 (MQ=255)
               aGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAg                                  <  1:8099758/49‑1 (MQ=255)
               aGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGc                                 >  1:6625238/1‑50 (MQ=255)
               aGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGc                                 >  1:268740/1‑50 (MQ=255)
                ggAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCt                                <  1:1988022/50‑1 (MQ=255)
                 gAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCtt                               <  1:3472214/50‑1 (MQ=255)
                 gAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCtt                               >  1:3613182/1‑50 (MQ=255)
                 gAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCtt                               <  1:7003044/50‑1 (MQ=255)
                  aGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTc                              <  1:3182134/50‑1 (MQ=255)
                  aGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTc                              >  1:4848115/1‑50 (MQ=255)
                  aGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTc                              <  1:827045/50‑1 (MQ=255)
                   gCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCt                             <  1:4471594/50‑1 (MQ=255)
                   gCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCt                             <  1:2830874/50‑1 (MQ=255)
                   gCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCt                             <  1:1445485/50‑1 (MQ=255)
                    cTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTg                            >  1:8028435/1‑50 (MQ=255)
                    cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg                            >  1:4186162/1‑50 (MQ=255)
                    cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg                            <  1:4669585/50‑1 (MQ=255)
                     tCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGc                           >  1:1435883/1‑50 (MQ=255)
                     tCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGc                           >  1:7451901/1‑50 (MQ=255)
                     tCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGc                           >  1:328220/1‑50 (MQ=255)
                     tCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGc                           >  1:5321762/1‑50 (MQ=255)
                      cAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                          >  1:8270831/1‑50 (MQ=255)
                      cAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                          >  1:5710018/1‑50 (MQ=255)
                       aaCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCCTCTGCtt                         >  1:709615/1‑50 (MQ=255)
                       aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                         <  1:1066373/50‑1 (MQ=255)
                       aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                         >  1:339160/1‑50 (MQ=255)
                       aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCt                          <  1:6937838/49‑1 (MQ=255)
                        aCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                         >  1:1907236/1‑49 (MQ=255)
                         ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt                        >  1:2514031/1‑49 (MQ=255)
                         ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTa                       <  1:5621899/50‑1 (MQ=255)
                          cTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTaa                      <  1:3814003/50‑1 (MQ=255)
                          cTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTaa                      <  1:5629008/50‑1 (MQ=255)
                          cTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTaa                      <  1:4333911/50‑1 (MQ=255)
                            cGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                    >  1:1971900/1‑50 (MQ=255)
                             ggTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg                    >  1:8084576/1‑49 (MQ=255)
                             ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                    <  1:391440/49‑1 (MQ=255)
                             ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAg                     <  1:4653407/48‑1 (MQ=255)
                             ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGc                   <  1:7169443/50‑1 (MQ=255)
                             ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGc                   <  1:2134623/50‑1 (MQ=255)
                              gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                    >  1:4388470/1‑48 (MQ=255)
                              gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                    >  1:405426/1‑48 (MQ=255)
                              gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                    >  1:276822/1‑48 (MQ=255)
                               tACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAg                 <  1:3810761/50‑1 (MQ=255)
                               tACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAg                 >  1:3174794/1‑50 (MQ=255)
                                aCGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                >  1:5985742/1‑50 (MQ=255)
                                  gCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                >  1:7357659/1‑48 (MQ=255)
                                  gCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                >  1:3562488/1‑48 (MQ=255)
                                   cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg             <  1:2654723/50‑1 (MQ=255)
                                   cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg             >  1:7471312/1‑50 (MQ=255)
                                    tGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg             >  1:4098082/1‑49 (MQ=255)
                                     gTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg             >  1:1119505/1‑48 (MQ=255)
                                     gTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg             >  1:5133116/1‑48 (MQ=255)
                                      tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg             >  1:8047154/1‑47 (MQ=255)
                                      tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc            <  1:1074324/48‑1 (MQ=255)
                                      tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc            <  1:6009029/48‑1 (MQ=255)
                                      tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc            <  1:4482485/48‑1 (MQ=255)
                                      tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCag          <  1:4115881/50‑1 (MQ=255)
                                      tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCag          <  1:5696712/50‑1 (MQ=255)
                                       aCCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa           >  1:4555318/1‑48 (MQ=255)
                                        ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag        >  1:7283730/1‑50 (MQ=255)
                                        ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag        <  1:5732188/50‑1 (MQ=255)
                                        ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCaga         <  1:4349771/49‑1 (MQ=255)
                                         cAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc       >  1:5245697/1‑50 (MQ=255)
                                         cAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc       >  1:5317303/1‑50 (MQ=255)
                                          aGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCg      <  1:250674/50‑1 (MQ=255)
                                           gATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg     <  1:7315384/50‑1 (MQ=255)
                                           gATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg     <  1:2075631/50‑1 (MQ=255)
                                            aTCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa    >  1:4784963/1‑50 (MQ=255)
                                            aTCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa    >  1:2146924/1‑50 (MQ=255)
                                               aCCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa  <  1:2385514/49‑1 (MQ=255)
                                                |                                               
AGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACA  >  NC_000913/3440146‑3440241

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