Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,569,344 (T)6→5 intergenic (‑16/+2) glgC ← / ← glgX glucose‑1‑phosphate adenylyltransferase/limit dextrin alpha‑1,6‑glucohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,569,3390T.97.1% 134.6 / ‑1.5 36intergenic (‑11/+7)glgC/glgXglucose‑1‑phosphate adenylyltransferase/limit dextrin alpha‑1,6‑glucohydrolase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (20/14);  total (21/15)
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.73e-01

CAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCT  >  NC_000913/3569260‑3569415
                                                                               |                                                                            
cAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCtt                                                                             >  2:687941/1‑81 (MQ=255)
 aGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCttt                                                                            <  2:687884/81‑1 (MQ=255)
 aGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCtt                                                                             <  2:687883/80‑1 (MQ=255)
 aGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCtt                                                                             <  1:687882/80‑1 (MQ=255)
  gggCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCtttt                                                                           <  1:687885/81‑1 (MQ=255)
  gggCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCtttt                                                                           >  2:687939/1‑81 (MQ=255)
   ggCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCttttt                                                                          <  1:687886/81‑1 (MQ=255)
    gCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCttttta                                                                         <  2:687887/81‑2 (MQ=255)
        cAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTatcat                                                                    >  2:687961/1‑81 (MQ=255)
          gATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCAtct                                                                  <  1:687888/81‑1 (MQ=255)
            tttCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATctct                                                                <  1:687889/81‑1 (MQ=255)
                 aTGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGaa                                                            >  2:687959/1‑80 (MQ=255)
                 aTGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAac                                                           <  1:687891/81‑1 (MQ=255)
                 aTGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAac                                                           <  2:687890/81‑1 (MQ=255)
                    gccgCTGGCGCGCCAACCTTAAGTGATCGTTCTTCTCTAAACTAACCATGACTCACTCCTTTTAAACATCTCTGGAacaca                                                        <  2:687893/78‑1 (MQ=255)
                     cAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAacacac                                                       <  2:687892/81‑1 (MQ=255)
                            cgcgcCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATcc                                                <  1:687894/81‑1 (MQ=255)
                               gcCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCgtg                                             >  1:687943/1‑81 (MQ=255)
                                ccAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCgtgt                                            >  1:687934/1‑81 (MQ=255)
                                       tAAGTAATCGTTCTTCTCAAAACTAACAATACCAACCCCCTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTc                                     <  1:687895/80‑1 (MQ=255)
                                        aaGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTcc                                    >  2:687942/1‑81 (MQ=255)
                                         aGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTccc                                   >  2:688013/1‑81 (MQ=255)
                                         aGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTccc                                   >  1:687938/1‑81 (MQ=255)
                                          gTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCt                                  <  1:687896/81‑1 (MQ=255)
                                          gTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCt                                  >  1:687953/1‑81 (MQ=255)
                                            gATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTg                                 >  1:687962/1‑81 (MQ=255)
                                               cGTTCTTCTCTAAACTAACCATGACTAACTCCTTTCTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGc                                >  1:687940/1‑78 (MQ=255)
                                                 ttcttcTCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGAc                           <  1:687897/81‑1 (MQ=255)
                                                    ttctCTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGc                        >  1:687915/1‑81 (MQ=255)
                                                     tctcTAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGcc                       >  1:687947/1‑81 (MQ=255)
                                                         tAAACTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTaa                   >  1:687944/1‑81 (MQ=255)
                                                          aaaCTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGAGCAGGTCCCTGCCAGACAGCCGTAAt                  >  2:687950/1‑81 (MQ=255)
                                                           aaCTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACTCACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTa                    <  2:687898/78‑1 (MQ=255)
                                                             cTAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCAc               >  1:687968/1‑81 (MQ=255)
                                                              tAACCATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACt              >  1:687952/1‑81 (MQ=255)
                                                                 ccATGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTggg           <  2:687899/81‑1 (MQ=255)
                                                                   aTGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGtt         <  2:687900/81‑1 (MQ=255)
                                                                   aTGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGtt         <  2:687901/81‑1 (MQ=255)
                                                                   aTGACTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGt          >  2:687956/1‑80 (MQ=255)
                                                                       cTAACTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATcc     >  2:687966/1‑81 (MQ=255)
                                                                         aaCTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCctc   >  2:688002/1‑81 (MQ=255)
                                                                          aCTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCctct  >  1:687977/1‑81 (MQ=255)
                                                                          aCTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCctct  >  2:687969/1‑81 (MQ=255)
                                                                          aCTCCTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCct    <  1:687902/79‑1 (MQ=255)
                                                                               |                                                                            
CAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCT  >  NC_000913/3569260‑3569415

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: