Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,570,114 C→T 9.7% A251T (GCC→ACC)  sad‑II ← NAD+‑dependent succinate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,570,1140CT 9.7% 112.9 / 2.4 41A251T (GCC→ACC) sad‑IINAD+‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (21/16);  new base T (2/2);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

AGCGGGTTGCCGACTTTCAACTGGCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCGGTCAGGGTCACGGCGGCGATGC  >  NC_002947/3569966‑3570257
                                                                                                                                                    |                                                                                                                                               
aGCGGGTTGCCGACTTTCAACTGGCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGAcggc                                                                                                                                                 <  7:88068/149‑1 (MQ=255)
 gCGGGTTGCCGACTTTCAACTGGCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTAGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGAcggcg                                                                                                                                                >  5:215299/1‑149 (MQ=255)
        gCCGACTTTCAACTGGCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATcc                                                                                                                                          <  2:273622/148‑1 (MQ=255)
         ccGACTTTCAACTGGCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAg                                                                                                                                        >  8:263928/1‑149 (MQ=255)
             cTTTCAACTGTCGAGTGGCGTCAACTAACTTGTGGGTGAACGCGTCGACTATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCTGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGTCGGCATCCAGGTCg                                                                                                                                    <  6:215299/149‑1 (MQ=255)
             cTTTCAACTGGCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCg                                                                                                                                    <  2:69898/149‑1 (MQ=255)
                tCAACTGGCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCgg                                                                                                                                   <  7:263928/147‑1 (MQ=255)
                   aCTGGCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCg                                                                                                                              <  1:6927/149‑1 (MQ=255)
                       gCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGc                                                                                                                            >  8:260637/1‑147 (MQ=255)
                             ggCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGa                                                                                                                    >  8:173068/1‑149 (MQ=255)
                                  cAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGa                                                                                                                 <  7:260637/147‑1 (MQ=255)
                                   aaCGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAgg                                                                                                              >  2:193875/1‑149 (MQ=255)
                                      gAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGa                                                                                                            <  6:137972/148‑1 (MQ=255)
                                      gAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGa                                                                                                            >  5:137972/1‑148 (MQ=255)
                                                    aCGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGAcggcg                                                                                                                                                >  5:276954/1‑98 (MQ=255)
                                                    aCGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGAcggcg                                                                                                                                                <  6:276954/98‑1 (MQ=255)
                                                        gTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTcc                                                                                          >  6:190516/1‑148 (MQ=255)
                                                         tCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCaa                                                                                        <  1:193875/149‑1 (MQ=255)
                                                                 gCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAacac                                                                                     >  4:84986/1‑144 (MQ=255)
                                                                 gCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACCGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAacac                                                                                     <  3:84986/144‑1 (MQ=255)
                                                                    ggCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGc                                                                                                     >  6:150721/1‑125 (MQ=255)
                                                                    ggCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGc                                                                                                     <  5:150721/125‑1 (MQ=255)
                                                                     gCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCActtctt                                                                            <  5:190516/149‑1 (MQ=255)
                                                                     gCTTCCACGATCAGGCGCTTGGCGGCGGCGCACAACTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGACTGGACGGCGGAATACAGGTCGGCATCGGCCAGCACGATGAAGTGATCGGAGCCTCCCAGCTCCAACAAGCActtct                                                                             >  3:208724/1‑148 (MQ=255)
                                                                       ttCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTccc                                                                                                >  1:173352/1‑127 (MQ=255)
                                                                       ttCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTccc                                                                                                <  2:173352/127‑1 (MQ=255)
                                                                           aCGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGTGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAgcgc                                                                      >  5:131122/1‑149 (MQ=255)
                                                                                 aGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGcc                                                                <  6:313239/149‑1 (MQ=255)
                                                                                 aGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGcc                                                                >  5:240815/1‑149 (MQ=255)
                                                                                  ggCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCg                                                               <  3:230230/149‑1 (MQ=255)
                                                                                                  gcACACCTGGCCGGTGTTCTGGTAGTGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCg                                               >  1:5244/1‑149 (MQ=255)
                                                                                                  gaacacCGGCCAGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGAt                                                                                                                   >  7:195262/3‑81 (MQ=255)
                                                                                                  gaacacCGGCCAGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGAt                                                                                                                   <  8:195262/79‑1 (MQ=255)
                                                                                                        cTGGCCGGTGTTCTGGTAGCGCCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGAACGATGAAGGGATCGGAGCCTCACAGCTACAACACGCACTTCTTCAGCGCGGAGACGGCCTGTGCTCCAATCGCCATGc                                         >  2:167136/1‑149 (MQ=255)
                                                                                                               gtgtTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGgcacg                                  >  8:190842/1‑149 (MQ=255)
                                                                                                                  ttCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGgcacgcac                               >  2:324974/1‑149 (MQ=255)
                                                                                                                       gTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGc                          >  8:184394/1‑149 (MQ=255)
                                                                                                                          gCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCgg                       >  7:263539/1‑149 (MQ=255)
                                                                                                                               cTATCACCGCGGCCTGGACGGTGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCGGTCAgg                  <  6:131122/149‑1 (MQ=255)
                                                                                                                               cTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACTATGAAGGGATCGGAGCCTCACAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCGGTCAgg                  <  4:208724/149‑1 (MQ=255)
                                                                                                                                     ccGCGGCCTGGACGGTGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCGGTCAGGGt                >  4:257400/1‑145 (MQ=255)
                                                                                                                                     ccGCGGCCTGGACGGTGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCGGTCAGGGt                <  3:257400/145‑1 (MQ=255)
                                                                                                                                     ccGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCGGTCAGGGTCAcg            >  8:316715/1‑149 (MQ=255)
                                                                                                                                              ggACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCGGTCAGGGTCACGGCGGCGATg   >  1:107253/1‑149 (MQ=255)
                                                                                                                                               gACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCGGTCAGGGTCACGGCGGCGATgc  >  1:160056/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AGCGGGTTGCCGACTTTCAACTGGCGAGTGGCTTCAACGAACTTGTGGGTGAACGCGTCGACGATGCTGGCTTCCACGATCAGGCGCTTGGCGGCGGCGCACACCTGGCCGGTGTTCTGGTAGCGGCCTATCACCGCGGCCTGGACGGCGGCATCCAGGTCGGCATCGGCCAGCACGATGAAGGGATCGGAGCCTCCCAGCTCCAACACGCACTTCTTCAGCGCGGCGCCGGCCTGTGCACCAATCGCCATGCCGGCACGCACGCTGCCGGTCAGGGTCACGGCGGCGATGC  >  NC_002947/3569966‑3570257

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