Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 321,742 T→C 11.1% intergenic (+27/+431) aguA → / ← PP_0267 agmatine deiminase/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947321,7420TC11.1% 105.8 / 5.4 45intergenic (+27/+431)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (21/19);  new base C (3/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

GCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGG  >  NC_002947/321602‑321887
                                                                                                                                            |                                                                                                                                                 
gCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGc                                                                                                                                          <  2:427636/150‑1 (MQ=255)
ccctggATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCAGGCGCCGGTCAAGCGCTGAGCCAGTAATGAATGCGCGGCCCGTCATGTGACGGGc                                                                                                                                          >  1:26054/5‑150 (MQ=255)
  aaGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCct                                                                                                                                        >  1:732509/1‑150 (MQ=255)
             ccGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCttt                                                                                                                                    >  2:424684/1‑143 (MQ=255)
             ccGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCttt                                                                                                                                    <  1:424684/143‑1 (MQ=255)
              cGACCATGAGTTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCCAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGcgcg                                                                                                                            >  1:199509/1‑150 (MQ=255)
                    tGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGAc                                                                                                                      >  1:110907/1‑150 (MQ=255)
                        gtggtgATGATCCCGGGCCGCAAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCAGTAATGAATGAGAGGCCCGTCATGTGACGGGCATCTTTTTTGCGCGGGCGACGcac                                                                                                                  >  2:176665/1‑148 (MQ=255)
                           gtgATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGACCGTAACATATGTGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCAGGCGGCGCGCGGa                                                                                                               >  2:551864/1‑150 (MQ=255)
                           gtgATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACACAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGa                                                                                                               >  2:666105/1‑150 (MQ=255)
                            tgatgaTCCCGGGCCGCGTGTATCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCACCCGGCGCCGGTCAAGCGCTGAGCCCGTAATTAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGaa                                                                                                              <  1:413971/150‑1 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCttt                                                                                                                                    >  1:322097/1‑124 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCttt                                                                                                                                    <  2:322097/124‑1 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCg                                                                                                          >  2:438975/1‑150 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCACATGACgg                                                                                                                                            <  2:408112/116‑1 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCACATGACgg                                                                                                                                            >  1:408112/1‑116 (MQ=255)
                                  tCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAg                                                                                                        >  1:89806/1‑150 (MQ=255)
                                   cccGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGc                                                                                                       <  2:89806/150‑1 (MQ=255)
                                     cGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTg                                                                                                     >  1:104034/1‑150 (MQ=255)
                                      gggCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTgg                                                                                                    <  1:454845/150‑1 (MQ=255)
                                             gAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGt                                                                                             >  2:37688/1‑150 (MQ=255)
                                                  gctCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGt                                                                                        >  1:457984/1‑150 (MQ=255)
                                                   ctctTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGAGACGGTCTAGTGCTGATCCCGTAAAAAAAGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTa                                                                                       <  1:355934/150‑1 (MQ=255)
                                                       tGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGtt                                                                                  <  2:385237/150‑1 (MQ=255)
                                                                 cAACATTCACTTCCTGACCCAACCGCATCCGGCGCCGGTCAATCCCTTAGCCCGTAATGAATTAGAGGCCCGTCATGTTACGGGCCTCTTTTTTACGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCa                                                                         <  1:176665/150‑1 (MQ=255)
                                                                 cAACATTCACTGCCTGACCCAACAGAAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAAAAAACGAGAGACCAGTCATGTGAAGGCACACTTTTTTGCGTGGTCGACGCCTTGTAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCa                                                                         <  2:273966/150‑1 (MQ=255)
                                                                         aCTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCaa                                                                 >  1:326782/1‑150 (MQ=255)
                                                                             ccTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATg                                                             >  2:675512/1‑150 (MQ=255)
                                                                                   ccAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGAACTGCCTCCGCCGATCTCCTTGAGGCCTGCCGCGGGGGGTCGCGGGCGACGACCGGAAGGAGAGCTGTCATACGTTTTGTATTGTTTTTGAAGCAGATAGCAAGATGGCACgg                                                       >  1:578496/1‑150 (MQ=255)
                                                                                            gcCGGCGCCGGTCGCGAGCTGCGAACGTAATCTATTAGACGCCCTTCTTTTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCcgc                                              <  2:457984/150‑1 (MQ=255)
                                                                                                     ggTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGAt                                     <  2:266320/150‑1 (MQ=255)
                                                                                                          aGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTa                                >  2:380838/1‑150 (MQ=255)
                                                                                                           gcgcTGAGCACGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCGGGAGGCGGAGCTGGCATGAGTTTTGTATTGTTTTCCCCGCAGAGAGCAACATAGCACGGTCAGGCCGCGCCTCTGATCTGTaa                               <  2:26054/150‑1 (MQ=255)
                                                                                                              cTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACaa                            <  1:534853/150‑1 (MQ=255)
                                                                                                                     cGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTAc                     <  1:337145/150‑1 (MQ=255)
                                                                                                                           gaatgaGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAAt               <  2:473154/150‑1 (MQ=255)
                                                                                                                                 gagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGtttt                                                                                >  2:533820/1‑79 (MQ=38)
                                                                                                                                 gagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGtttt                                                                                <  1:533820/79‑1 (MQ=255)
                                                                                                                                  agGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGc        >  2:399010/1‑150 (MQ=255)
                                                                                                                                  agGCCCGTCACATGGCGGGCCTTGTTTTGGGGCGTGCGACTGCCGGAAGCCGAGCTGGCGTACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATCGCACGGCCAGCCCTCGTCTCTGATCTTTAACAATACCTACGTAAATTAGCCGc        >  2:163958/1‑150 (MQ=255)
                                                                                                                                      ccGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTCCGTAAATTAGCCGCTCAc    >  1:118583/1‑150 (MQ=255)
                                                                                                                                      ccGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCAc    <  1:37688/150‑1 (MQ=255)
                                                                                                                                       cGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTaa                               >  2:273139/1‑122 (MQ=255)
                                                                                                                                       cGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTaa                               <  1:273139/122‑1 (MQ=255)
                                                                                                                                        gTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACgg  <  1:399010/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGG  >  NC_002947/321602‑321887

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