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,8110CT96.8% 471.5 / 3.1 123V347M (GTG→ATG) chaANa(+)/K(+):H(+) antiporter ChaA
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base T (51/68);  total (51/72)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 123 total aligned reads displayed.

CATCAAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCACGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTAAAGAG  >  NC_000913/1270764‑1270857
                                               |                                              
catcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGt                                              >  1:2298283/1‑50 (MQ=255)
 atcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGtt                                             <  2:875515/50‑1 (MQ=255)
 atcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGtt                                             <  1:4908838/50‑1 (MQ=255)
 atcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGtt                                             <  1:2789992/50‑1 (MQ=255)
 atcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCACGtt                                             <  2:4522181/50‑1 (MQ=255)
 atcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCACGt                                              <  2:1253415/49‑1 (MQ=255)
  tcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGt                                              >  1:3233530/1‑48 (MQ=255)
  tcaAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTa                                            <  2:3786603/50‑1 (MQ=255)
   caAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAg                                           >  1:1366763/1‑50 (MQ=255)
   caAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAg                                           >  1:2971453/1‑50 (MQ=255)
   caAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAg                                           >  2:439915/1‑50 (MQ=255)
    aaaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGt                                          <  2:4999137/50‑1 (MQ=255)
    aaaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGt                                          <  2:2119067/50‑1 (MQ=255)
    aaaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGt                                          <  2:3972234/50‑1 (MQ=255)
    aaaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGt                                          <  2:1927392/50‑1 (MQ=255)
     aaTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTa                                         <  1:4732544/50‑1 (MQ=255)
      aTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTAc                                        <  2:3501848/50‑1 (MQ=255)
      aTAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTAc                                        <  1:2747967/50‑1 (MQ=255)
       tAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACg                                       >  2:4881509/1‑50 (MQ=255)
       tAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACg                                       <  1:3270016/50‑1 (MQ=255)
         ggcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTc                                     >  2:5073264/1‑50 (MQ=255)
         ggcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTc                                     <  1:4266696/50‑1 (MQ=255)
         ggcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTc                                     >  1:4053676/1‑50 (MQ=255)
         ggcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTc                                     <  1:2982714/50‑1 (MQ=255)
          gcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTcc                                    <  2:4966647/50‑1 (MQ=255)
          gcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTcc                                    <  2:3113359/50‑1 (MQ=255)
          gcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTcc                                    <  1:4933719/50‑1 (MQ=255)
          gcggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTcc                                    <  1:4857585/50‑1 (MQ=255)
           cggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCg                                   >  2:2566174/1‑50 (MQ=255)
           cggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCg                                   >  1:4840055/1‑50 (MQ=255)
            ggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCgg                                  >  1:3135437/1‑50 (MQ=255)
            ggcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCgg                                  >  1:580132/1‑50 (MQ=255)
            ggcAAACAGTGCCAGATGCGCTGCGCCATGGAGCATGTTAGTACGTCCgg                                  <  1:2899849/50‑1 (MQ=39)
             gcAAACAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGt                                 <  1:1954664/50‑1 (MQ=255)
               aaaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTgg                               <  2:2207878/50‑1 (MQ=255)
               aaaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  2:1717299/1‑49 (MQ=255)
               aaaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  1:1809907/1‑49 (MQ=255)
                aaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  1:1514048/1‑48 (MQ=255)
                aaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  2:4811229/1‑48 (MQ=255)
                aaCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGga                              <  1:1026257/50‑1 (MQ=255)
                 aCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTg                                >  2:2852916/1‑47 (MQ=255)
                 aCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgag                             <  1:454344/50‑1 (MQ=255)
                 aCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgag                             <  2:1980780/50‑1 (MQ=255)
                 aCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgag                             <  1:1121914/50‑1 (MQ=255)
                 aCAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgag                             <  2:5071145/50‑1 (MQ=255)
                  cAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgaga                            <  1:1722211/50‑1 (MQ=255)
                  cAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgaga                            <  1:339126/50‑1 (MQ=255)
                  cAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgag                             >  1:4341522/1‑49 (MQ=255)
                  cAGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGgag                             >  2:1635647/1‑49 (MQ=255)
                   aGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGaa                           <  2:1539573/50‑1 (MQ=255)
                   aGTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGaa                           <  1:3495129/50‑1 (MQ=255)
                    gTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  1:5271458/1‑50 (MQ=255)
                    gTGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          <  1:291203/50‑1 (MQ=255)
                     tGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  1:5336331/1‑49 (MQ=255)
                     tGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  2:2203949/1‑49 (MQ=255)
                     tGCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  2:366212/1‑49 (MQ=255)
                      gCCAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  1:1110527/1‑48 (MQ=255)
                       ccAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgag                       <  1:821349/50‑1 (MQ=255)
                       ccAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgag                       <  1:3660692/50‑1 (MQ=255)
                       ccAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGga                        <  1:2203282/49‑1 (MQ=255)
                        cAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAg                          >  2:3274016/1‑46 (MQ=255)
                        cAGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      <  2:2887154/50‑1 (MQ=255)
                         aGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      >  2:830945/1‑49 (MQ=255)
                         aGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      >  1:325076/1‑49 (MQ=255)
                         aGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      >  1:3081267/1‑49 (MQ=255)
                         aGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGgaga                      >  1:1066647/1‑49 (MQ=255)
                         aGATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGat                     <  2:3216334/50‑1 (MQ=255)
                          gATGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGata                    <  2:3940335/50‑1 (MQ=255)
                           aTGCGCTGCGCCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGatat                   >  1:4450641/1‑50 (MQ=255)
                            tgcgctgcgcCATTGAGCATGTTAGTCCGTCCGGTGGAGAAGGAGATATg                  <  2:3360415/50‑1 (MQ=39)
                            tgcgctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATg                  <  1:138511/50‑1 (MQ=255)
                            tgcgctgcgcCATTGAGCACGTTAGTACGTCCGGTGGAGAAGGAGATATg                  <  2:2790038/50‑1 (MQ=255)
                              cgctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATgg                 >  2:3857011/1‑49 (MQ=255)
                              cgctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATgg                 >  1:1055685/1‑49 (MQ=255)
                              cgctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGc                <  1:1346805/50‑1 (MQ=255)
                              cgctgcgcCATTGAGCACGTTAGTACGTCCGGTGGAGAAGGAGATATgg                 <  2:1341795/49‑1 (MQ=255)
                               gctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGca               <  1:3636161/50‑1 (MQ=255)
                               gctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGca               >  1:641424/1‑50 (MQ=255)
                               gctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGca               >  1:476902/1‑50 (MQ=255)
                                ctgcgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGca               >  2:1429169/1‑49 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGcaca             <  1:3473932/48‑1 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGcaca             <  2:601624/48‑1 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGcaca             <  1:1914829/48‑1 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGcaca             <  1:2308522/48‑1 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAg            >  2:3901218/1‑49 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAg            <  2:1975291/49‑1 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAg            >  1:1836792/1‑49 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAg            >  1:935529/1‑49 (MQ=255)
                                   cgcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAg            >  2:1189601/1‑49 (MQ=255)
                                    gcCATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCa          >  1:5384659/1‑50 (MQ=255)
                                     ccATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCAc         <  1:2691442/50‑1 (MQ=255)
                                     ccATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCAc         >  1:4235650/1‑50 (MQ=255)
                                      cATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACt        <  2:4826366/50‑1 (MQ=255)
                                      cATTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACt        <  1:4095781/50‑1 (MQ=255)
                                       aTTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACt        >  1:2595171/1‑49 (MQ=255)
                                       aTTGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTa       <  2:3349612/50‑1 (MQ=255)
                                        ttGAGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTaa      <  2:685031/50‑1 (MQ=255)
                                           aGCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTAAaga   <  2:3161139/50‑1 (MQ=255)
                                            gCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTAAagag  <  1:5158100/50‑1 (MQ=255)
                                            gCATGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTAAagag  <  2:3608041/50‑1 (MQ=255)
                                               |                                              
CATCAAATAGGCGGCAAACAGTGCCAGATGCGCTGCGCCATTGAGCACGTTAGTACGTCCGGTGGAGAAGGAGATATGGCACAGCACTAAAGAG  >  NC_000913/1270764‑1270857

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