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

ACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTT  >  NC_002947/1555247‑1555536
                                                                                                                                                |                                                                                                                                                  
aCTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCg                                                                                                                                               <  2:296749/149‑1 (MQ=255)
    cGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAC‑ATTGCCGGCg                                                                                                                                            <  5:352487/148‑1 (MQ=255)
     gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAC‑ATTGCCGGCGtt                                                                                                                                          <  6:278730/149‑1 (MQ=255)
     gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGtt                                                                                                                                          >  1:147499/1‑149 (MQ=255)
       aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGtttt                                                                                                                                        <  6:192112/149‑1 (MQ=255)
        aaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGttttt                                                                                                                                       >  1:41480/1‑149 (MQ=255)
           gcgcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTttgtt                                                                                                                                    <  2:41480/149‑1 (MQ=255)
                   gCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAg                                                                                                                            <  6:21245/149‑1 (MQ=255)
                      aGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCaa                                                                                                                         <  1:42307/149‑1 (MQ=255)
                            ggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCg                                                                                                                                               >  4:269500/1‑121 (MQ=255)
                            ggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCg                                                                                                                                               <  3:269500/121‑1 (MQ=255)
                                   gcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTtgttgt                                                                                                                                  <  4:266669/127‑1 (MQ=255)
                                   gcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTtgttgt                                                                                                                                  >  3:266669/1‑127 (MQ=255)
                                    ccGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCgg                                                                                                           <  1:236590/149‑1 (MQ=255)
                                       acgaGGCAGGCAGGCGCGTGTTGGA‑GCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa                                                                                                       >  2:237485/1‑149 (MQ=255)
                                            gcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAg                                                                                                   <  8:49454/149‑1 (MQ=255)
                                                    gcgcGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCTATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagcag                                                                                           <  7:96708/149‑1 (MQ=255)
                                                      gcgTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAggg                                                                                         >  6:214992/1‑149 (MQ=255)
                                                         tgtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTttgtt                                                                                                                                    >  6:238018/1‑103 (MQ=255)
                                                         tgtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTttgtt                                                                                                                                    <  5:238018/103‑1 (MQ=255)
                                                              gAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCa                                                                                  <  8:64869/148‑1 (MQ=255)
                                                                 gCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCtt                                                                              <  1:237485/149‑1 (MQ=255)
                                                                              aTGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcgc                                                                 >  5:29511/1‑149 (MQ=255)
                                                                                  tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCgccgcc                                                             >  3:206283/1‑149 (MQ=255)
                                                                                    gACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                            >  7:123693/1‑148 (MQ=255)
                                                                                            cAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGc                                                    <  8:123693/148‑1 (MQ=255)
                                                                                             aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCtt                                                  >  3:219319/1‑149 (MQ=255)
                                                                                                aaCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTc                                                >  4:77013/1‑148 (MQ=255)
                                                                                                   gCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCt                                                                               >  8:347558/1‑114 (MQ=255)
                                                                                                   gCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCt                                                                               <  7:347558/114‑1 (MQ=255)
                                                                                                          aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAAGTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTc                                                >  1:60459/1‑139 (MQ=255)
                                                                                                          aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTc                                                <  2:60459/138‑1 (MQ=255)
                                                                                                          aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAg                                     <  2:56572/149‑1 (MQ=255)
                                                                                                          aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCACAAGTGGCTTTCCACTCTTTCAg                                     >  8:335180/1‑149 (MQ=255)
                                                                                                                        aGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCt                                <  5:368947/140‑1 (MQ=255)
                                                                                                                        aGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCt                                >  6:368947/1‑140 (MQ=255)
                                                                                                                        aGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgc                       >  4:315225/1‑149 (MQ=255)
                                                                                                                         gATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgcc                      <  4:206283/149‑1 (MQ=255)
                                                                                                                           tAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccgc                    <  3:77013/149‑1 (MQ=255)
                                                                                                                             aaCAAAAACGCCGGCAC‑ATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCt                  >  7:256598/1‑149 (MQ=255)
                                                                                                                              aCAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAg                                     >  6:142009/1‑129 (MQ=255)
                                                                                                                              aCAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt                 >  2:16172/1‑149 (MQ=255)
                                                                                                                                aaaaaCGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg               >  3:285644/1‑149 (MQ=255)
                                                                                                                                 aaaaCGCCGGCAA‑TTTTCCGTCGTTTTTGTTGTGCGAAGAAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGTTTGCCCAGATTTCAAACACGTAGCCGCCCAAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGACGACTTTgg              >  4:183597/1‑149 (MQ=255)
                                                                                                                                  aaaCGCCGGCAC‑ATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGt             >  4:316311/1‑149 (MQ=255)
                                                                                                                                   aaCGCCGGCAC‑ATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGtt            >  4:119824/1‑149 (MQ=255)
                                                                                                                                    aCGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt           <  6:29511/149‑1 (MQ=255)
                                                                                                                                        cGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAt       <  4:285644/149‑1 (MQ=255)
                                                                                                                                         ggCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa        <  7:17069/147‑1 (MQ=255)
                                                                                                                                         ggCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa        >  8:17069/1‑147 (MQ=255)
                                                                                                                                             a‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATTCTGATAATCAGTCCAGCAGGTTTGCCCATCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCAATCTTTCAGGTTCTTTTGTTTGCCGCCTTTGGTTTCGATCACtt  <  3:183597/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTT  >  NC_002947/1555247‑1555536

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