Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,570,596 G→A 8.7% H104H (CAC→CAT)  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,5960GA 8.7% 125.3 / 3.1 46H104H (CAC→CAT) sad‑IINAD+‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (14/28);  new base A (2/2);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TACCCGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCGGCCAGGGCCAGCAGGTAC  >  NC_002947/3570448‑3570740
                                                                                                                                                    |                                                                                                                                                
tACCCGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGt                                                                                                                                                 <  2:487490/150‑1 (MQ=255)
     gCCAGAATCGCCGGCACGGCGCCGCGCAGCACGTGCCAGACCGGGAAGTTCCCCGGCATCCCGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCg                                                                                                                                            <  2:624757/150‑1 (MQ=255)
     gCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCg                                                                                                                                            >  1:739650/1‑150 (MQ=255)
     gCCAGAATCGCCGGCACGCCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCg                                                                                                                                            <  1:440671/150‑1 (MQ=255)
       cAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGAGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGc                                                                                                                                          >  1:40100/1‑150 (MQ=255)
          aaTCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTa                                                                                                                                       <  1:980054/150‑1 (MQ=255)
             cGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCa                                                                                                                                    <  2:285489/150‑1 (MQ=255)
              gCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAg                                                                                                                                   >  2:708691/1‑150 (MQ=255)
                cGGCACGGCGCCGCGCAGCACCTGCCAGACCGAGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCg                                                                                                                                 <  2:40100/150‑1 (MQ=255)
                     cGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCaaa                                                                                                                            <  2:799720/150‑1 (MQ=255)
                      ggCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAg                                                                                                                           <  2:710233/150‑1 (MQ=255)
                                  caCCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACtt                                                                                                               >  2:949774/1‑150 (MQ=255)
                                      tGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTg                                                                                                           <  2:584785/150‑1 (MQ=255)
                                       gCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGa                                                                                                          <  1:444525/150‑1 (MQ=255)
                                        ccAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGAc                                                                                                         <  1:487076/150‑1 (MQ=255)
                                           gACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGGCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACttt                                                                                                      >  2:1227048/1‑149 (MQ=255)
                                              cGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTac                                                                                                                     >  1:242301/1‑132 (MQ=255)
                                              cGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTac                                                                                                                     <  2:242301/132‑1 (MQ=255)
                                              cGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGa                                                                                                          <  2:753187/143‑1 (MQ=255)
                                              cGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGa                                                                                                          >  1:753187/1‑143 (MQ=255)
                                                ggAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACcgc                                                                                                 >  2:229385/1‑150 (MQ=255)
                                                  aaGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAATTAGCTACACACTTGCTGACTTCACCGCGa                                                                                               <  1:1197871/150‑1 (MQ=255)
                                                  aaGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGa                                                                                               <  1:949774/150‑1 (MQ=255)
                                                              aTCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCCTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATc                                                                                   <  1:727904/150‑1 (MQ=255)
                                                                     ccACGACCGGGCCCACCGGGTGGTATTCGATACGGGCTTTTTCGCCCTTGGTCGGCTCGGCGGCGAGCCTGTCCGTGCAGTGCTCGGCGTACCAGCTCCAAAGGCTTACTCACTTGCTGACTTCACCGTGAGTCTGGGCGTTCGGCTTGc                                                                            <  2:372610/150‑1 (MQ=255)
                                                                      cAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTCCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGcc                                                                           >  2:335019/1‑150 (MQ=255)
                                                                              ggCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTcgc                                                                   <  1:569434/150‑1 (MQ=255)
                                                                                cccAGCGGCGGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCgg                                                                 <  1:1019367/150‑1 (MQ=255)
                                                                                     cgggcggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGAt                                                            <  1:906282/150‑1 (MQ=255)
                                                                                       ggcggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTgg                                                                                         <  1:143506/119‑1 (MQ=255)
                                                                                       ggcggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTgg                                                                                         >  2:143506/1‑119 (MQ=255)
                                                                                        gcggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGc                                                                                >  1:492845/1‑127 (MQ=255)
                                                                                        gcggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGc                                                                                <  2:492845/127‑1 (MQ=255)
                                                                                                 gATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGCCGAGCACGGCCCGGCCATGCTCgc                                                                                                                                           >  1:769401/1‑58 (MQ=18)
                                                                                                 gATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGCCGAGCACGGCCCGGCCATGCTCgc                                                                                                                                           <  2:769401/59‑2 (MQ=18)
                                                                                                    aCGGGCTTTTCCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGc                                             <  1:335019/150‑1 (MQ=255)
                                                                                                          tttttCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGc                                       <  2:1185889/150‑1 (MQ=255)
                                                                                                          tttttCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGc                                       <  1:252683/150‑1 (MQ=255)
                                                                                                                       gTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAgg                          <  1:708691/150‑1 (MQ=255)
                                                                                                                         cggctcggcGGCACGACTGGCCGGGCCGTGCCCGGCGTACCAGCGGCAAAGGCGTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGgt                        <  1:1227048/150‑1 (MQ=255)
                                                                                                                           gctcggcGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGgtgt                      >  1:889064/1‑150 (MQ=255)
                                                                                                                             tcggcGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGATGATGATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCg                    >  1:1221166/1‑150 (MQ=255)
                                                                                                                               ggcggcGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCggc                  <  2:739650/150‑1 (MQ=255)
                                                                                                                                    cgAGCACGGCCCGGCCATGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAAtt                                                                       >  1:1023425/1‑92 (MQ=255)
                                                                                                                                    cgAGCACGGCCCGGCCATGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAAtt                                                                       <  2:1023425/92‑1 (MQ=255)
                                                                                                                                        cATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCggccagggccag         >  2:321854/1‑150 (MQ=255)
                                                                                                                                            gccgggccgTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGaaa                                               <  2:563589/108‑1 (MQ=255)
                                                                                                                                            gccgggccgTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGaaa                                               >  1:563589/1‑108 (MQ=255)
                                                                                                                                               gggccgTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCGGCCAGGGCCAGCAGGTAc  <  1:338817/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TACCCGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCGGCCAGGGCCAGCAGGTAC  >  NC_002947/3570448‑3570740

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑