Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,916,216 G→A G387E (GGG→GAG)  barA → sensor histidine kinase BarA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,916,2160GA91.9% 93.1 / ‑0.1 37G387E (GGG→GAG) barAsensor histidine kinase BarA
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (16/18);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGGGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAACGCGATCAAT  >  NC_000913/2915969‑2916440
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                
cGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGGGtt                                                                                                                                                                                                                               >  1:393498/1‑251 (MQ=255)
              tGGTGTTATTGGCTTTACCCGCCTGACGCTGACAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTCTTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGTTAAAGAGTTAGAACTGACGCTc                                                                                                                                                                                                                  >  2:466639/1‑250 (MQ=255)
                   ttATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATAtt                                                                                                                                                                                                            >  2:332158/1‑251 (MQ=255)
                                  gCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCt                                                                                                                                                                                             <  2:422333/251‑1 (MQ=255)
                                        cGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAAc                                                                                                                                                                                       <  2:113008/251‑1 (MQ=255)
                                           tGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTg                                                                                                                                                                                    >  1:409931/1‑251 (MQ=255)
                                              aaaCAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATc                                                                                                                                                                                 <  1:29482/251‑1 (MQ=255)
                                                   gAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCg                                                                                                                                                                             >  1:258957/1‑250 (MQ=255)
                                                       tAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCg                                                                                                                                                                        >  1:340101/1‑251 (MQ=255)
                                                                        gATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTATCGAAACTGGAAGCAGGTAAGCTTATTCTGGAAAGTATTCCATTCACACTACGCAGCACGCTGGATGACGTCGTTACTCTGCTGGCACATTCTTCTCACGATAATGAGTTAGAACTGACGCTCCATATTACAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAat                                                                                                                                                       >  1:229870/1‑251 (MQ=255)
                                                                         aTCACCTGAATACGATTGAACGTTCGGCAAATAAGTTGTTGGCAATTATTAATGATGTTCTCGACGTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTAGTCCATTCCCACTACGCAGCACGCTGGATGAAGTCTTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGTGACCCGCTGCGATTACAGCAAatc                                                                                                                                                      <  1:23594/251‑1 (MQ=255)
                                                                                gAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTa                                                                                                                                               <  1:119788/251‑1 (MQ=255)
                                                                                                tCGGCAAATAATGTGCTGGCAATTAGTAATGATGTTCTCGTCTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCGGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATg                                                                                                                                <  1:325917/250‑1 (MQ=255)
                                                                                                   gCAAATAATTTGCTGGCAATTATTAATGTTGTTCTCTACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAAGCATCACTAACCTGGTGGGGAATGCAAt                                                                                                                            <  1:356261/251‑1 (MQ=255)
                                                                                                    cAAATAAGTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAAtt                                                                                                                           <  1:444989/251‑1 (MQ=255)
                                                                                                     aaaTAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCAATTCCCACTACGCAGCACGCTGGATGAAGTAGTTAATCTGCTGGCACATTCTTCTCACGATAAAGGGATAGAACTGACGCTAAATATTAAACGCGACGAGCCTGATAACGTGACCGGAGACCCGCTGCGCTTCCAGCAAATCATCACTAACCTGGTGGGGAATGAAAtt                                                                                                                           >  2:418867/1‑250 (MQ=255)
                                                                                                                  gCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGAATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATAATCACTAACCTGGTGGGGAATGCAATTAAATTCCCCgaga                                                                                                              >  2:327544/1‑250 (MQ=255)
                                                                                                                        attaataaggttttcgaaatctcGAAACTGGCAGCAGGAAAGCTGTTTCTGGATTGAATTCCATTCCCACTACGCAGCAAGCTGGAGGAAGTCGTGACTCTTCTGGCACATTCTTCTCACTCTAAAGAGTTAGAACTGGCGCTCCATATTAAAAGCGACGTGCCTGATGACGTGATCGGCGACCCGCTGCGGTTACAGCAAATCATCACTAACCTGGTGGGGCGTGCACTTAACTTCACCgagactgaca                                                                                                        >  2:75151/20‑244 (MQ=255)
                                                                                                                        attaATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCaa                                                                                                       <  1:79207/251‑1 (MQ=255)
                                                                                                                                      cGACTTCCCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTgg                                                                                         <  2:69597/251‑1 (MQ=255)
                                                                                                                                                          cAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCt                                                                                                                                                                                             >  2:193562/1‑131 (MQ=255)
                                                                                                                                                          cAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCt                                                                                                                                                                                             <  1:193562/131‑1 (MQ=255)
                                                                                                                                                                     aTTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCa                                                                                                                    >  1:239318/1‑193 (MQ=255)
                                                                                                                                                                     aTTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCa                                                                                                                    <  2:239318/193‑1 (MQ=255)
                                                                                                                                                                          ggAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAg                                                      <  1:332158/250‑1 (MQ=255)
                                                                                                                                                                             aaGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGGCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTg                                                  <  2:88341/251‑1 (MQ=255)
                                                                                                                                                                                 attccattccCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGGGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAg                                               <  2:195872/250‑1 (MQ=255)
                                                                                                                                                                                  ttccattccCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCCGATTGAAGTg                                             >  1:92924/1‑251 (MQ=255)
                                                                                                                                                                                    ccattccCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCa                                           <  2:482777/251‑1 (MQ=255)
                                                                                                                                                                                           cACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCgg                                    >  2:224912/1‑251 (MQ=255)
                                                                                                                                                                                                cGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATAc                               >  2:184349/1‑251 (MQ=255)
                                                                                                                                                                                                     caCGCTGGATGAAGGCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCGCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCa                          <  2:472373/251‑1 (MQ=255)
                                                                                                                                                                                                       cGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCAtt                        >  2:98522/1‑251 (MQ=255)
                                                                                                                                                                                                              tGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGCGGGGAATGCAATTAAATTCACCGCGAGTGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTc                 >  1:286357/1‑251 (MQ=255)
                                                                                                                                                                                                                       tACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAACACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAAcgcg        >  2:72278/1‑251 (MQ=255)
                                                                                                                                                                                                                        cctctGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCTACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAACGCGa       <  1:214950/250‑1 (MQ=255)
                                                                                                                                                                                                                             gctgGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAACGCGATCAAt  <  1:224912/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                
CGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGGGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAACGCGATCAAT  >  NC_000913/2915969‑2916440

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