Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,052 T→C 9.1% L256P (CTG→CCG)  gapB → glyceraldehyde‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,456,0520TC 9.1% 119.1 / 2.5 45L256P (CTG→CCG) gapBglyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base C (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACACCATCACTGCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGC  >  NC_002947/2455904‑2456190
                                                                                                                                                    |                                                                                                                                          
aacaCCATCACTGCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCt                                                                                                                                            >  6:218752/1‑149 (MQ=255)
         aCTGCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGccc                                                                                                                                   <  5:218752/149‑1 (MQ=255)
             ccAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCggg                                                                                                                                <  7:146295/148‑1 (MQ=255)
             ccAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGt                                                                                                                               <  7:267128/149‑1 (MQ=255)
                 cGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCCGGGCAGGCCAGGtgctg                                                                                                                           <  7:32977/149‑1 (MQ=255)
                         tGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGtg                                                                                                                              >  4:267616/1‑138 (MQ=255)
                         tGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGtg                                                                                                                              <  3:267616/138‑1 (MQ=255)
                             ccAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATcc                                                                                                               >  3:283844/1‑149 (MQ=255)
                             ccAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGACTGCCCGGGTGCTGCCCGCGTGATcc                                                                                                               >  5:269746/1‑149 (MQ=255)
                                        aCGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCt                                                                                                    >  3:146927/1‑149 (MQ=255)
                                             aaGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGG‑CAGCACCCGGGCAGGCCAGGTGCTGgccca                                                                                                                      >  4:16967/1‑123 (MQ=38)
                                             aaGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGG‑CAGCACCCGGGCAGGCCAGGTGCTGgccca                                                                                                                      <  3:16967/125‑3 (MQ=38)
                                               gagCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAgg                                                                                             <  4:216231/149‑1 (MQ=255)
                                                      aGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCAGCCCGCCTGATCCTCACCGCACCTGGCAAGGGCGCGCt                                                                                      >  2:45020/1‑149 (MQ=255)
                                                       gCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGgcgcg                                                                                        >  2:174124/1‑146 (MQ=255)
                                                       gCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGgcgcg                                                                                        <  1:174124/146‑1 (MQ=255)
                                                         gAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTgaa                                                                                   <  8:21573/149‑1 (MQ=255)
                                                           aGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTgaaga                                                                                 <  1:23764/149‑1 (MQ=255)
                                                               gACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACAt                                                                             >  1:232450/1‑149 (MQ=255)
                                                                               gCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAAc                                                             >  6:160666/1‑149 (MQ=255)
                                                                                     aCAATGCGCTCATCATCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGt                                                       >  7:133872/1‑149 (MQ=255)
                                                                                           cgcTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCAGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGACCATCGTGCACGGTATCAACCACGGTGatat                                                  >  7:21030/1‑148 (MQ=255)
                                                                                            gcTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCa                                                >  3:127205/1‑149 (MQ=255)
                                                                                             cTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATc                                                                >  4:57483/1‑132 (MQ=255)
                                                                                             cTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCAGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATc                                                                <  3:57483/132‑1 (MQ=255)
                                                                                               cATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTg                                             >  3:4278/1‑149 (MQ=255)
                                                                                               cATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTg                                             >  5:73149/1‑149 (MQ=255)
                                                                                                aTCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACt                                              <  2:227921/147‑1 (MQ=255)
                                                                                                aTCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGc                                            <  5:160666/149‑1 (MQ=255)
                                                                                                aTCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGc                                            >  5:160758/1‑149 (MQ=255)
                                                                                                  cgtcgACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCcgc                                                                                                                      >  1:23552/1‑73 (MQ=255)
                                                                                                  cgtcgACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCcgc                                                                                                                      <  2:23552/73‑1 (MQ=255)
                                                                                                      gACATCACGGTCGTATTGCGTGACGCCGTCGGCCTGGGCCAGCACGTGGCCTGGGCGGGTGCTGGCCGGGTGATCCTGACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTTTCAAGCACGGTGATATCACTGCCGATGa                                      <  8:21030/149‑1 (MQ=255)
                                                                                                      gACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGa                                      <  4:127205/149‑1 (MQ=255)
                                                                                                       acaacaCCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGAc                                     <  4:4278/149‑1 (MQ=255)
                                                                                                         aacaCCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACaa                                   <  8:133872/149‑1 (MQ=255)
                                                                                                                     tGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTcggc                       <  4:283844/149‑1 (MQ=255)
                                                                                                                             cGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGcctcct               >  6:88814/1‑149 (MQ=255)
                                                                                                                              gCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTg              <  4:146927/149‑1 (MQ=255)
                                                                                                                                cGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGca            <  2:232450/149‑1 (MQ=255)
                                                                                                                                          gggcAGCACCCGGGCAGGCCAGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGc  >  5:117083/4‑149 (MQ=255)
                                                                                                                                           gCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGc  >  1:283391/1‑148 (MQ=255)
                                                                                                                                             cAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTg              <  5:283674/134‑1 (MQ=255)
                                                                                                                                             cAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTg              >  6:283674/1‑134 (MQ=255)
                                                                                                                                              aGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCaa                                                                                               >  7:174881/1‑52 (MQ=255)
                                                                                                                                              aGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCaa                                                                                               <  8:174881/52‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
AACACCATCACTGCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGC  >  NC_002947/2455904‑2456190

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