Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,291,084 A→C K195N (AAA→AAC narP → response regulator in two‑component regulatory system with NarQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,291,0840AC98.7% 261.7 / ‑3.8 74K195N (AAA→AACnarPresponse regulator in two‑component regulatory system with NarQ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (41/32);  total (42/32)
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.62e-01

ACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAG  >  NC_000913/2290987‑2291179
                                                                                                 |                                                                                               
aCGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCt                                                                                               <  1:543665/100‑1 (MQ=255)
  gAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCa                                                                                             <  2:226634/100‑1 (MQ=255)
    gCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAAt                                                                                           >  2:69694/1‑100 (MQ=255)
     cTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATg                                                                                          <  1:1405476/100‑1 (MQ=255)
       ggCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTc                                                                                        >  1:741251/1‑100 (MQ=255)
       ggCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTc                                                                                        >  2:993683/1‑100 (MQ=255)
        gCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTcc                                                                                       >  2:1002997/1‑100 (MQ=255)
         cacaGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCg                                                                                      >  2:812070/1‑100 (MQ=255)
          acaGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGc                                                                                     >  1:744627/1‑100 (MQ=255)
          acaGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGc                                                                                     >  1:742948/1‑100 (MQ=255)
               ggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCAcg                                                                                >  1:198896/1‑100 (MQ=255)
                gCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCAcgc                                                                               <  2:1052195/100‑1 (MQ=255)
                   gTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTg                                                                            >  2:707492/1‑100 (MQ=255)
                    tCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTgg                                                                           >  2:737468/1‑100 (MQ=255)
                    tCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTgg                                                                           >  2:965246/1‑100 (MQ=255)
                     cAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTggc                                                                          >  1:318453/1‑100 (MQ=255)
                     cAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTggc                                                                          >  2:1066248/1‑100 (MQ=255)
                     cAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTggc                                                                          >  2:816653/1‑100 (MQ=255)
                          aaaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGcca                                                                     >  2:1472734/1‑100 (MQ=255)
                           aaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGccac                                                                    >  2:321531/1‑100 (MQ=255)
                           aaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGccac                                                                    <  2:351321/100‑1 (MQ=255)
                            aCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGccacc                                                                   <  1:1560496/100‑1 (MQ=255)
                            aCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGccacc                                                                   >  2:799542/1‑100 (MQ=255)
                             cAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGccacca                                                                  >  1:291303/1‑100 (MQ=255)
                              aGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCAt                                                                 >  2:1473815/1‑100 (MQ=255)
                               gATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCAtt                                                                <  2:854022/100‑1 (MQ=255)
                               gATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCAtt                                                                >  1:1488437/1‑100 (MQ=255)
                                aTTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTACCCTCAATGTCCGCTCACGCGTGGCGGCCACCATTc                                                               <  2:1146627/100‑1 (MQ=255)
                                aTTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTc                                                               >  1:192372/1‑100 (MQ=255)
                                  tGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTg                                                             >  2:1240726/1‑100 (MQ=255)
                                   gCCTCGGTGTTGAATATGTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCt                                                              <  2:54533/98‑1 (MQ=255)
                                    ccgatctgtTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGtt                                                           <  2:1290867/94‑1 (MQ=255)
                                      tCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTcc                                                         <  1:581476/100‑1 (MQ=255)
                                          tgtTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGca                                                     <  1:816658/100‑1 (MQ=255)
                                           gtTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGttagatcgg                                                    >  1:1290876/1‑93 (MQ=255)
                                             tGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaaca                                                  <  1:799547/100‑1 (MQ=255)
                                              gAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACGTCAATGTCCGCTCACCCGTGGCGGACACCATTCTGTTCCTGcaacaa                                                 >  2:1061347/1‑100 (MQ=255)
                                              gAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaacaa                                                 >  2:735423/1‑100 (MQ=255)
                                               aaTATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaacaa                                                 >  2:1428421/1‑99 (MQ=255)
                                                atatTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAAcg                                               >  2:58448/1‑100 (MQ=255)
                                                  atTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAACGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAAcgcg                                             <  2:553623/100‑1 (MQ=255)
                                                   tttCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCGCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCgg                                            >  1:601184/1‑100 (MQ=255)
                                                      ccGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGc                                         <  2:33697/100‑1 (MQ=255)
                                                       cGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGca                                        <  1:176876/100‑1 (MQ=255)
                                                       cGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGc                                         <  2:1198568/99‑1 (MQ=255)
                                                       cGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTCCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGca                                        >  1:1105786/1‑100 (MQ=255)
                                                         aGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGcaca                                      <  2:1444866/100‑1 (MQ=255)
                                                          gCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACaa                                     >  1:600843/1‑100 (MQ=255)
                                                           cAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCATGGCTGCCACCATTCTGTTCCTGCAACCACGCGGGGcacca                                     >  2:1094155/1‑97 (MQ=255)
                                                           cAGACAGTAAAAGGACATATTCGCAATCTGCTGCGTAACCTCAATCTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAAt                                    <  2:1070705/100‑1 (MQ=255)
                                                             gACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaa                                  >  1:1434005/1‑100 (MQ=255)
                                                             gACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATa                                   >  2:979262/1‑99 (MQ=255)
                                                              aCAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaaa                                 >  2:1201864/1‑100 (MQ=255)
                                                               cAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAAt                                    <  2:599025/96‑1 (MQ=255)
                                                               cAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaaaa                                <  1:1270584/100‑1 (MQ=255)
                                                                 gTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAAt                              >  2:1384834/1‑100 (MQ=255)
                                                                 gTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAAt                              >  1:875114/1‑100 (MQ=255)
                                                                  tAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATa                             <  1:1003004/100‑1 (MQ=255)
                                                                   aaaaGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAg                            <  1:1121126/100‑1 (MQ=255)
                                                                     aaGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGcc                          <  1:69694/100‑1 (MQ=255)
                                                                      aGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGccc                         >  1:1341434/1‑100 (MQ=255)
                                                                        tACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGa                       <  1:1240734/100‑1 (MQ=255)
                                                                           atatTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATgg                    <  1:321533/100‑1 (MQ=255)
                                                                               tCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGAttt                >  2:451891/1‑100 (MQ=255)
                                                                                cGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTa               <  1:216414/100‑1 (MQ=255)
                                                                                cGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTa               <  1:1431807/100‑1 (MQ=255)
                                                                                cGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTa               >  2:33969/1‑100 (MQ=255)
                                                                                 gCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTat              <  2:1070349/100‑1 (MQ=255)
                                                                                      ctgctgCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTACCATCg         >  2:537166/1‑100 (MQ=255)
                                                                                       tgctgcGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCgg        <  1:1299807/100‑1 (MQ=255)
                                                                                         ctgcGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGc      >  1:1560379/1‑100 (MQ=255)
                                                                                           gcGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTg    >  2:1444277/1‑100 (MQ=255)
                                                                                            cGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAACAATAGCCCGATGGATTTATCATCGGGCTga   <  2:448200/100‑1 (MQ=255)
                                                                                             gTAACCTCATTGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCTGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTgag  <  2:313029/100‑1 (MQ=255)
                                                                                                 |                                                                                               
ACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAG  >  NC_000913/2290987‑2291179

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