Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,570,586 T→C 8.3% M108V (ATG→GTG)  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,5860TC 8.3% 133.8 / 3.1 48M108V (ATG→GTG) sad‑IINAD+‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (15/29);  new base C (2/2);  total (17/31)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

TACCCGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCGGCCAGGGCCAG  >  NC_002947/3570448‑3570733
                                                                                                                                          |                                                                                                                                                   
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)
            tCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATggccgg                                                                                                                                               <  2:942250/133‑1 (MQ=255)
            tCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATggccgg                                                                                                                                               >  1:942250/1‑133 (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                                                                                                   >  1:753187/1‑143 (MQ=255)
                                              cGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGa                                                                                                   <  2:753187/143‑1 (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                                                                                  >  2:143506/1‑119 (MQ=255)
                                                                                       ggcggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTgg                                                                                  <  1:143506/119‑1 (MQ=255)
                                                                                        gcggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGc                                                                         >  1:492845/1‑127 (MQ=255)
                                                                                        gcggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGc                                                                         <  2:492845/127‑1 (MQ=255)
                                                                                                 gATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGCCGAGCACGGCCCGGCCATGCTCgc                                                                                                                                    <  2:769401/59‑2 (MQ=18)
                                                                                                 gATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGCCGAGCACGGCCCGGCCATGCTCgc                                                                                                                                    >  1:769401/1‑58 (MQ=18)
                                                                                                    aCGGGCTTTTCCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGc                                      <  1:335019/150‑1 (MQ=255)
                                                                                                          tttttCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGc                                <  1:252683/150‑1 (MQ=255)
                                                                                                          tttttCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGc                                <  2:1185889/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                                                                <  2:1023425/92‑1 (MQ=255)
                                                                                                                                    cgAGCACGGCCCGGCCATGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAAtt                                                                >  1:1023425/1‑92 (MQ=255)
                                                                                                                                        cATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCggccagggccag  >  2:321854/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TACCCGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCGGCCAGGGCCAG  >  NC_002947/3570448‑3570733

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