Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,364,196 A→G 9.1% E485G (GAA→GGA)  puuC → gamma‑glutamyl‑gamma‑aminobutyraldehyde dehydrogenase; succinate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,364,1960AG 9.1% 122.9 / 3.5 44E485G (GAA→GGA) puuCgamma‑glutamyl‑gamma‑aminobutyraldehyde dehydrogenase; succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (17/23);  new base G (2/2);  total (19/25)
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.30e-01

GCGCGCACCGCATGAGCCGACGCCTGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCTGCGACGTTTGCGTGGTTGGCGGCGGCTATAC  >  NC_000913/1364049‑1364344
                                                                                                                                                   |                                                                                                                                                    
gcgcgcACCGCATGAGCCGACGCCTGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGTCGATATGAACGTGCCGTTTGGAGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATTCCCTTGAAAAATTAactgaa                                                                                                                                                     <  1:254693/149‑1 (MQ=255)
      aCCGCATGAGCCGACGCCTGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGaaa                                                                                                                                               <  4:53350/149‑1 (MQ=255)
      aCCGCATGAGCCGACGCCTGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGaaa                                                                                                                                               <  2:178874/149‑1 (MQ=255)
       ccGCATGAGCCGACGCCTGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGaaaa                                                                                                                                              <  5:225704/149‑1 (MQ=255)
       ccGCATGAGCCGACGCCTGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGaaaa                                                                                                                                              <  5:25764/149‑1 (MQ=255)
                      ccTGAAAGCCGGTTCCGGCTTCTTCAATAACTACACAGACGGCGATATTACAGTGCCGTTTGGCGGCTATAAGCAGAGCTGCAACGGTCGAGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGcc                                                                                                                               <  5:250764/149‑1 (MQ=255)
                       cTGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGGACTGAAAACCATCTGGATAAGCCt                                                                                                                              >  1:46112/1‑149 (MQ=255)
                       cTGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGGACTGAAAACCATCTGGATAAGCCt                                                                                                                              >  3:278109/1‑149 (MQ=255)
                        tGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTg                                                                                                                             >  5:301631/1‑149 (MQ=255)
                          aaaGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGGACTGAAAACCATCTGGATAAGCCTGGa                                                                                                                           <  2:46112/149‑1 (MQ=255)
                          aaaGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGGACTGAAAACCATCTGGATAAGCCTGGa                                                                                                                           <  4:278109/149‑1 (MQ=255)
                               cGGTTCCGTCTTCGTCAACACCTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGACCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGATGCCt                                                                                                                      <  2:93057/149‑1 (MQ=255)
                                    ccGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGaaa                                                                                                                  <  5:241169/148‑1 (MQ=255)
                                      gTCTTCGTCAATAACTACCACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGCAAAATTCACTGAACTGAAAACCATCTGGATACTCCTGGAGGCCTGAAATGa                                                                                                               >  2:236296/1‑149 (MQ=255)
                                           cGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGaaa                                                                                                                  >  6:245409/1‑141 (MQ=255)
                                           cGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGaaa                                                                                                                  <  5:245409/141‑1 (MQ=255)
                                                aTAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATAc                                                                                                     >  5:139461/1‑149 (MQ=255)
                                                   aCTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACca                                                                                                   <  7:138736/148‑1 (MQ=255)
                                                    cTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGCCCGAACATACcag                                                                                                  >  3:164480/1‑148 (MQ=255)
                                                       cAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACcagca                                                                                                <  7:220095/147‑1 (MQ=255)
                                                          cgacgGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTtac                                                                                           <  3:221803/149‑1 (MQ=255)
                                                                      gACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGt                                                                               >  5:145115/1‑149 (MQ=255)
                                                                          gTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTtac                                                                                           <  5:222092/133‑1 (MQ=255)
                                                                          gTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTtac                                                                                           >  6:222092/1‑133 (MQ=255)
                                                                                tttGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTg                                                                              >  6:23763/1‑140 (MQ=255)
                                                                                tttGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTg                                                                              <  5:23763/140‑1 (MQ=255)
                                                                                  tGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGaataaata                                                                    <  1:279745/148‑1 (MQ=255)
                                                                                    gcggcTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATg                                                                  <  1:286574/148‑1 (MQ=255)
                                                                                         tataAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCCCCAt                                                            >  7:115819/1‑149 (MQ=255)
                                                                                           taAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTc                                                          <  6:301631/149‑1 (MQ=255)
                                                                                             aGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGa                                                        <  2:245207/149‑1 (MQ=255)
                                                                                               cAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGaca                                                      >  2:95147/1‑149 (MQ=255)
                                                                                                             acgcgACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTc                                         >  5:8966/2‑148 (MQ=255)
                                                                                                                 gACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCa                                                             <  7:167361/124‑1 (MQ=255)
                                                                                                                 gACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCa                                                             >  8:167361/1‑124 (MQ=255)
                                                                                                                    aaaTCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCt                                 >  4:251765/1‑149 (MQ=255)
                                                                                                                    aaaTCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCt                                 <  3:251765/149‑1 (MQ=255)
                                                                                                                             cATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCTGCGACGttt                        <  6:8966/149‑1 (MQ=255)
                                                                                                                                  ccTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCTGCGACGTTTGCGTg                   <  4:164480/149‑1 (MQ=255)
                                                                                                                                   cTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCTGCGACGTTTGCGTgg                  >  2:113657/1‑149 (MQ=255)
                                                                                                                                   cTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCTGCGACGTTTGCGTgg                  >  7:243099/1‑149 (MQ=255)
                                                                                                                                     tGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCAc                                   >  6:167537/1‑130 (MQ=255)
                                                                                                                                     tGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCAc                                   <  5:167537/130‑1 (MQ=255)
                                                                                                                                      gAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCTGCGACGTTTGCGTGGTTg               >  5:244077/1‑149 (MQ=255)
                                                                                                                                              cACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCTGCGACGTTTGCGTGGTTggcggcgg        <  1:95147/148‑1 (MQ=255)
                                                                                                                                                   aactgaaAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCTGCGACGTTTGCGTGGTTGGCGGCGGCTATAc  >  3:171926/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGCGCACCGCATGAGCCGACGCCTGAAAGCCGGTTCCGTCTTCGTCAATAACTACAACGACGGCGATATGACCGTGCCGTTTGGCGGCTATAAGCAGAGCGGCAACGGTCGCGACAAATCCCTGCATGCCCTTGAAAAATTCACTGAACTGAAAACCATCTGGATAAGCCTGGAGGCCTGAAATGACCGAACATACCAGCAGTTACTACGCCGCCAGTGCGAATAAATATGCACCATTCGACACGCTGAATGAGTCAATCACCTGCGACGTTTGCGTGGTTGGCGGCGGCTATAC  >  NC_000913/1364049‑1364344

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑