Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,555,432 Δ1 bp 13.3% intergenic (+28/+46) sbcB → / ← PP_1366 exodeoxyribonuclease I/transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,555,4320G.13.3% 106.5 / 16.3 45intergenic (+28/+46)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base . (4/2);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

CGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATC  >  NC_002947/1555303‑1555574
                                                                                                                                  |                                                                                                                                              
cgcCAGCAAGCCTGGGAGGGCGCCTACGAGGCAGGTAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGc                                                                                                                            <  2:192971/150‑1 (MQ=255)
cgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCCGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACGTGTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGc                                                                                                                            <  1:73663/150‑1 (MQ=255)
           cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGACGCAAGAgg                                                                                                                 >  2:100159/1‑150 (MQ=255)
           cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg                                                                                                                 <  1:363789/150‑1 (MQ=255)
                        gacgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATa                                                                                                    >  1:830468/1‑150 (MQ=255)
                         acgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa                                                                                                   <  1:702716/150‑1 (MQ=255)
                            aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                                <  2:830468/150‑1 (MQ=255)
                              gcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGt                                                                                              >  1:975515/1‑150 (MQ=255)
                               caggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTc                                                                                             <  1:485733/150‑1 (MQ=255)
                                 ggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCACTCGGCTGACACCCAATAAACAAAAACGCCGGTA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCca                                                                                           <  1:134681/150‑1 (MQ=255)
                                    aggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagca                                                                                        <  2:975515/150‑1 (MQ=255)
                                                      tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACAT‑TGCCGGCGttt                                                                                                                                     >  2:406437/1‑87 (MQ=38)
                                                      tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑CATTGCCGGCGttt                                                                                                                                     <  1:406437/87‑1 (MQ=38)
                                                      tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCaa                                                                      <  1:1054519/150‑1 (MQ=255)
                                                      tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCaa                                                                      <  1:1054390/150‑1 (MQ=255)
                                                                    tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                        >  2:621453/1‑150 (MQ=255)
                                                                       aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCaa                                                     >  1:833204/1‑150 (MQ=255)
                                                                       aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCaa                                                     >  2:299231/1‑150 (MQ=255)
                                                                                  aaCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCa                                          <  1:212088/150‑1 (MQ=255)
                                                                                      cTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCActct                                      >  2:914783/1‑150 (MQ=255)
                                                                                         ttCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGCCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTc                                   <  1:131576/150‑1 (MQ=255)
                                                                                            aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                <  1:645416/150‑1 (MQ=255)
                                                                                            aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                >  1:403643/1‑150 (MQ=255)
                                                                                             aTCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAggg                               >  2:937422/1‑150 (MQ=255)
                                                                                             aTCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCAATCTTTCatgt                               >  2:630957/1‑147 (MQ=255)
                                                                                                   tGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGAGTTGACCAGCTTACAACCACTTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgt                         >  1:1098123/1‑150 (MQ=255)
                                                                                                            aTAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccgc                >  2:205847/1‑150 (MQ=255)
                                                                                                               aaCAAAAACGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt             >  1:348650/1‑150 (MQ=255)
                                                                                                                aCAAAAACGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt            >  2:784777/1‑150 (MQ=255)
                                                                                                                aCAAAAACGCCGGCA‑CATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt            >  2:787042/1‑150 (MQ=255)
                                                                                                                aCAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt            >  1:414793/1‑150 (MQ=255)
                                                                                                                 caacaaCCCCCGCA‑ACGCGCCGGCGTTTTTGTTGTCCGCAGCAAGAGGGCGTATGCGGATAATCAGGCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAACTGGCTTTCCTCTCATTCAGGGTCTTGTGGTTGCCGCCTTTg           <  2:414793/150‑1 (MQ=255)
                                                                                                                 cAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg           >  2:624401/1‑150 (MQ=255)
                                                                                                                  aaaaaCGCCGGCA‑CATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg          >  1:837055/1‑150 (MQ=255)
                                                                                                                   aaaaCGCCGGCA‑ATGTGACACGGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGt         <  1:899565/150‑1 (MQ=255)
                                                                                                                    aaaCGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTc                                                               <  1:359459/95‑1 (MQ=38)
                                                                                                                    aaaCGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTc                                                               >  2:359459/1‑95 (MQ=38)
                                                                                                                     aaCGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt       <  2:403643/150‑1 (MQ=255)
                                                                                                                      aCGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAg                                                                                               <  2:887196/61‑1 (MQ=25)
                                                                                                                      aCGCCGGCA‑CATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAg                                                                                               >  1:887196/1‑61 (MQ=25)
                                                                                                                       cGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCg     >  2:654041/1‑150 (MQ=255)
                                                                                                                        gCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa    >  2:605409/1‑150 (MQ=255)
                                                                                                                        gCCGGCA‑ATGTGCCGGCGTTTTTATTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGTTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa    <  1:860269/150‑1 (MQ=255)
                                                                                                                          cGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATc  >  1:14636/1‑150 (MQ=255)
                                                                                                                          cGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGAATGATGATATTCTGCCCAGCGGGGTTACCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTATGGTTGCCGCCTTTGGTTTCGATc  >  1:22324/1‑150 (MQ=255)
                                                                                                                                  |                                                                                                                                              
CGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATC  >  NC_002947/1555303‑1555574

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