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,1940CA95.7% 247.2 / 5.1 70G279V (GGT→GTT) rpoARNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (37/30);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGA  >  NC_000913/3440147‑3440239
                                               |                                             
ggTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAcc                                              >  1:904501/1‑49 (MQ=255)
  ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGat                                           <  1:4681622/50‑1 (MQ=255)
  ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGat                                           <  1:5112676/50‑1 (MQ=255)
  ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGa                                            >  1:4818165/1‑49 (MQ=255)
   tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata                                          >  1:4349613/1‑50 (MQ=255)
   tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                          >  1:3814132/1‑50 (MQ=255)
   tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                          >  1:3931847/1‑50 (MQ=255)
   tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                          >  1:3707084/1‑50 (MQ=255)
    aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata                                          <  1:2037689/49‑1 (MQ=255)
    aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                          <  1:615750/49‑1 (MQ=255)
    aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                          <  1:1740767/49‑1 (MQ=255)
    aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                          >  1:1595373/1‑49 (MQ=255)
     ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGatat                                         <  1:4941528/49‑1 (MQ=255)
     ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                          >  1:3229384/1‑48 (MQ=255)
      gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAg                                       <  1:2934546/50‑1 (MQ=255)
      gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAg                                       <  1:2043796/50‑1 (MQ=255)
        gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                     >  1:1870220/1‑50 (MQ=255)
        gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                     <  1:4247227/50‑1 (MQ=255)
         ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTgg                                    <  1:2282071/50‑1 (MQ=255)
         ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTgg                                    <  1:4933949/50‑1 (MQ=255)
         ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                     >  1:4167682/1‑49 (MQ=255)
          tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                   >  1:3947216/1‑50 (MQ=255)
          tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                   >  1:4948198/1‑50 (MQ=255)
          tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                   >  1:3845673/1‑50 (MQ=255)
             aaGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATa                                 <  1:405531/49‑1 (MQ=255)
             aaGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAg                                >  1:190015/1‑50 (MQ=255)
              aGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGc                               >  1:4436282/1‑50 (MQ=255)
              aGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGc                               >  1:2126986/1‑50 (MQ=255)
              aGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGc                               <  1:3320206/50‑1 (MQ=255)
               ggAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCt                              >  1:2592147/1‑50 (MQ=255)
               ggAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCt                              >  1:1190483/1‑50 (MQ=255)
                 aGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTc                            <  1:3094205/50‑1 (MQ=255)
                  gCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCt                           <  1:2807718/50‑1 (MQ=255)
                   cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg                          >  1:1774632/1‑50 (MQ=255)
                   cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg                          >  1:3561291/1‑50 (MQ=255)
                   cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg                          >  1:1624072/1‑50 (MQ=255)
                   cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg                          >  1:2390179/1‑50 (MQ=255)
                      aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                       >  1:2960193/1‑50 (MQ=255)
                      aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                       >  1:793117/1‑50 (MQ=255)
                      aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                       >  1:1396223/1‑50 (MQ=255)
                       aCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt                       >  1:1921142/1‑49 (MQ=255)
                       aCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt                      <  1:2088122/50‑1 (MQ=255)
                        ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt                      >  1:168041/1‑49 (MQ=255)
                        ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTa                     <  1:5011315/50‑1 (MQ=255)
                        ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTa                     <  1:885595/50‑1 (MQ=255)
                         cTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTaa                    <  1:1217759/50‑1 (MQ=255)
                            ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                  >  1:3792290/1‑49 (MQ=255)
                            ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                  <  1:2912575/49‑1 (MQ=255)
                            ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                  <  1:4690584/49‑1 (MQ=255)
                             gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                  >  1:4753916/1‑48 (MQ=255)
                               aCGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt              >  1:4279349/1‑50 (MQ=255)
                                cGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt              <  1:3601995/49‑1 (MQ=255)
                                cGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt              <  1:2851971/49‑1 (MQ=255)
                                 gCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt              >  1:4697761/1‑48 (MQ=255)
                                  cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg           >  1:2664602/1‑50 (MQ=255)
                                  cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg           <  1:2221849/50‑1 (MQ=255)
                                  cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg           <  1:173405/50‑1 (MQ=255)
                                     tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg           >  1:894701/1‑47 (MQ=255)
                                     tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc          <  1:114677/48‑1 (MQ=255)
                                     tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc          <  1:4263760/48‑1 (MQ=255)
                                     tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc          <  1:4173665/48‑1 (MQ=255)
                                      aCCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa         >  1:2353051/1‑48 (MQ=255)
                                       ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag      >  1:2590833/1‑50 (MQ=255)
                                       ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag      >  1:255280/1‑50 (MQ=255)
                                       ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag      >  1:1032281/1‑50 (MQ=255)
                                        cAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc     >  1:284413/1‑50 (MQ=255)
                                          gATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg   <  1:2156110/50‑1 (MQ=255)
                                          gATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg   <  1:905751/50‑1 (MQ=255)
                                           aTCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa  <  1:2812032/50‑1 (MQ=255)
                                            tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa  >  1:4411051/1‑49 (MQ=255)
                                            tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa  >  1:3647611/1‑49 (MQ=255)
                                               |                                             
GGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGA  >  NC_000913/3440147‑3440239

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