Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,555,390 G→T 11.8% 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,3900GT11.8% 137.6 / 10.0 51intergenic (+28/+46)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base T (4/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

TTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAT  >  NC_002947/1555249‑1555531
                                                                                                                                             |                                                                                                                                             
ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGc                                                                                                                                        >  2:129057/1‑149 (MQ=255)
   gAAAAGGCGCGCCAGCAAGCCTGGGAGGGTGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtt                                                                                                                                     >  4:182268/1‑149 (MQ=255)
     aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCGTGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGtttt                                                                                                                                   <  5:357478/149‑1 (MQ=255)
       aGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttg                                                                                                                                 <  1:129057/149‑1 (MQ=255)
       aGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttg                                                                                                                                 <  2:269047/149‑1 (MQ=255)
        cacgcgcctgcctgcctcgtcggcgccctcccaggcttgctggcgcgTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGcc                                                                                                                                           <  1:185474‑M2/91‑1 (MQ=255)
        cacgcgcctgcctgcctcgtcggcgccctcccaggcttgctggcgcgTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGcc                                                                                                                                           >  2:185474‑M2/48‑138 (MQ=255)
           GCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTtgttgt                                                                                                                             <  1:227922‑M1/149‑1 (MQ=255)
           GCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTtgttgt                                                                                                                             <  2:127745‑M1/149‑1 (MQ=255)
             GCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGc                                                                                                                           <  7:65692‑M1/149‑1 (MQ=255)
               CAGCAAGCCTGGGAGGGCTCCGACGAGGCAGGCAGGCGCGTGTTGTAGGCCTGGCAACTGCATGCGCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTTTTGTGCGa                                                                                                                         >  3:352167‑M1/1‑149 (MQ=255)
                aGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCg                                                                                                                          >  2:48718/1‑147 (MQ=255)
                   aaGCCTGGGAGGGTGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCa                                                                                                                     <  3:182268/149‑1 (MQ=255)
                   aaGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGc                                                                                                                      >  6:74704/1‑148 (MQ=255)
                    aGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCaa                                                                                                                    <  4:239145/149‑1 (MQ=255)
                      ccTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAaga                                                                                                                  <  3:267900/149‑1 (MQ=255)
                          ggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg                                                                                                                <  2:201949/147‑1 (MQ=255)
                                       gaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAAt                                                                                                 >  1:113357/1‑149 (MQ=255)
                                         ggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAAt                                                                                                 >  3:276310/1‑147 (MQ=255)
                                          gctggatggcgcgTGTTGTTGGCCTGGCAACTCCATGCTCTACATTTGCAGTAACGATATTCAATCGGCTTATGCCATATAAAAAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                               <  2:345545‑M2/135‑1 (MQ=255)
                                                  gcgcgTGTTGTAGGCCTGGAAACTGCATGCTCGACAATTGCAGTAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGc                                                                                                                           <  6:243655‑M2/107‑1 (MQ=255)
                                                  gcgcgTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGc                                                                                                                           >  5:243655‑M2/6‑112 (MQ=255)
                                                    gcgTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAggg                                                                                    <  8:65692‑M2/146‑1 (MQ=255)
                                                        gtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAtgt                                                                                                                                              <  1:127935/87‑1 (MQ=255)
                                                        gtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAtgt                                                                                                                                              >  2:127935/1‑87 (MQ=255)
                                                           ggAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                        <  1:150144/122‑1 (MQ=255)
                                                           ggAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                        >  2:150144/1‑122 (MQ=255)
                                                             aGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGc                                                                           <  8:45826/149‑1 (MQ=255)
                                                              ggCCTGGCAACTGCATGCTCGACAATTGCGGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCt                                                                          >  3:344788/1‑149 (MQ=255)
                                                                  tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTc                                                                       >  4:60842/1‑148 (MQ=255)
                                                                    gCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCaa                                                                     <  3:60842/148‑1 (MQ=255)
                                                                        cTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCAc                                                                 >  4:229337/1‑148 (MQ=255)
                                                                                   aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGt                                                                                             <  4:89452/109‑1 (MQ=255)
                                                                                   aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGt                                                                                             >  3:89452/1‑109 (MQ=255)
                                                                                       ttGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTg                                                 >  7:222974/1‑149 (MQ=255)
                                                                                         gCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTgg                                                <  4:344788/148‑1 (MQ=255)
                                                                                            ggAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                            <  4:276310/149‑1 (MQ=255)
                                                                                                     ttCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCtt                                    >  6:96089/1‑148 (MQ=255)
                                                                                                        aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa                                                                                                  <  5:38514/83‑1 (MQ=255)
                                                                                                        aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa                                                                                                  >  6:38514/1‑83 (MQ=255)
                                                                                                                 aTGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtg                         <  2:227922/147‑1 (MQ=255)
                                                                                                                        cacaaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccg                >  3:107758/4‑149 (MQ=255)
                                                                                                                         acaaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccgc               >  8:214764/3‑149 (MQ=255)
                                                                                                                            aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt            >  8:255329/1‑149 (MQ=255)
                                                                                                                              aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg          <  1:48718/149‑1 (MQ=255)
                                                                                                                               aaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg         <  7:255329/149‑1 (MQ=255)
                                                                                                                                 aaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGccc                                                                              >  7:263243/1‑78 (MQ=38)
                                                                                                                                 aaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGccc                                                                              <  8:263243/78‑1 (MQ=38)
                                                                                                                                  aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt      >  5:259865/1‑149 (MQ=255)
                                                                                                                                  aCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt      >  6:12268/1‑149 (MQ=255)
                                                                                                                                     ccGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa   >  5:322934/1‑149 (MQ=255)
                                                                                                                                      cGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGGAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAt  <  3:229337/149‑1 (MQ=255)
                                                                                                                                      cGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAt  <  5:12268/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
TTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAT  >  NC_002947/1555249‑1555531

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