Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,555,389 T→A 9.8% intergenic (+27/+47) 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,3890TA 9.8% 141.0 / 6.3 51intergenic (+27/+47)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base A (3/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTT  >  NC_002947/1555241‑1555536
                                                                                                                                                     |                                                                                                                                                   
tGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCTATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑t                                                                                                                                                     >  8:96708/1‑149 (MQ=255)
      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                                >  6:368947/1‑140 (MQ=255)
                                                                                                                              aGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCt                                <  5:368947/140‑1 (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)
                                                                                                                                                     |                                                                                                                                                   
TGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAA‑TGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTT  >  NC_002947/1555241‑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: ‑