Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,052 T→C 11.7% 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,0520TC11.7% 130.9 / 7.5 51L256P (CTG→CCG) gapBglyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (27/18);  new base C (2/4);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

ACACCATCACTGCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGCCATCGTGCC  >  NC_002947/2455905‑2456199
                                                                                                                                                   |                                                                                                                                                   
acacCATCACTGCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTTGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTg                                                                                                                                                    <  1:81835/149‑1 (MQ=255)
 cacCATCACTGCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTgg                                                                                                                                                   <  3:211375/149‑1 (MQ=255)
    cATCACTGCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCt                                                                                                                                                >  7:173330/1‑149 (MQ=255)
           gCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCggg                                                                                                                                         >  6:344327/1‑149 (MQ=255)
            ccAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCTACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCgg                                                                                                                                          <  8:306575/147‑1 (MQ=255)
            ccAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCgg                                                                                                                                          <  1:27511/147‑1 (MQ=255)
             cAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGtg                                                                                                                                       >  3:113577/1‑149 (MQ=255)
              aaCGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGtgc                                                                                                                                      >  6:116560/1‑149 (MQ=255)
              aaCGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCCGGGCAGGCCAGGtgc                                                                                                                                      <  6:24842/149‑1 (MQ=255)
                 ctcAACCTCATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCCGGGCAGGCCAGGtgctgg                                                                                                                                   <  7:307605/147‑2 (MQ=255)
                    aaCCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCCGGGCAGGCCAGGtgctgg                                                                                                                                   >  8:303167/1‑145 (MQ=255)
                    aaCCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCCGGGCAGGCCAGGtgctgg                                                                                                                                   <  7:303167/146‑2 (MQ=255)
                    aaCCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCCGGGCAGGCCAGGtgctgg                                                                                                                                   >  8:307605/1‑145 (MQ=255)
                      ccTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGG‑CAGCACCCGGGCAGGCCAGGTGCTGgcccagg                                                                                                                             <  2:103773/149‑5 (MQ=38)
                                 ttATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCAcc                                                                                                                   >  7:261392/1‑149 (MQ=255)
                                     cTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCAc                                                                                                               >  4:198404/1‑149 (MQ=255)
                                        cGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCAc                                                                                                               <  7:90567/146‑1 (MQ=255)
                                        cGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCAc                                                                                                               >  8:90567/1‑146 (MQ=255)
                                           cAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGa                                                                                                                           <  7:73022/131‑1 (MQ=255)
                                           cAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGa                                                                                                                           >  8:73022/1‑131 (MQ=255)
                                           cAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAAAGCGCTCATCGACGACAACACCGGCGTAAGCACTGACGCCGACGGCCTGGGCAAGCACCTGGCCTGACAGGGTGCTGCCCGCGGGATCCTCACAGAACCTgga                                                                                                          >  8:98644/1‑147 (MQ=255)
                                                gCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGgc                                                                                                    <  4:113577/149‑1 (MQ=255)
                                                              gACTACACCCAGTACGGCATCGACAAAGCCCTCATCGAAAACAAAAACGACGTATGCCGTGACGACCACCGCCTGGGCCAGCACATGTCCTTCACGGGTGCTGACCGCGTGATCCTCAACCAACATGACAAGGACGCGCTGAAGAACa                                                                                       >  4:309000/1‑148 (MQ=255)
                                                                   cacCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATc                                                                                                                         <  4:217371/109‑1 (MQ=255)
                                                                   cacCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATc                                                                                                                         >  3:217371/1‑109 (MQ=255)
                                                                    acCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCa                                                                                >  4:266362/1‑149 (MQ=255)
                                                                            cGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATc                                                                         <  3:266362/148‑1 (MQ=255)
                                                                             ggCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCaa                                                                       >  4:38086/1‑149 (MQ=255)
                                                                                aTCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCa                                                                    >  5:32144/1‑149 (MQ=255)
                                                                                    aCAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACTTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGt                                                                >  7:48565/1‑149 (MQ=255)
                                                                                    aCAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGt                                                                >  8:40395/1‑149 (MQ=255)
                                                                                    aCAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGt                                                                <  8:261392/149‑1 (MQ=255)
                                                                                                tcgtcgACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGcc                                                    >  8:263084/1‑149 (MQ=255)
                                                                                                     gACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACTTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGa                                               <  8:48565/149‑1 (MQ=255)
                                                                                                            ccGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGatc                                        <  7:69024/149‑1 (MQ=255)
                                                                                                             cGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGatca                                       >  4:42170/1‑149 (MQ=255)
                                                                                                               gCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGatcatc                                     >  1:188802/1‑149 (MQ=255)
                                                                                                               gCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGATGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAAACACTGTGATATCAATGCCGATGACAAGataat                                      >  2:253780/1‑148 (MQ=255)
                                                                                                                cGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCAtct                                    <  7:40395/149‑1 (MQ=255)
                                                                                                                  tATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTcg                                  <  3:198404/149‑1 (MQ=255)
                                                                                                                  tATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTcg                                  >  4:324162/1‑149 (MQ=255)
                                                                                                                        gTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGcc                            >  8:221576/1‑149 (MQ=255)
                                                                                                                         tGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACg                                                                              <  2:223949/98‑1 (MQ=255)
                                                                                                                         tGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACg                                                                              >  1:223949/1‑98 (MQ=255)
                                                                                                                             gCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTcggcg                               >  8:138976/1‑141 (MQ=255)
                                                                                                                             gCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTcggcg                               <  7:138976/141‑1 (MQ=255)
                                                                                                                                 aCGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGcacc                   <  2:188802/149‑1 (MQ=255)
                                                                                                                                              gCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGAAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGCCATCg      >  2:73530/1‑149 (MQ=255)
                                                                                                                                                aCCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGAGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGCCATCGTg    >  4:292159/1‑149 (MQ=255)
                                                                                                                                                  cTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGCCATCGTGcc  >  8:238772/1‑149 (MQ=255)
                                                                                                                                                  cTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGCCATCGTGcc  <  7:238772/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACACCATCACTGCCAACGGCAACCTGATCCAGGTTATCTACGCCAAGAGCCCTAGCGAAGTCGACTACACCCAGTACGGCATCGACAATGCGCTCATCGTCGACAACACCGGCGTATGGCGTGACGCCGACGGCCTGGGCCAGCACCTGGCCTGCCCGGGTGCTGCCCGCGTGATCCTCACCGCACCTGGCAAGGGCGCGCTGAAGAACATCGTGCACGGTATCAACCACGGTGATATCACTGCCGATGACAAGATCATCTCGGCGGCCTCCTGCACCACCAACGCCATCGTGCC  >  NC_002947/2455905‑2456199

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