Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 864,299 A→G 7.1% I220V (ATT→GTT)  fsaA → fructose‑6‑phosphate aldolase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913864,2990AG 7.1% 166.3 / 3.2 56I220V (ATT→GTT) fsaAfructose‑6‑phosphate aldolase 1
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base G (2/2);  total (27/29)
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.76e-01

GCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTCATC  >  NC_000913/864151‑864446
                                                                                                                                                    |                                                                                                                                                   
gcgTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGa                                                                                                                                                     <  7:118324/149‑1 (MQ=255)
 cgTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGAt                                                                                                                                                    <  4:129657/149‑1 (MQ=255)
         gcTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGGCGATTTAACCgg                                                                                                                                            <  2:95814/149‑1 (MQ=255)
                  cTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCt                                                                                                                                     <  7:85799/147‑1 (MQ=255)
                  cTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTg                                                                                                                                    >  6:192809/1‑148 (MQ=255)
                       tGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCt                                                                                                                              >  6:159812/1‑149 (MQ=255)
                       tGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCt                                                                                                                              >  4:181790/1‑149 (MQ=255)
                         gCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTc                                                                                                                             <  2:173649/148‑1 (MQ=255)
                         gCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCa                                                                                                                            <  8:1068/149‑1 (MQ=255)
                             gATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGGTTTAACCGGACGTTCTGCATCCTCATaaa                                                                                                                        >  5:143468/1‑149 (MQ=255)
                              aTGTGAATCAATTACTCTGCCACTGGATTTGGCACAACAGATGATTAGCTATCCTGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAAt                                                                                                                       <  2:17917/149‑1 (MQ=255)
                                             tctGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGc                                                                                                        >  6:248220/1‑149 (MQ=255)
                                                  cACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGt                                                                                                   >  3:183527/1‑149 (MQ=255)
                                                   aCTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCg                                                                                                       <  4:18464/144‑1 (MQ=255)
                                                   aCTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTg                                                                                                  <  2:27918/149‑1 (MQ=255)
                                                   aCTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGAAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCg                                                                                                       >  3:18464/1‑144 (MQ=255)
                                                                 aacaGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAg                                                                                    <  4:264391/149‑1 (MQ=255)
                                                                 aacaGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAg                                                                                    >  5:69354/1‑149 (MQ=255)
                                                                   caGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGc                                                                                  <  6:69354/149‑1 (MQ=255)
                                                                    aGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGcc                                                                                 >  1:261009/1‑149 (MQ=255)
                                                                    aGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGcc                                                                                 >  1:287876/1‑149 (MQ=255)
                                                                     gatgatTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCTCGTCTTATCAGGCCt                                                                                >  1:273788/1‑149 (MQ=255)
                                                                      atgatTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTg                                                                               <  6:287475/149‑1 (MQ=255)
                                                                      atgatTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTg                                                                               <  6:143794/149‑1 (MQ=255)
                                                                      atgatTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTg                                                                               <  5:37681/149‑1 (MQ=255)
                                                                         atTAGCTATCCGGCGGTTGATGCCGCGGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGa                                                                             <  1:158240/148‑1 (MQ=255)
                                                                          ttAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAgg                                                                           >  6:142038/1‑149 (MQ=255)
                                                                            aGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGCCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTg                                                                         >  8:326874/1‑149 (MQ=255)
                                                                             gCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTgg                                                                        <  5:142038/149‑1 (MQ=255)
                                                                             gCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTgg                                                                        >  7:36269/1‑149 (MQ=255)
                                                                              cTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGt                                                                                            >  6:281965/1‑128 (MQ=255)
                                                                                   cggcggTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTTCTGATGACGTGCCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAAtt                                                                   <  7:326874/148‑1 (MQ=255)
                                                                                     gcggTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACt                                                                >  6:290383/1‑149 (MQ=255)
                                                                                        gTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGcc                                                             >  7:164647/1‑149 (MQ=255)
                                                                                                       gCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTc                                                                                           >  6:111430/1‑104 (MQ=255)
                                                                                                       gCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTc                                                                                           <  5:111430/104‑1 (MQ=255)
                                                                                                         gAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCAcccccc                                            >  6:50250/1‑149 (MQ=255)
                                                                                                         gAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACAcccc                                            >  1:49376/1‑149 (MQ=255)
                                                                                                         gAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACAcccc                                            >  2:214337/1‑149 (MQ=255)
                                                                                                                gAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTa                                      <  5:192809/148‑1 (MQ=255)
                                                                                                                gAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTAc                                     >  8:201907/1‑149 (MQ=255)
                                                                                                                 aGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGGTTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTAc                                     <  6:143468/148‑1 (MQ=255)
                                                                                                                 aGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACg                                    >  5:2276/1‑149 (MQ=255)
                                                                                                                  gCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACg                                    >  3:32248/1‑148 (MQ=255)
                                                                                                                     ggACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATc                                <  3:276540/149‑1 (MQ=255)
                                                                                                                           gCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTc                          <  8:164647/149‑1 (MQ=255)
                                                                                                                                gAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCTCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTc                     <  2:273788/149‑1 (MQ=255)
                                                                                                                                 aGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTcg                    <  5:290383/149‑1 (MQ=255)
                                                                                                                                  gCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTcgc                   >  3:129928/1‑149 (MQ=255)
                                                                                                                                       tGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGaaa               <  1:277499/148‑1 (MQ=255)
                                                                                                                                         gCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACt             <  4:183527/148‑1 (MQ=255)
                                                                                                                                          cAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGa           >  7:81425/1‑149 (MQ=255)
                                                                                                                                          cAGAACGTCCGGTTAAATCGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTcgc                   <  7:192937/141‑1 (MQ=255)
                                                                                                                                          cAGAACGTCCGGTTAAATCGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTcgc                   >  8:192937/1‑141 (MQ=255)
                                                                                                                                             aaCGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAg        <  5:159812/149‑1 (MQ=255)
                                                                                                                                                   gATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTcatc  <  2:287876/149‑1 (MQ=255)
                                                                                                                                                   gATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTcatc  <  3:166358/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCTATCCGGCGGTTGATGCCGCTGTGGCGAAGTTTGAGCAGGACTGGCAGGGAGCGTTTGGCAGAACGTCGATTTAACCGGACGTTCTGCATCCTCATAAATTGCTGATGACGTGGCGGAGTGCCGCGTCTTATCAGGCCTGGAGGTGGCAATTACTGCCCACACACCTCACACCCCGGATTACGCATCAGTTTCATTTCGCGAAACTGACAGGTCATC  >  NC_000913/864151‑864446

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