Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,241,973 Δ1 bp 9.1% 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,9730A. 9.1% 133.0 / 7.8 44intergenic (+32/+26)grxD/bfr‑IImonothiol glutaredoxin/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base . (2/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGA  >  NC_002947/1241833‑1242114
                                                                                                                                             |                                                                                                                                             
cGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCNTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGtttt                                                                                                                                      <  1:973242/150‑1 (MQ=255)
    gAGCTGTTTGGCGACAGCGCCATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGTCTTCCAAGTCCAAGGCTTCTGAAGCTTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTa                                                                                                                                  <  1:61445/150‑1 (MQ=255)
     aGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTAt                                                                                                                                 <  2:944719/150‑1 (MQ=255)
      gCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATg                                                                                                                                >  1:82650/1‑150 (MQ=255)
          gTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGTTTTTTTAt                                                                                                                                 <  2:434996/145‑1 (MQ=255)
          gTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGTTTTTTTAt                                                                                                                                 >  1:434996/1‑145 (MQ=255)
          gTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGc                                                                                                                            >  2:264760/1‑150 (MQ=255)
               cggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATa                                                                                                                       <  2:621847/150‑1 (MQ=255)
                   aGCGACATCATCCAGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTa                                                                                                                   <  1:264760/150‑1 (MQ=255)
                       aCATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTa                                                                                                                                  >  2:635177/1‑131 (MQ=255)
                       aCATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTa                                                                                                                                  <  1:635177/131‑1 (MQ=255)
                                       gTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACt                                                                                               >  2:965615/1‑150 (MQ=255)
                                        ttCGAGAAGGGTGAGCTGCAAACCCTGATCAGGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTg                                                                                              >  1:26770/1‑150 (MQ=255)
                                        ttCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTTCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGGTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATATGCGAcgg                                                                                              >  1:1053650/1‑148 (MQ=255)
                                                      gCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                <  1:965615/150‑1 (MQ=255)
                                                             aCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGAt                                                                         <  2:1145377/150‑1 (MQ=255)
                                                               ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtctc                                                                                  >  2:313021/1‑139 (MQ=255)
                                                               ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtctc                                                                                  <  1:313021/139‑1 (MQ=255)
                                                                cTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                      >  2:190129/1‑150 (MQ=255)
                                                                     cAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCa                                                                                            >  1:1114260/1‑123 (MQ=255)
                                                                     cAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCa                                                                                            <  2:1114260/123‑1 (MQ=255)
                                                                            gcgGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGccc                                                          >  1:254692/1‑150 (MQ=255)
                                                                              ggCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCa                                                        >  1:793647/1‑150 (MQ=255)
                                                                              ggCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCa                                                        <  2:82650/150‑1 (MQ=255)
                                                                              ggCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCa                                                        >  2:1051533/1‑150 (MQ=255)
                                                                               gCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAg                                                       >  2:265964/1‑150 (MQ=255)
                                                                                   ccaaggccaaggCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   <  1:1051533/150‑1 (MQ=255)
                                                                                   ccaaggccaaggCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   <  2:287579/150‑1 (MQ=255)
                                                                                    caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCgg                                                                                                                                           <  2:923087/61‑1 (MQ=255)
                                                                                    caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCgg                                                                                                                                           >  1:923087/1‑61 (MQ=255)
                                                                                    caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg                                                  <  2:26770/150‑1 (MQ=255)
                                                                                    caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCAATCAGCCCCAGTTggg                                                  >  1:805780/1‑150 (MQ=255)
                                                                                        gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                              >  2:356681/1‑150 (MQ=255)
                                                                                              gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   <  1:568827/139‑1 (MQ=255)
                                                                                              gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   >  2:568827/1‑139 (MQ=255)
                                                                                                 tCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTc                                     <  1:265964/150‑1 (MQ=255)
                                                                                                    gAAGCCTGATCTTCCGGCAAGCATAAAAACACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                  <  2:1053650/150‑1 (MQ=255)
                                                                                                    gAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                  <  1:687464/150‑1 (MQ=255)
                                                                                                     aaGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTg                                                                     >  2:614720/1‑114 (MQ=255)
                                                                                                     aaGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTg                                                                     <  1:614720/114‑1 (MQ=255)
                                                                                                          tGATCTTCCTGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                      >  2:1282307/1‑108 (MQ=255)
                                                                                                          tGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                      <  1:1282307/108‑1 (MQ=255)
                                                                                                                               aaaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCACCCCAAGATGAGAAAAAAGCCAGTCGATGTGCTCCTCCTCGGACTCGCGGATATc       >  2:140945/1‑150 (MQ=255)
                                                                                                                                    ccGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGa  >  2:250499/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                             
CGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGA  >  NC_002947/1241833‑1242114

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