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,1940CA80.8% 262.4 / 64.6 104G279V (GGT→GTT) rpoARNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (9/11);  new base A (46/38);  total (55/49)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGT  >  NC_000913/3440149‑3440243
                                             |                                                 
ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGa                                                >  1:7386721/1‑49 (MQ=255)
 tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                              >  1:4737971/1‑50 (MQ=255)
  aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTTCCAGATCACCGata                                              <  1:3982944/49‑1 (MQ=255)
  aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata                                              <  1:642141/49‑1 (MQ=255)
  aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                              <  1:4289840/49‑1 (MQ=255)
  aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                              >  1:561102/1‑49 (MQ=255)
   ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGatat                                             <  1:5089900/49‑1 (MQ=255)
   ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata                                              >  1:5012624/1‑48 (MQ=255)
   ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGAtata                                            <  1:7433305/50‑1 (MQ=255)
   ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGAtata                                            <  1:5776539/50‑1 (MQ=255)
   ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGAtata                                            <  1:7090095/50‑1 (MQ=255)
    gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAg                                           <  1:4644914/50‑1 (MQ=255)
      gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTg                                         >  1:780459/1‑50 (MQ=255)
      gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                         >  1:1119441/1‑50 (MQ=255)
      gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                         >  1:4782135/1‑50 (MQ=255)
      gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                         <  1:3101402/50‑1 (MQ=255)
       ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTgg                                        <  1:3976191/50‑1 (MQ=255)
        tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                       >  1:7197372/1‑50 (MQ=255)
        tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                       >  1:4808084/1‑50 (MQ=255)
         ttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGa                                       <  1:690038/49‑1 (MQ=255)
         ttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa                                       <  1:1821217/49‑1 (MQ=255)
          tAAGGGGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGa                                       >  1:2079653/1‑48 (MQ=255)
          tAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGAt                                      <  1:505002/49‑1 (MQ=255)
           aaGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATa                                     <  1:4052779/49‑1 (MQ=255)
           aaGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATa                                     <  1:6396261/49‑1 (MQ=255)
           aaGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAg                                    >  1:34076/1‑50 (MQ=255)
           aaGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAg                                    >  1:3070744/1‑50 (MQ=255)
            aGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGAGGATAGc                                   >  1:4742292/1‑50 (MQ=255)
            aGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAg                                    <  1:2907466/49‑1 (MQ=255)
            aGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGc                                   <  1:7163088/50‑1 (MQ=255)
             ggAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCt                                  >  1:1976996/1‑50 (MQ=255)
             ggAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCt                                  <  1:1422194/50‑1 (MQ=255)
             ggAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCt                                  <  1:7372213/50‑1 (MQ=255)
              gAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCtt                                 <  1:3249438/50‑1 (MQ=255)
              gAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCtt                                 >  1:472136/1‑50 (MQ=255)
               aGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTc                                <  1:2603153/50‑1 (MQ=255)
               aGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTc                                >  1:1839211/1‑50 (MQ=255)
               aGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTc                                >  1:5529286/1‑50 (MQ=255)
               aGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTc                                >  1:1753583/1‑50 (MQ=255)
                gCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCt                               <  1:1617286/50‑1 (MQ=255)
                gCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCt                               <  1:5392403/50‑1 (MQ=255)
                gCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCt                               <  1:706372/50‑1 (MQ=255)
                gCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCt                               <  1:3383779/50‑1 (MQ=255)
                 cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg                              >  1:7029935/1‑50 (MQ=255)
                 cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg                              >  1:1567323/1‑50 (MQ=255)
                  tCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                             >  1:7247029/1‑50 (MQ=255)
                  tCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                             >  1:1769931/1‑50 (MQ=255)
                   cAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCt                            >  1:1109581/1‑50 (MQ=255)
                   cAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCt                            >  1:2496578/1‑50 (MQ=255)
                    aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                           >  1:6092474/1‑50 (MQ=255)
                    aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                           <  1:3234485/50‑1 (MQ=255)
                    aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                           >  1:6320884/1‑50 (MQ=255)
                    aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                           >  1:5308452/1‑50 (MQ=255)
                    aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCt                            <  1:817962/49‑1 (MQ=255)
                     aCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCttt                          <  1:6480250/50‑1 (MQ=255)
                     aCCTCGGTACGCTGTACCAGATCAACGATATAGCGGATAGCTTCTGCtt                           >  1:3351489/1‑49 (MQ=39)
                      ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt                          >  1:686975/1‑49 (MQ=255)
                      ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTa                         <  1:7124656/50‑1 (MQ=255)
                       cTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCttt                          >  1:6862574/1‑48 (MQ=255)
                       cTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt                          >  1:311748/1‑48 (MQ=255)
                        tCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAg                       >  1:5572453/1‑50 (MQ=255)
                         cGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                      <  1:5311768/50‑1 (MQ=255)
                          ggTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg                      <  1:3312137/49‑1 (MQ=255)
                          ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                      <  1:6788705/49‑1 (MQ=255)
                          ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                      <  1:2634644/49‑1 (MQ=255)
                          ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                      <  1:6285663/49‑1 (MQ=255)
                          ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                      >  1:6094225/1‑49 (MQ=255)
                          ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                      >  1:4393776/1‑49 (MQ=255)
                          ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAg                       <  1:6060057/48‑1 (MQ=255)
                          ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGc                     <  1:6886715/50‑1 (MQ=255)
                           gTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa                    <  1:1131683/50‑1 (MQ=255)
                           gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                      >  1:6159799/1‑48 (MQ=255)
                           gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                      >  1:5737798/1‑48 (MQ=255)
                           gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCa                    <  1:4918492/50‑1 (MQ=255)
                            tACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAg                   >  1:6075565/1‑50 (MQ=255)
                            tACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAg                   >  1:1217960/1‑50 (MQ=255)
                            tACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAg                   >  1:5712204/1‑50 (MQ=255)
                            tACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAg                   >  1:2517688/1‑50 (MQ=255)
                              cGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                  <  1:5757786/49‑1 (MQ=255)
                              cGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                  <  1:1777534/49‑1 (MQ=255)
                               gCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                  >  1:7014260/1‑48 (MQ=255)
                               gCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTa                <  1:164170/50‑1 (MQ=255)
                                cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg               <  1:637472/50‑1 (MQ=255)
                                cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg               <  1:3750078/50‑1 (MQ=255)
                                  gTACCAGATCACCGATATAGCGGATAGCTTCTGCTTTAAGGCAGTTAg               >  1:3598890/1‑48 (MQ=255)
                                   tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc              <  1:1000697/48‑1 (MQ=255)
                                    aCCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa             >  1:3012358/1‑48 (MQ=255)
                                     ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag          >  1:46322/1‑50 (MQ=255)
                                     ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag          >  1:346115/1‑50 (MQ=255)
                                     ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCaga           <  1:1939728/49‑1 (MQ=255)
                                      cAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc         >  1:6274128/1‑50 (MQ=255)
                                      cAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc         <  1:3727881/50‑1 (MQ=255)
                                      cAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc         >  1:4983438/1‑50 (MQ=255)
                                       aGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCg        >  1:84688/1‑50 (MQ=255)
                                         aTCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa      >  1:2938916/1‑50 (MQ=255)
                                         aTCAACGATATAGCGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa      >  1:5621308/1‑50 (MQ=255)
                                          tCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGAc     <  1:5630130/50‑1 (MQ=255)
                                          tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa      >  1:2047377/1‑49 (MQ=255)
                                          tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa      >  1:4744272/1‑49 (MQ=255)
                                             ccGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGt  >  1:916458/1‑50 (MQ=255)
                                             |                                                 
TTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGT  >  NC_000913/3440149‑3440243

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