Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,241,968:1 +T 8.9% intergenic (+27/+31) 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,9681.T 8.9% 163.5 / 11.1 45intergenic (+27/+31)grxD/bfr‑IImonothiol glutaredoxin/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base . (21/20);  new base T (2/2);  total (23/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.83e-01

AACTGTGGGTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGA  >  NC_002947/1241820‑1242114
                                                                                                                                                     |                                                                                                                                                  
aaCTGTGGGTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATGAAAAAACCcgcg                                                                                                                                                   <  2:1044375/150‑1 (MQ=255)
  cTGTGGGTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGAt                                                                                                                                                 >  1:621847/1‑150 (MQ=255)
             cGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCNTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATCGCGGGtttt                                                                                                                                      <  1:973242/150‑1 (MQ=255)
                 gAGCTGTTTGGCGACAGCGCCATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGTCTTCCAAGTCCAAGGCTTCTGAAGCTTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTa                                                                                                                                  <  1:61445/150‑1 (MQ=255)
                  aGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTAt                                                                                                                                 <  2:944719/150‑1 (MQ=255)
                   gCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATg                                                                                                                                >  1:82650/1‑150 (MQ=255)
                       gTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATCGCGGGTTTTTTTAt                                                                                                                                 <  2:434996/145‑1 (MQ=255)
                       gTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATCGCGGGTTTTTTTAt                                                                                                                                 >  1:434996/1‑145 (MQ=255)
                       gTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGc                                                                                                                            >  2:264760/1‑150 (MQ=255)
                            cggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATa                                                                                                                       <  2:621847/150‑1 (MQ=255)
                                aGCGACATCATCCAGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTa                                                                                                                   <  1:264760/150‑1 (MQ=255)
                                    aCATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTa                                                                                                                                  >  2:635177/1‑131 (MQ=255)
                                    aCATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTa                                                                                                                                  <  1:635177/131‑1 (MQ=255)
                                                    gTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACt                                                                                               >  2:965615/1‑150 (MQ=255)
                                                     ttCGAGAAGGGTGAGCTGCAAACCCTGATCAGGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTg                                                                                              >  1:26770/1‑150 (MQ=255)
                                                     ttCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTTCGGCAAGCATAAAAAAACCCGCGATCAGCGGGGTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATATGCGAcgg                                                                                              >  1:1053650/1‑148 (MQ=255)
                                                                   gCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                <  1:965615/150‑1 (MQ=255)
                                                                          aCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGAt                                                                         <  2:1145377/150‑1 (MQ=255)
                                                                            ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtctc                                                                                  >  2:313021/1‑139 (MQ=255)
                                                                            ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtctc                                                                                  <  1:313021/139‑1 (MQ=255)
                                                                             cTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                      >  2:190129/1‑150 (MQ=255)
                                                                                  cAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCa                                                                                            >  1:1114260/1‑123 (MQ=255)
                                                                                  cAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCa                                                                                            <  2:1114260/123‑1 (MQ=255)
                                                                                         gcgGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGccc                                                          >  1:254692/1‑150 (MQ=255)
                                                                                           ggCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCa                                                        >  1:793647/1‑150 (MQ=255)
                                                                                           ggCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCa                                                        <  2:82650/150‑1 (MQ=255)
                                                                                           ggCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCa                                                        >  2:1051533/1‑150 (MQ=255)
                                                                                            gCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAg                                                       >  2:265964/1‑150 (MQ=255)
                                                                                                ccaaggccaaggCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   <  1:1051533/150‑1 (MQ=255)
                                                                                                ccaaggccaaggCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   <  2:287579/150‑1 (MQ=255)
                                                                                                 caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCgg                                                                                                                                           >  1:923087/1‑61 (MQ=255)
                                                                                                 caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCgg                                                                                                                                           <  2:923087/61‑1 (MQ=255)
                                                                                                 caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg                                                  <  2:26770/150‑1 (MQ=255)
                                                                                                 caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCAATCAGCCCCAGTTggg                                                  >  1:805780/1‑150 (MQ=255)
                                                                                                     gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                              >  2:356681/1‑150 (MQ=255)
                                                                                                           gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   <  1:568827/139‑1 (MQ=255)
                                                                                                           gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                   >  2:568827/1‑139 (MQ=255)
                                                                                                              tCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTc                                     <  1:265964/150‑1 (MQ=255)
                                                                                                                 gAAGCCTGATCTTCCGGCAAGCATAAAAACACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                  <  2:1053650/150‑1 (MQ=255)
                                                                                                                 gAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                  <  1:687464/150‑1 (MQ=255)
                                                                                                                  aaGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTg                                                                     >  2:614720/1‑114 (MQ=255)
                                                                                                                  aaGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTg                                                                     <  1:614720/114‑1 (MQ=255)
                                                                                                                       tGATCTTCCTGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                      >  2:1282307/1‑108 (MQ=255)
                                                                                                                       tGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                      <  1:1282307/108‑1 (MQ=255)
                                                                                                                                            aaaaCCCGCTGATCGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCACCCCAAGATGAGAAAAAAGCCAGTCGATGTGCTCCTCCTCGGACTCGCGGATATc       >  2:140945/1‑150 (MQ=255)
                                                                                                                                                 ccGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGa  >  2:250499/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AACTGTGGGTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGA  >  NC_002947/1241820‑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: