Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,270,811 C→T V347M (GTG→ATG)  chaA ← Na(+)/K(+):H(+) antiporter ChaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,270,8110CT99.0% 370.4 / ‑6.0 96V347M (GTG→ATG) chaANa(+)/K(+):H(+) antiporter ChaA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (52/43);  minor base A (0/1);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCACGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTAAAGAG  >  NC_000913/1270764‑1270857
                                               |                                              
catcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGt                                              >  1:6725211/1‑50 (MQ=255)
 atcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGt                                              >  1:6610082/1‑49 (MQ=255)
 atcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGt                                              >  1:4134716/1‑49 (MQ=255)
 atcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGt                                              <  1:8220475/49‑1 (MQ=255)
 atcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGt                                              <  1:4430624/49‑1 (MQ=255)
  tcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGt                                              >  1:2094273/1‑48 (MQ=255)
  tcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTa                                            <  1:200248/50‑1 (MQ=255)
  tcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTa                                            <  1:3538306/50‑1 (MQ=255)
   caAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAg                                           >  1:6648465/1‑50 (MQ=255)
   caAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAg                                           >  1:7432068/1‑50 (MQ=255)
    aaaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAg                                           >  1:5065427/1‑49 (MQ=255)
    aaaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAg                                           <  1:2710094/49‑1 (MQ=255)
    aaaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAg                                           <  1:76483/49‑1 (MQ=255)
    aaaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGt                                          <  1:5130739/50‑1 (MQ=255)
     aaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTa                                         <  1:356625/50‑1 (MQ=255)
     aaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTa                                         >  1:2016490/1‑50 (MQ=255)
     aaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTa                                         >  1:2788625/1‑50 (MQ=255)
      aTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTAc                                        <  1:4124557/50‑1 (MQ=255)
       tAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACg                                       <  1:5064123/50‑1 (MQ=255)
       tAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACg                                       >  1:840718/1‑50 (MQ=255)
       tAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACg                                       >  1:5880622/1‑50 (MQ=255)
       tAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACg                                       <  1:2339154/50‑1 (MQ=255)
       tAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACg                                       <  1:7807680/50‑1 (MQ=255)
        aGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGt                                      >  1:7404381/1‑50 (MQ=255)
         ggcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGt                                      <  1:3686104/49‑1 (MQ=255)
         ggcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGt                                      <  1:4039682/49‑1 (MQ=255)
         ggcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTc                                     <  1:6269848/50‑1 (MQ=255)
         ggcggcAAACAGTGCCAGATGCGCAGCGCCATTGAGCATGTTAGTACGTc                                     <  1:2267847/50‑1 (MQ=38)
           cggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCg                                   >  1:2692667/1‑50 (MQ=255)
           cggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCg                                   >  1:2630025/1‑50 (MQ=255)
            ggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCgg                                  >  1:626356/1‑50 (MQ=255)
            ggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCgg                                  >  1:8348103/1‑50 (MQ=255)
            ggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCgg                                  >  1:1389014/1‑50 (MQ=255)
             gcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGt                                 <  1:7878798/50‑1 (MQ=255)
             gcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGt                                 >  1:6239641/1‑50 (MQ=255)
              cAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  1:6636479/1‑50 (MQ=255)
              cAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  1:6540345/1‑50 (MQ=255)
              cAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  1:5571857/1‑50 (MQ=255)
               aaaCAGTGCCAGATGCGCTGCGCCCTTGAGCATGTTAGTACGTCCGGTg                                >  1:5322276/1‑49 (MQ=39)
               aaaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTgg                               <  1:5786354/50‑1 (MQ=255)
               aaaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  1:4251344/1‑49 (MQ=255)
               aaaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  1:5054100/1‑49 (MQ=255)
                aaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGga                              <  1:5116369/50‑1 (MQ=255)
                 aCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  1:3756089/1‑47 (MQ=255)
                 aCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  1:4905032/1‑47 (MQ=255)
                  cAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgaga                            <  1:3806313/50‑1 (MQ=255)
                  cAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgag                             >  1:6729738/1‑49 (MQ=255)
                   aGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGaa                           <  1:7226948/50‑1 (MQ=255)
                    gTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          <  1:4732531/50‑1 (MQ=255)
                     tGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAgg                         <  1:1149128/50‑1 (MQ=255)
                     tGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAgg                         <  1:5901308/50‑1 (MQ=255)
                      gCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGGGAAg                          >  1:1695026/1‑48 (MQ=39)
                       ccAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  1:560416/1‑47 (MQ=255)
                       ccAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  1:5558725/1‑47 (MQ=255)
                       ccAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgag                       <  1:8138588/50‑1 (MQ=255)
                       ccAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgag                       <  1:3810492/50‑1 (MQ=255)
                        cAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  1:8267932/1‑46 (MQ=255)
                        cAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  1:8444456/1‑46 (MQ=255)
                        cAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      <  1:550521/50‑1 (MQ=255)
                        cAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      <  1:2440683/50‑1 (MQ=255)
                        cAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      <  1:6618200/50‑1 (MQ=255)
                        cAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      <  1:8146426/50‑1 (MQ=255)
                         aGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      >  1:3235324/1‑49 (MQ=255)
                          gATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGata                    >  1:6222837/1‑50 (MQ=255)
                          gATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGata                    <  1:3713591/50‑1 (MQ=255)
                           aTGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGatat                   <  1:6645090/50‑1 (MQ=255)
                           aTGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGatat                   <  1:5685814/50‑1 (MQ=255)
                           aTGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGatat                   >  1:2574442/1‑50 (MQ=255)
                           aTGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGatat                   <  1:2398368/50‑1 (MQ=255)
                            tgcgctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATg                  >  1:7016143/1‑50 (MQ=255)
                            tgcgctgcgcCATTGAGCATGTTAGTACGTCCGGCGGAGAAGGAGATATg                  >  1:2961746/1‑50 (MQ=255)
                             gcgctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATgg                 >  1:2695426/1‑50 (MQ=255)
                              cgctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATgg                 >  1:2991773/1‑49 (MQ=255)
                               gctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGca               >  1:4666934/1‑50 (MQ=255)
                                ctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGcac              <  1:5144845/50‑1 (MQ=255)
                                 tgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGca               >  1:1702335/1‑48 (MQ=255)
                                 tgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGca               >  1:6602883/1‑48 (MQ=255)
                                 tgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGca               >  1:4667729/1‑48 (MQ=255)
                                  gcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGCGATATGGCACAg            >  1:8617267/1‑50 (MQ=255)
                                  gcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAg            >  1:2276112/1‑50 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGcaca             <  1:3281085/48‑1 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAg            >  1:521186/1‑49 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAg            >  1:4912346/1‑49 (MQ=255)
                                   cgcCATTGAACAAGTTAGTACGTCCGGTGGAGAAGGAGATATGGcaca             <  1:4500279/48‑1 (MQ=39)
                                    gcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCa          >  1:6156901/1‑50 (MQ=255)
                                       aTTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACt        >  1:510328/1‑49 (MQ=255)
                                       aTTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTa       <  1:4542165/50‑1 (MQ=255)
                                        ttGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTaa      <  1:4468700/50‑1 (MQ=255)
                                        ttGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTaa      <  1:3093786/50‑1 (MQ=255)
                                          gAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTaaa     <  1:9101/49‑1 (MQ=255)
                                          gAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTaaa     <  1:6879288/49‑1 (MQ=255)
                                          gAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTaaa     <  1:4745506/49‑1 (MQ=255)
                                           aGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTaaa     >  1:2926046/1‑48 (MQ=255)
                                           aGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTaaa     >  1:6091680/1‑48 (MQ=255)
                                           aGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTAAag    <  1:3634551/49‑1 (MQ=255)
                                            gCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTAAagag  >  1:5319958/1‑50 (MQ=255)
                                               |                                              
CATCAAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCACGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTAAAGAG  >  NC_000913/1270764‑1270857

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