Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 321,740 T→C 8.7% 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,7400TC 8.7% 125.3 / 3.0 46intergenic (+27/+431)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (23/19);  new base C (2/2);  total (25/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 = 7.67e-01

ATCCTGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGG  >  NC_002947/321594‑321885
                                                                                                                                                  |                                                                                                                                                 
aTCCTGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCAtgt                                                                                                                                                 <  6:295002/149‑1 (MQ=255)
 tCCTGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCAtgt                                                                                                                                                 >  8:160489/1‑148 (MQ=255)
   cTGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGAc                                                                                                                                              <  6:22140/149‑1 (MQ=255)
      gCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACgg                                                                                                                                            <  8:194065/148‑1 (MQ=255)
         aGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACgg                                                                                                                                            <  1:5077/145‑1 (MQ=255)
         aGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACgg                                                                                                                                            >  2:5077/1‑145 (MQ=255)
                tCCCAGACCATGTGGTGGTGATGAGCCCGGTCCGGTTGCTTGTCTTGTGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCGGGTCAAGCGGTGAGCCCGTAATGAATGAGAGGCCCTTCATGTGACGGGCCTCttttt                                                                                                                                  <  8:116787/148‑1 (MQ=255)
                  cccGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTgc                                                                                                                               >  5:296648/1‑149 (MQ=255)
                      aCCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCgg                                                                                                                           <  6:296648/149‑1 (MQ=255)
                         aTGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGc                                                                                                                         >  8:152893/1‑148 (MQ=255)
                                       atCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCACCAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGGGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGcc                                                                                                           >  5:356906/1‑148 (MQ=255)
                                        tCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGcc                                                                                                           >  5:98433/1‑147 (MQ=255)
                                        tCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGa                                                                                                         <  6:356906/149‑1 (MQ=255)
                                         cccGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGAc                                                                                                                                              <  6:299061/111‑1 (MQ=255)
                                         cccGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGAc                                                                                                                                              >  5:299061/1‑111 (MQ=255)
                                                     gctgctCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGtt                                                                                            <  6:98433/149‑1 (MQ=255)
                                                                      gcAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAg                                                                           >  6:283058/1‑149 (MQ=255)
                                                                       cAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGc                                                                          <  7:152893/149‑1 (MQ=255)
                                                                       cAACATTCACTGCCAGACCCAACAGCAGCCGGAGCCGGTCAACCGCTGAGCCAGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCATTTTTGCTCGTGCGACTCCCGGAAGCCGAGCTGGCATAAGATTTGTATTTTTTTTCAAg                                                                           >  6:362517/1‑148 (MQ=255)
                                                                          cATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGa                                                                       <  5:283058/149‑1 (MQ=255)
                                                                               aCTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGc                                                                   >  7:248569/1‑148 (MQ=255)
                                                                                         ccAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTAtt                                                                                     >  8:323709/1‑120 (MQ=255)
                                                                                         ccAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACg                                                        >  4:87420/1‑149 (MQ=255)
                                                                                           aaCAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGt                                                      <  8:248569/149‑1 (MQ=255)
                                                                                            acttctgcCGGCGCGGTTCAAGCGCTGAGCCCGTAATGAATGAGTGTCCCGTCATGTGACGGGCCTCTTTTTTGCGCGTGCTACGCACGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCTAGCAGATAGCAAGATGGCACGGTc                                                     <  5:362517/143‑1 (MQ=255)
                                                                                              agcagcCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAg                                                   <  7:160489/149‑1 (MQ=255)
                                                                                                   ccGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCcgc                                              >  4:266160/1‑149 (MQ=255)
                                                                                                     ggCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGAt                                                                      <  8:247468/123‑1 (MQ=255)
                                                                                                     ggCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGAt                                                                      >  7:247468/1‑123 (MQ=255)
                                                                                                            gTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGAt                                     >  2:312132/1‑149 (MQ=255)
                                                                                                                  cgcTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTa                                >  1:44204/1‑148 (MQ=255)
                                                                                                                     tGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCgcgc                                            <  1:251702/133‑1 (MQ=255)
                                                                                                                     tGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCgcgc                                            >  2:251702/1‑133 (MQ=255)
                                                                                                                     tGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACa                             >  2:12388/1‑148 (MQ=255)
                                                                                                                         cccGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATAcc                        >  3:287526/1‑149 (MQ=255)
                                                                                                                          ccGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCGCGGCCTGACCGTGccatct                                       >  5:302676/1‑129 (MQ=255)
                                                                                                                          ccGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCGCGGCCTGACCGTGccatct                                       <  6:302676/133‑5 (MQ=255)
                                                                                                                          ccGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATAcc                        >  8:321295/1‑148 (MQ=255)
                                                                                                                                    tgaGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTa             <  7:321295/149‑1 (MQ=255)
                                                                                                                                       gagGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGt                                                      >  6:165429/1‑105 (MQ=255)
                                                                                                                                       gagGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGt                                                      <  5:165429/105‑1 (MQ=255)
                                                                                                                                       gagGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGcc          >  2:5581/1‑149 (MQ=255)
                                                                                                                                       gagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGt                                 <  2:331782/126‑1 (MQ=255)
                                                                                                                                       gagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGt                                 >  1:331782/1‑126 (MQ=255)
                                                                                                                                        agGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGcc                                           >  4:232408/1‑115 (MQ=255)
                                                                                                                                        agGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGcc                                           <  3:232408/115‑1 (MQ=255)
                                                                                                                                          gCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGAGCCTCTGATCTGTAACAATACATACGTAAATTAGCCGc        >  8:150123/1‑148 (MQ=255)
                                                                                                                                               tCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACgg  <  2:44204/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATCCTGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGG  >  NC_002947/321594‑321885

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