Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,291,084 A→C 51.4% 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,0840AC51.4% 12.2 / 98.1 73K195N (AAA→AACnarPresponse regulator in two‑component regulatory system with NarQ
Reads supporting (aligned to +/- strand):  ref base A (20/15);  new base C (23/14);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

ACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATT  >  NC_000913/2290987‑2291182
                                                                                                 |                                                                                                  
aCGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACt                                                                                                  <  2:101742/100‑1 (MQ=255)
  gAGGTGGCGGAGGGGGTGTCAAATAAACAGAGTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCt                                                                                                  <  2:2146694/98‑1 (MQ=255)
    gCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAAt                                                                                              <  1:578891/100‑1 (MQ=255)
     cTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATg                                                                                             <  2:2317731/100‑1 (MQ=255)
     cTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATg                                                                                             <  2:368879/100‑1 (MQ=255)
         cacaGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCg                                                                                         >  2:318950/1‑100 (MQ=255)
             tgggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCa                                                                                     >  2:273287/2‑100 (MQ=255)
             ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCa                                                                                     >  1:302372/1‑100 (MQ=255)
             ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCa                                                                                     <  2:947050/100‑1 (MQ=255)
             ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCa                                                                                     >  2:1977532/1‑100 (MQ=255)
                gCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCAcgc                                                                                  >  2:1723293/1‑100 (MQ=255)
                gCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAcgc                                                                                  >  2:1635645/1‑100 (MQ=255)
                gCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAcgc                                                                                  >  1:452288/1‑100 (MQ=255)
                 cTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCAcgcg                                                                                 <  2:2327587/100‑1 (MQ=255)
                  tGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGt                                                                                >  2:529159/1‑100 (MQ=255)
                   gTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCAcgc                                                                                  >  2:1551388/1‑97 (MQ=255)
                   gTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTg                                                                               <  1:2239172/100‑1 (MQ=255)
                    tCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGGCCGCTCACGCGTg                                                                               <  2:2160291/99‑1 (MQ=255)
                     cAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTggc                                                                             >  2:2254922/1‑100 (MQ=255)
                     cAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACACGTagc                                                                             >  1:33556/1‑100 (MQ=255)
                     cAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCATggc                                                                             >  1:1311731/1‑100 (MQ=255)
                      aaaTAAACAGATTGCCTCAGAGTTGAATATTACCGAGAAGACAGTAAAAGTACATATTAGCAATCTGATGCGTAATCTCAATGTCCGCTCACGCGTggcg                                                                            >  2:2059391/1‑100 (MQ=255)
                       aaTAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAATTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcgg                                                                           <  1:577571/100‑1 (MQ=255)
                        aTAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcggc                                                                          >  1:1853203/1‑100 (MQ=255)
                         tAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGcc                                                                         >  1:1700685/1‑100 (MQ=255)
                         tAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGcc                                                                         <  2:194922/100‑1 (MQ=255)
                         tAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGcc                                                                         >  2:1587863/1‑100 (MQ=255)
                          aaaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGcca                                                                        <  1:1695992/100‑1 (MQ=255)
                           aaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGccac                                                                       <  2:1551204/100‑1 (MQ=255)
                           aaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGccac                                                                       >  2:1506674/1‑100 (MQ=255)
                           aaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccac                                                                       >  1:2347081/1‑100 (MQ=255)
                           aaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccac                                                                       <  2:1051339/100‑1 (MQ=255)
                              aGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCAt                                                                    >  2:286629/1‑100 (MQ=255)
                               gATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCAtt                                                                   >  2:1963346/1‑100 (MQ=255)
                                aTTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTc                                                                  <  2:1087707/100‑1 (MQ=255)
                                     cTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTc                                                             >  1:831547/1‑100 (MQ=255)
                                     cTCGGTGTTGAATATTTCCAAGCAGACAGTAAAAGTACATATTCGCAATCTGCTCCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTc                                                             <  2:2263259/100‑1 (MQ=255)
                                      tctgtgtTGAATATTTCCGAGCAGACATTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTcc                                                            <  2:2134262/97‑1 (MQ=255)
                                          tgtTGAATATTTCCGAGCAGACAGTAAAAGTATATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGca                                                        <  2:2240115/100‑1 (MQ=255)
                                          tgtTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTg                                                          <  2:2100244/98‑1 (MQ=255)
                                            ttGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaac                                                      >  1:2156330/1‑100 (MQ=255)
                                             tGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaaca                                                     >  1:2248310/1‑100 (MQ=255)
                                             tGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaaca                                                     >  1:274991/1‑100 (MQ=255)
                                              gAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaacaa                                                    >  1:19859/1‑100 (MQ=255)
                                               aaTATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTCGCGGCCACCATTCTGTTCCTGCaacaac                                                   >  1:823294/1‑100 (MQ=255)
                                               aaTATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCaacaac                                                   >  2:985253/1‑100 (MQ=255)
                                                   tttCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCgg                                               >  1:1872554/1‑100 (MQ=255)
                                                        gAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGcac                                          >  2:66944/1‑100 (MQ=255)
                                                             gACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaa                                     >  2:1417349/1‑100 (MQ=255)
                                                             gACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaa                                     >  1:1366621/1‑100 (MQ=255)
                                                             gACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaa                                     >  1:172448/1‑100 (MQ=255)
                                                             gACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaa                                     <  2:1585196/100‑1 (MQ=255)
                                                               cAGTAAAAGTACCTATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaaa                                    <  2:2405784/99‑1 (MQ=255)
                                                               cAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaaaa                                   >  2:1269241/1‑100 (MQ=255)
                                                                 gTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGACCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAAt                                 >  2:2168011/1‑100 (MQ=255)
                                                                  tAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATa                                <  2:1637489/100‑1 (MQ=255)
                                                                   aaaaGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAg                               >  2:743586/1‑100 (MQ=255)
                                                                   aaaaGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAg                               <  1:252564/100‑1 (MQ=255)
                                                                     aaGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGcc                             <  1:1351468/100‑1 (MQ=255)
                                                                      aGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGccc                            >  2:540598/1‑100 (MQ=255)
                                                                         aCATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGAt                         >  1:679016/1‑100 (MQ=255)
                                                                            tatTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGa                      >  1:1585174/1‑100 (MQ=255)
                                                                            tatTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGa                      >  2:1609904/1‑100 (MQ=255)
                                                                              ttCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGAtt                    <  2:1136977/100‑1 (MQ=255)
                                                                              ttCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGa                      >  2:448909/1‑98 (MQ=255)
                                                                              ttCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGAtt                    >  1:1530762/1‑100 (MQ=255)
                                                                               tCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGAttt                   <  2:1872532/100‑1 (MQ=255)
                                                                                cGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTa                  >  2:1998392/1‑100 (MQ=255)
                                                                                 gCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTat                 >  1:1354549/1‑100 (MQ=255)
                                                                                    aTCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTACTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTatc                <  1:946675/98‑1 (MQ=255)
                                                                                      ctgcggCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTACCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCg            <  2:1181864/100‑1 (MQ=255)
                                                                                         ctgcGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGc         <  2:1867451/100‑1 (MQ=255)
                                                                                             gTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTgag     >  2:929606/1‑100 (MQ=255)
                                                                                               aaaCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGAt   >  2:1644749/1‑100 (MQ=255)
                                                                                                aaCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGAtt  >  1:1192937/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
ACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATT  >  NC_000913/2290987‑2291182

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