Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 322,140 T→C 12.5% intergenic (+427/+31) 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_002947322,1400TC12.5% 138.0 / 9.9 56intergenic (+427/+31)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base C (3/4);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

GGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCC  >  NC_002947/321993‑322282
                                                                                                                                                   |                                                                                                                                              
ggAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTc                                                                                                                                               >  4:68316/1‑149 (MQ=255)
  aaGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgc                                                                                                                                             <  6:315364/149‑1 (MQ=255)
  aaGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgc                                                                                                                                             <  2:179833/149‑1 (MQ=255)
     cATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                          <  8:297920/149‑1 (MQ=255)
      aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgcg                                                                                                                                            >  6:207839/1‑146 (MQ=255)
       tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                        <  7:295468/149‑1 (MQ=255)
       tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                        <  5:180101/149‑1 (MQ=255)
       tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGCCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                        <  6:153566/149‑1 (MQ=255)
         aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAgg                                                                                                                                      <  6:317783/149‑1 (MQ=255)
           ggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCa                                                                                                                                    <  6:292916/149‑1 (MQ=255)
                    gaccgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAgg                                                                                                                            >  1:330724/1‑148 (MQ=255)
                         accCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACtt                                                                                                                      <  5:207839/149‑1 (MQ=255)
                          cccACGACGTCTGCTTGCCTGGCGATTCCCCCAAGGGTTTATCAGCGAAGCAGAGATGCAGGCGGTGATGGCGGTTTGGCACTTATCGATTTCACCTGTACCTGCCTCTTCGCGGGCATGCTCTCGAATAGGCAGGTACAGGCAACtt                                                                                                                      <  4:144745/148‑1 (MQ=255)
                            cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCTAAGAGGCAGGTACAGGCAACTTaaa                                                                                                                   <  3:64330/149‑1 (MQ=255)
                            cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaa                                                                                                                    >  4:207933/1‑148 (MQ=255)
                             acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                   >  1:153280/1‑132 (MQ=255)
                             acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                   <  2:153280/132‑1 (MQ=255)
                                  gTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGTTCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                              >  7:51014/1‑148 (MQ=255)
                                      gCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTAAAGGCAACTTAAAATAAGAAACt                                                                                                         >  7:177149/1‑149 (MQ=255)
                                          gCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAg                                                                                                     >  4:17047/1‑149 (MQ=255)
                                               gCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAcc                                                                                                >  6:208749/1‑149 (MQ=255)
                                                  aTACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaaa                                                                                                                  >  3:85429/1‑128 (MQ=255)
                                                  aTACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaaa                                                                                                                  <  4:85429/128‑1 (MQ=255)
                                                       cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacca                                                                                        >  6:116719/1‑149 (MQ=255)
                                                         cAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACACGCAACTTAAAATCAGAAACTGTAGCTACCCGTcaccac                                                                                       >  5:251111/1‑148 (MQ=255)
                                                          aaGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACaa                                                                                     <  3:17047/149‑1 (MQ=255)
                                                            gggTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAgg                                                                                   <  3:207933/149‑1 (MQ=255)
                                                                  aTCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcg                                                                              <  3:68316/148‑1 (MQ=255)
                                                                     aGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                      <  6:76668/121‑1 (MQ=255)
                                                                     aGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                      >  5:76668/1‑121 (MQ=255)
                                                                      gCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCaa                                                                                                                         >  2:129253/1‑101 (MQ=255)
                                                                      gCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCaa                                                                                                                         <  1:129253/101‑1 (MQ=255)
                                                                       cGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATc                                                                        >  5:23696/1‑149 (MQ=255)
                                                                               agaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCt                                                                <  5:208749/149‑1 (MQ=255)
                                                                                    cAGGCGCTGATGGCGGTGTGGCACTGTTCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGa                                                             <  8:51014/147‑1 (MQ=255)
                                                                                       gcgcTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTgc                                                        <  8:258381/149‑1 (MQ=255)
                                                                                         gcTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTgcgc                                                      <  2:330724/149‑1 (MQ=255)
                                                                                           tGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                      <  2:204520/99‑1 (MQ=255)
                                                                                           tGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                      >  1:204520/1‑99 (MQ=255)
                                                                                                                cGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgtag                                >  5:123770/1‑148 (MQ=255)
                                                                                                                         cTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                             >  3:246645/1‑94 (MQ=255)
                                                                                                                         cTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                             <  4:246645/94‑1 (MQ=255)
                                                                                                                            tACCTTCCTGTTCGCGGTCCTGCTCTCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAATCTGTAGCTACCCGTCACCACAAGGCTACTGGGATCTCCCACCTGGAGCTGCGCCGCACTGGTCTCCTAGATGTAGTAGGTATTGTCa                    <  6:314233/148‑1 (MQ=255)
                                                                                                                            tACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAg                   <  6:123770/149‑1 (MQ=255)
                                                                                                                               cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa                <  6:23696/149‑1 (MQ=255)
                                                                                                                               cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa                >  8:274182/1‑149 (MQ=255)
                                                                                                                               cTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa                >  2:5297/1‑149 (MQ=255)
                                                                                                                                 gCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATg              >  3:137470/1‑149 (MQ=255)
                                                                                                                                  ccTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATg              <  5:116719/148‑1 (MQ=255)
                                                                                                                                  ccTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGt             >  8:254758/1‑149 (MQ=255)
                                                                                                                                  ccTCTTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATg              <  1:5297/148‑1 (MQ=255)
                                                                                                                                    tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg           >  2:306711/1‑149 (MQ=255)
                                                                                                                                     ctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc          >  5:273065/1‑149 (MQ=255)
                                                                                                                                      ttCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCt         >  6:246796/1‑149 (MQ=255)
                                                                                                                                      ttCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACACGTAACCACAAGGCTACGCGCATCGCCCACATGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTTTCAGCAATGTTGCt         >  6:66761/1‑149 (MQ=255)
                                                                                                                                             gCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATgcc  >  2:161081/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCC  >  NC_002947/321993‑322282

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: ‑