Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,569,459 A→T 12.9% intergenic (+18/+20) 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,4590AT12.9% 117.1 / 7.0 46intergenic (+18/+20)PP_3150/sad‑IIhypothetical protein/NAD+‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base T (4/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAG  >  NC_002947/3569318‑3569602
                                                                                                                                             |                                                                                                                                               
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                                                                                                                                 >  7:226182/1‑149 (MQ=255)
         cAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCAt                                                                                                                                 >  2:302374/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)
                                             cGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCC‑GGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTg                                                                                            >  8:94592/1‑149 (MQ=255)
                                              gCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTg                                                                                            >  8:303736/1‑149 (MQ=255)
                                               ccTGTTCGCGGGCGCGGCATGAATTGATCGTGCC‑GGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGt                                                                                          <  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)
                                                                                                                                             |                                                                                                                                               
GGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAG  >  NC_002947/3569318‑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: ‑