Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,516,507 T→C 6.9% G310G (GGT→GGC)  patD → gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,516,5070TC 6.9% 188.2 / 2.7 58G310G (GGT→GGC) patDgamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base C (2/2);  total (30/28)
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.97e-01

TTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATT  >  NC_000913/1516359‑1516652
                                                                                                                                                    |                                                                                                                                                 
ttGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGt                                                                                                                                                   <  1:152873/149‑1 (MQ=255)
           gTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGa                                                                                                                                        >  3:80557/1‑149 (MQ=255)
           gTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGa                                                                                                                                        >  1:398760/1‑149 (MQ=255)
            tCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGAc                                                                                                                                       <  4:40720/149‑1 (MQ=255)
              gAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGa                                                                                                                                     >  3:69258/1‑149 (MQ=255)
                            aTTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTg                                                                                                                       <  1:296844/149‑1 (MQ=255)
                                   tATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCttt                                                                                                                >  2:285007/1‑149 (MQ=255)
                                   tATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCttt                                                                                                                >  2:268465/1‑149 (MQ=255)
                                     ttACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTaa                                                                                                              >  3:339914/1‑149 (MQ=255)
                                     ttACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTaa                                                                                                              <  2:45654/149‑1 (MQ=255)
                                           tGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGc                                                                                                        <  2:196824/149‑1 (MQ=255)
                                               ggACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGCGCACCAGAtt                                                                                                                                         <  4:12065/112‑2 (MQ=255)
                                               ggACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGCGCACCAGAtt                                                                                                                                         >  3:12065/1‑111 (MQ=255)
                                                gACAGGATTGTGCTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcg                                                                                                   >  1:251074/1‑149 (MQ=255)
                                                gACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAGCGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcg                                                                                                   >  3:345584/1‑149 (MQ=255)
                                                gACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcg                                                                                                   >  4:128285/1‑149 (MQ=255)
                                                   aGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGa                                                                                                                                     <  2:26161/112‑1 (MQ=255)
                                                   aGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGa                                                                                                                                     >  1:26161/1‑112 (MQ=255)
                                                         gTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGGGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGaa                                                                                          <  3:159369/149‑1 (MQ=255)
                                                         gTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGaa                                                                                          <  2:430164/149‑1 (MQ=255)
                                                          tACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAAc                                                                                         <  4:268043/149‑1 (MQ=255)
                                                          tACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAAc                                                                                         >  4:125154/1‑149 (MQ=255)
                                                             tGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAAcgcg                                                                                      >  3:189370/1‑149 (MQ=255)
                                                                ggCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCg                                                                                   <  1:261434/149‑1 (MQ=255)
                                                                   ttGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCa                                                                                >  2:31655/1‑149 (MQ=255)
                                                                          aTCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGt                                                                         >  3:173002/1‑149 (MQ=255)
                                                                          aTCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGt                                                                         <  4:80557/149‑1 (MQ=255)
                                                                                gcgcAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTagaaga                                                                   <  3:125154/149‑1 (MQ=255)
                                                                                    aaaaaGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGGAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCg                                                               >  3:406511/1‑149 (MQ=255)
                                                                                     aaaaGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGGAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGa                                                              <  4:406511/149‑1 (MQ=255)
                                                                                      aaaGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGaa                                                             <  1:268465/149‑1 (MQ=255)
                                                                                       aaGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGaaa                                                            <  4:173002/149‑1 (MQ=255)
                                                                                       aaGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGaaa                                                            <  1:268415/149‑1 (MQ=255)
                                                                                        aGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGa                                                                                                                            >  3:27429/1‑84 (MQ=255)
                                                                                        aGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGa                                                                                                                            <  4:27429/84‑1 (MQ=255)
                                                                                            aTTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGAc                                                       >  2:271893/1‑149 (MQ=255)
                                                                                             tttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACa                                                      >  2:246234/1‑149 (MQ=255)
                                                                                               tACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCa                                                                           <  3:170109/126‑1 (MQ=255)
                                                                                               tACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAgg                                                    <  1:285007/149‑1 (MQ=255)
                                                                                               tACGATACGCTGGAGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCa                                                                           >  4:170109/1‑126 (MQ=255)
                                                                                                       gCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCa                                            >  3:417212/1‑149 (MQ=255)
                                                                                                         tggtggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCaaa                                          <  2:251074/149‑1 (MQ=255)
                                                                                                           gtggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGt                                        <  4:339914/149‑1 (MQ=255)
                                                                                                               aaaaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATc                                    >  3:205427/1‑149 (MQ=255)
                                                                                                                aaaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCa                                   >  3:229470/1‑149 (MQ=255)
                                                                                                                          gctgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGa                         <  1:31655/149‑1 (MQ=255)
                                                                                                                            tgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGaaa                       >  3:317638/1‑149 (MQ=255)
                                                                                                                                     aaCGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTagaaga                                                                   >  4:99349/1‑96 (MQ=255)
                                                                                                                                     aaCGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTagaaga                                                                   <  3:99349/96‑1 (MQ=255)
                                                                                                                                      aCGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAggg             >  2:467545/1‑149 (MQ=255)
                                                                                                                                      aCGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAggg             <  1:246234/149‑1 (MQ=255)
                                                                                                                                        gTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGaaaa                      >  2:188765/1‑138 (MQ=255)
                                                                                                                                        gTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGaaaa                      <  1:188765/138‑1 (MQ=255)
                                                                                                                                         tcgtcaTCTGGCGCACCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACTTCAAAGTGa                                      <  4:50365/116‑1 (MQ=255)
                                                                                                                                         tcgtcaTCTGGCGCACCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGa                                      >  3:50365/6‑121 (MQ=255)
                                                                                                                                           aaaaTCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATg        <  2:33286/149‑1 (MQ=255)
                                                                                                                                               tCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCa                                                                           <  3:351038/78‑1 (MQ=255)
                                                                                                                                               tCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCa                                                                           >  4:351038/1‑78 (MQ=255)
                                                                                                                                                 tGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTAtt  >  3:75302/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATT  >  NC_000913/1516359‑1516652

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