Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,569,457 A→T 10.9% intergenic (+16/+22) PP_3150 → / ← sad‑II hypothetical protein/NAD+‑dependent succinate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,569,4570AT10.9% 124.4 / 5.7 46intergenic (+16/+22)PP_3150/sad‑IIhypothetical protein/NAD+‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base T (3/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

AGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAG  >  NC_002947/3569311‑3569602
                                                                                                                                                  |                                                                                                                                                 
aGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCa                                                                                                                                                   >  5:96852/1‑147 (MQ=255)
aGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCa                                                                                                                                                   >  2:26077/1‑147 (MQ=255)
       gggTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAGGAACCTGGGGCGGCAAAGCAGccc                                                                                                                                          >  1:53103/1‑149 (MQ=255)
       gggTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGccc                                                                                                                                          >  1:306229/1‑149 (MQ=255)
          tAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAg                                                                                                                                       <  1:249500/149‑1 (MQ=255)
            aCCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCTGCTTTGCGGCCCCAgct                                                                                                                                     <  1:11552/149‑3 (MQ=255)
              cGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGATTCATGGACTGGGGCAGGGACCCGGTGCGGGTTTGATCCGCATGGGGCCGAAATGCAGCCCCTggtac                                                                                                                                   >  1:6576/1‑146 (MQ=255)
                cAATGAACGGGCGCGTATGCTGGCTGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCAt                                                                                                                                 <  4:277863/149‑1 (MQ=255)
                cAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCAt                                                                                                                                 >  2:302374/1‑149 (MQ=255)
                cAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCAt                                                                                                                                 >  7:226182/1‑149 (MQ=255)
                 aaTGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCAt                                                                                                                                 >  7:193982/1‑148 (MQ=255)
                  aTGAACGGGCGCGTATGCTGGCGGGGTTTTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGa                                                                                                                               <  4:224494/149‑1 (MQ=255)
                   tGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGaa                                                                                                                              <  5:110336/149‑1 (MQ=255)
                                aTGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTc                                                                                                                  >  3:149812/1‑148 (MQ=255)
                                 tGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCa                                                                                                                <  2:100616/149‑1 (MQ=255)
                                    tGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGcc                                                                                                             <  2:306229/149‑1 (MQ=255)
                                         gggTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACg                                                                                                        >  1:115554/1‑149 (MQ=255)
                                               tataCCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTgg                                                                                                  <  6:96852/149‑1 (MQ=255)
                                                ataCCGCCTGTTTGCGGGGGCGGGATTAATTTTTCCTTCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTgg                                                                                                  <  2:6576/148‑1 (MQ=255)
                                                    cGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTg                                                                                            >  8:94592/1‑149 (MQ=255)
                                                     gCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTg                                                                                            >  8:303736/1‑149 (MQ=255)
                                                      ccTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGt                                                                                          <  7:94592/149‑1 (MQ=255)
                                                       cTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATTAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGt                                                                                          <  5:90838/149‑1 (MQ=255)
                                                            cgcgGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAAc                                                                                      <  2:115554/148‑1 (MQ=255)
                                                               gggCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTcgc                                                                                  <  8:226182/149‑1 (MQ=255)
                                                                        aTGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCa                                                                          >  2:53868/1‑148 (MQ=255)
                                                                            aTTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTg                                                                     >  2:59622/1‑149 (MQ=255)
                                                                             ttGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGc                                                                    <  7:303736/149‑1 (MQ=255)
                                                                              tGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGc                                                                    >  4:274694/1‑148 (MQ=255)
                                                                                aTCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGa                                                                  <  4:149812/148‑1 (MQ=255)
                                                                                 tCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATAAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGa                                                                  >  1:318815/1‑147 (MQ=255)
                                                                                    tGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCt                                                             <  7:171571/149‑1 (MQ=255)
                                                                                               aTGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGAAAGCTCACGGCCGaaa                                                  <  8:193982/149‑1 (MQ=255)
                                                                                                           gTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCttcttc                                      <  1:302374/149‑1 (MQ=255)
                                                                                                                 gACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGcca                                >  3:62281/1‑149 (MQ=255)
                                                                                                                  aCCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGcca                                <  1:59622/148‑1 (MQ=255)
                                                                                                                             gATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCa                                     >  4:208994/1‑132 (MQ=255)
                                                                                                                             gATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCa                                     <  3:208994/132‑1 (MQ=255)
                                                                                                                              aTTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCAcgc                   >  3:200865/1‑149 (MQ=255)
                                                                                                                                   agcagcTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCggggg              >  6:230307/1‑149 (MQ=255)
                                                                                                                                     acccTGGGGCTGCTTTGCGGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCAc                                              <  6:294288/112‑1 (MQ=255)
                                                                                                                                     acccTGGGGCTGCTTTGCGGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCAc                                              >  5:294288/4‑115 (MQ=255)
                                                                                                                                      cccTGGGGCTGCTTTGCGGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTcc           >  1:98378/3‑149 (MQ=255)
                                                                                                                                       gcTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTcc           >  8:245303/1‑148 (MQ=255)
                                                                                                                                       ccTGGGGCTGCTTTGCGGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTcc           >  2:63019/2‑148 (MQ=255)
                                                                                                                                           gggCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTGACGCCACCAAACGCCACGCGGGGGTCCGaa        <  1:26077/147‑1 (MQ=255)
                                                                                                                                            ggCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAgcg     <  1:53868/149‑1 (MQ=255)
                                                                                                                                               cGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAg  >  7:151701/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAG  >  NC_002947/3569311‑3569602

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: