Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,241,967 Δ1 bp 8.8% intergenic (+32/+26) grxD → / ← bfr‑II monothiol glutaredoxin/bacterioferritin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,241,9670A. 8.8% 176.2 / 11.0 57intergenic (+32/+26)grxD/bfr‑IImonothiol glutaredoxin/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base . (2/3);  total (26/31)
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.36e-01

GTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGA  >  NC_002947/1241822‑1242108
                                                                                                                                                  |                                                                                                                                             
gTAGCCGGTGAGCTGTTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCg                                                                                                                                            <  8:243338/149‑1 (MQ=255)
          aGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGttttttt                                                                                                                                   <  3:271283/148‑1 (MQ=255)
            cTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATg                                                                                                                                <  2:360979/149‑1 (MQ=255)
               gTTGGCGGCAGCGAGATCATGCTTGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGTTTTTTTATgc                                                                                                                               <  6:253464/147‑2 (MQ=39)
                          cGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTAc                                                                                                                  >  2:315643/1‑149 (MQ=255)
                          cGACATCATGCTGGAAAAGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTAc                                                                                                                  >  6:34959/1‑149 (MQ=255)
                           gACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCATGATCTTCCGGCAAGCATAAAAAAACACGC‑GATCAGCGGGTTTTTTTATGGGGAGAATAGATAc                                                                                                                  >  3:259365/1‑148 (MQ=255)
                                    cTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCa                                                                                                        >  6:199101/1‑149 (MQ=255)
                                     tGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATgg                                                                                                                               <  8:100626/125‑1 (MQ=255)
                                     tGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATgg                                                                                                                               >  7:100626/1‑125 (MQ=255)
                                     tGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCa                                                                                                        >  4:38958/1‑148 (MQ=255)
                                      ggAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATc                                                                                                      >  6:235778/1‑149 (MQ=255)
                                       gAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCAt                                                                                                       <  5:34959/147‑1 (MQ=255)
                                       gAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCt                                                                                                     >  3:249550/1‑149 (MQ=255)
                                           tGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGa                                                                                                 <  5:199101/149‑1 (MQ=255)
                                             ttCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACt                                                                                               <  5:235778/149‑1 (MQ=255)
                                                  ttgatgTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCACGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAg                                                                                           <  8:194585/143‑1 (MQ=255)
                                                  gaAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATggg                                                                                                                              <  6:329727/113‑1 (MQ=255)
                                                         gAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtct                                                                                   >  3:319864/1‑149 (MQ=255)
                                                               cAAACCCTGATCAAGGACGCTGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATa                                                                              <  6:310485/148‑1 (MQ=255)
                                                               cAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGTCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATAc                                                                             <  1:114993/149‑1 (MQ=255)
                                                                 aaCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATAcc                                                                            >  3:342848/1‑148 (MQ=255)
                                                                    ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTa                                                                                                                   <  7:66462/106‑1 (MQ=255)
                                                                    ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTa                                                                                                                   >  8:66462/1‑106 (MQ=255)
                                                                            aGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTCCAGGTAGTTCTCGATACCGATCTTGTCGAt                                                                >  2:243714/1‑149 (MQ=255)
                                                                                cgcgGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTAGTCGATCAGc                                                            <  1:315643/149‑1 (MQ=255)
                                                                                  cgGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGccc                                                          <  4:342848/149‑1 (MQ=255)
                                                                                      tGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGTGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGt                                                                    >  5:109994/1‑135 (MQ=255)
                                                                                      tGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGTGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGt                                                                    <  6:109994/135‑1 (MQ=255)
                                                                                        ccaaggccaaggCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTg                                                    >  6:301873/1‑149 (MQ=255)
                                                                                         caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   >  1:140698/1‑149 (MQ=255)
                                                                                         caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   >  6:113012/1‑149 (MQ=255)
                                                                                          aaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg                                                  <  3:38958/149‑1 (MQ=255)
                                                                                            ggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGtt                                                <  1:300344/149‑1 (MQ=255)
                                                                                              ccaaggcTTCTGAAGCCTGATCTTCAGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                              <  4:52773/149‑1 (MQ=255)
                                                                                               caaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                              >  5:87517/1‑148 (MQ=255)
                                                                                                aaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                              >  6:191171/1‑147 (MQ=255)
                                                                                                  ggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccagc                                          <  5:113012/149‑1 (MQ=255)
                                                                                                   gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACTCGC‑GATAAGCGGTTTTTTTTATTGTGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTTGGTTTccagc                                          <  4:259365/148‑1 (MQ=255)
                                                                                                     ttCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccagccag                                       >  3:243282/1‑149 (MQ=255)
                                                                                                          aaGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                  >  4:166137/1‑149 (MQ=255)
                                                                                                                 aTCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGctc                            <  4:249550/148‑1 (MQ=255)
                                                                                                                 aTCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGctc                            <  6:87517/148‑1 (MQ=255)
                                                                                                                       cGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTc                      <  3:166137/148‑1 (MQ=255)
                                                                                                                        ggCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCg                     >  6:88382/1‑148 (MQ=255)
                                                                                                                              cATAAAAAAACCCGC‑GATCAGCGGGTTTTTATATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGa              <  8:222948/149‑1 (MQ=255)
                                                                                                                               aTAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAg             <  4:319864/149‑1 (MQ=255)
                                                                                                                                  aaaaaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCAATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGa           >  2:211487/1‑148 (MQ=255)
                                                                                                                                   aaaaaCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   <  7:261246/107‑1 (MQ=255)
                                                                                                                                   aaaaaCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   >  8:261246/1‑107 (MQ=255)
                                                                                                                                     aaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTc                      >  2:61462/1‑134 (MQ=255)
                                                                                                                                     aaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTc                      <  1:61462/134‑1 (MQ=255)
                                                                                                                                         ccGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                      >  6:362013/1‑82 (MQ=255)
                                                                                                                                         ccGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                      <  5:362013/82‑1 (MQ=255)
                                                                                                                                         ccGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCg   <  5:88382/149‑1 (MQ=255)
                                                                                                                                         ccGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTACTCGGACTCGAGGATATCt      >  7:54328/1‑146 (MQ=255)
                                                                                                                                          cgc‑gATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGACCATCTGCGACTGAAGGTAGTTCTCGATACGGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGa  <  1:341588/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGA  >  NC_002947/1241822‑1242108

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