Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,455,455 T→A 8.9% I57N (ATT→AAT)  gapB → glyceraldehyde‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,455,4550TA 8.9% 125.9 / 3.5 45I57N (ATT→AAT) gapBglyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (27/14);  new base A (2/2);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

CCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAGCTGGCCAACAAGTTCAA  >  NC_002947/2455313‑2455602
                                                                                                                                              |                                                                                                                                                   
cccGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCaa                                                                                                                                               <  2:359166/149‑1 (MQ=255)
     ccAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGgcgc                                                                                                                                          <  4:292273/149‑1 (MQ=255)
             ttGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCttt                                                                                                                                  >  3:242679/1‑149 (MQ=255)
                     tGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCa                                                                                                                          >  7:214352/1‑149 (MQ=255)
                     tGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTAGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCa                                                                                                                          >  7:108510/1‑149 (MQ=255)
                        aTTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCg                                                                                                                       >  3:156755/1‑149 (MQ=255)
                                 gAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACcga                                                                                                              >  3:316642/1‑149 (MQ=255)
                                  aaGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACcgac                                                                                                             <  1:25636/149‑1 (MQ=255)
                                           cTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGaa                                                                                                    <  7:291768/149‑1 (MQ=255)
                                            tGAAGCGATGATCCCGCTTATCGGTCAGCTATACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCAACCGCCAAACCGACGAAACCGAAc                                                                                                   >  6:100567/1‑149 (MQ=255)
                                                        cccGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGCTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCAc                                                                                        >  7:140566/1‑148 (MQ=255)
                                                                   gTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAACCCGAACTGTCCGTCCACGAGACATTcccc                                                                            >  5:380645/1‑149 (MQ=255)
                                                                         tctACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCt                                                                       >  2:188643/1‑148 (MQ=255)
                                                                             ccGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTcccc                                                                            <  2:34422/139‑1 (MQ=255)
                                                                             ccGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTcccc                                                                            >  1:34422/1‑139 (MQ=255)
                                                                              cGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGc                                                                 >  2:281092/1‑149 (MQ=255)
                                                                               gCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAg                                                                                                                                              <  7:104981/71‑1 (MQ=255)
                                                                               gCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAg                                                                                                                                              >  8:104981/1‑71 (MQ=255)
                                                                                aaaTAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCTACGAAACCGAACTGTACGTCCCCGAGACATTCACCCTGCTCAAGGCCa                                                               >  7:348306/2‑149 (MQ=255)
                                                                                         tgtgGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATgagcgagc                                                      <  4:242679/149‑1 (MQ=255)
                                                                                              tgAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGaaa                                                 <  8:108510/149‑1 (MQ=255)
                                                                                                  cTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAGCTgg                                             >  6:257850/1‑149 (MQ=255)
                                                                                                  cTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGAATCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTc                                                                               >  6:52005/1‑115 (MQ=255)
                                                                                                  cTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGAATCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTc                                                                               >  6:46568/1‑115 (MQ=255)
                                                                                                  cTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGAATCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTc                                                                               <  5:46568/115‑1 (MQ=255)
                                                                                                  cTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGAATCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTc                                                                               <  5:52005/115‑1 (MQ=255)
                                                                                                   tCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTggg                                            >  7:244372/1‑149 (MQ=255)
                                                                                                   tCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTggg                                            >  4:273887/1‑149 (MQ=255)
                                                                                                   tCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTggg                                            >  8:235177/1‑149 (MQ=255)
                                                                                                     gATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGgcg                                          >  3:218669/1‑149 (MQ=255)
                                                                                                       tttATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGcgcc                                        <  7:235177/149‑1 (MQ=255)
                                                                                                          aTGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGcgccgc                                      >  4:11588/1‑148 (MQ=255)
                                                                                                               cGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGt                                <  8:25446/149‑1 (MQ=255)
                                                                                                                      tGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGGTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTg                         >  1:310532/1‑149 (MQ=255)
                                                                                                                      tGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTg                         >  5:92234/1‑149 (MQ=255)
                                                                                                                        aTCAACCGTTCGGTTATCTCGATTCTCAAGGCTCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTgg                        <  8:140566/148‑1 (MQ=255)
                                                                                                                            aCCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGa                                                   <  1:370080/117‑1 (MQ=255)
                                                                                                                            aCCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGa                                                   >  2:370080/1‑117 (MQ=255)
                                                                                                                            aCCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAg                   >  3:286246/1‑149 (MQ=255)
                                                                                                                               gTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAGCTg                <  4:316642/149‑1 (MQ=255)
                                                                                                                                        tctcGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAGCTGGCcaacaa        <  4:156755/148‑1 (MQ=255)
                                                                                                                                         ctcGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAGCTGGCCAACAAGt      >  7:206672/1‑149 (MQ=255)
                                                                                                                                         ctcGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCCAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAGCTGGCCAACAAGt      >  7:177570/1‑149 (MQ=255)
                                                                                                                                          tcGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAGCTGGCCAACAAGtt     >  5:231370/1‑149 (MQ=255)
                                                                                                                                             aTTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAGCTGGCCAACAAGTTc    >  1:228976/1‑147 (MQ=255)
                                                                                                                                              ttCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAGCTGGCCAACAAGTTCaa  <  3:11588/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCGCAATAACAATGTGGTGAGCTCGATTTATGGCCGCAGCCTGATCAACCGTTCGGTTATCTCGATTCTCAAGGCGCACCGCTTTGCGCGCCACCGCCAAACCGACGAAACCGAACTGTCCGTCCACGAGACATTCCCCCTGCTCAAGGCCATGAGCGAGCTGAAACTGGGCGCCGCATCGGTCGACCTGGGCAAGCTGGCCAACAAGTTCAA  >  NC_002947/2455313‑2455602

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