Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,455,360 G→A 13.1% A84V (GCG→GTG)  fbp ← fructose‑1,6‑bisphosphatase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,455,3600GA13.1% 54.7 / 3.4 31A84V (GCG→GTG) fbpfructose‑1,6‑bisphosphatase I
Reads supporting (aligned to +/- strand):  ref base G (16/11);  new base A (2/2);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

GTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATG  >  NC_000913/4455218‑4455502
                                                                                                                                              |                                                                                                                                              
gTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGa                                                                                                                                          <  1:107276/149‑1 (MQ=255)
    aTGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCCACAATATc                                                                                                                                      >  2:115106/1‑149 (MQ=255)
            aaTGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCATAGGCAATGCCCGCAACGATATCGCGTGcct                                                                                                                              >  1:242782/1‑147 (MQ=255)
                 tACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCCCAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGATGCAATGCCCGCAACGAGATCACGTGCTTTCAGt                                                                                                                         >  1:715670/1‑149 (MQ=255)
                     gACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCgg                                                                                                                     >  2:155951/1‑149 (MQ=255)
                                         tGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCACAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCg                                                                                                 >  2:494204/1‑149 (MQ=255)
                                         tGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCACAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCg                                                                                                 >  2:46395/1‑149 (MQ=255)
                                                    ccATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAg                                                                                      <  2:97160/149‑1 (MQ=255)
                                                                     cAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACGATCTCATATTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCCGTGCGGCTGTCAGTTTTTGATTCGCGCACAAGATGCGTATCTGCTGACCCTctc                                                                     >  2:656401/1‑147 (MQ=255)
                                                                      aGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTGTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAAcgcggcc                                                                    >  2:104557/1‑143 (MQ=255)
                                                                       gCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGCGCGGCTTACAGTTTTTCATTAGCGAACAAGTCGAGTATCTGCTGAACCTCcccc                                                                   >  1:521853/1‑149 (MQ=255)
                                                                                    tGCCTGGTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGGAACGATATCGCGTGCTTTCAGTGCGGCTTTCCGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAg                                                      <  1:174300/149‑1 (MQ=255)
                                                                                     gCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGCGCTTTCATCGCCGCTTTCAGCTATACAATAACGAACAACTCGAGATACTAACTATCCTCGCCATGCAAGTTCTCAGc                                                     >  1:253106/1‑149 (MQ=255)
                                                                                       gtgtTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCACAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCAc                                                   <  1:494204/149‑1 (MQ=255)
                                                                                       gtgtTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCACAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCAc                                                   <  1:46395/149‑1 (MQ=255)
                                                                                            cacaGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTGTCAGTTCGGCGATCCGTTTTTCACCAGCGAACAAGTCAAGTTTCTTCTGATCCTCCCCCTGCACGTTCTCTGCCCCGCTg                                              >  1:409061/1‑149 (MQ=255)
                                                                                                       aaGCCGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACAATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGat                                   <  1:144569/149‑1 (MQ=255)
                                                                                                       aaGACGAATATCCCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCGCGTTTTCAGCACCGCTGGCACCCAGGat                                   <  1:523143/149‑1 (MQ=255)
                                                                                                               aaTCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCTGTGCGGCTTTCAGCTTTTCATTAGCCAACAAGTCGAGTTTCTGCTTAACCTCGCCCTGCACGTTCTCAGCATCGCTGGCACGCAGGATATCAACCa                           >  1:443160/1‑149 (MQ=255)
                                                                                                                  ctcATCTTCTTCCTCAGAGGCAATGCCCGCAACGATCTCGCGTCCTTCCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCCCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTc                        <  2:211241/149‑1 (MQ=255)
                                                                                                                   tcAGCTTCTTCTTCAGAGGCAATGCCCGGGGGGATATCGCGTGCTTTCAGTGTGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTcc                       <  2:129931/149‑1 (MQ=255)
                                                                                                                     attttcttctccAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTg                     <  1:115106/146‑1 (MQ=255)
                                                                                                                      tcttcttctTCAGAGGCAATGTACGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGc                    >  2:775685/1‑149 (MQ=255)
                                                                                                                           ttcttcAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACAGCGGGCACCCAGGATATCAACCAGTCCTGCTttgt               >  2:437853/1‑149 (MQ=255)
                                                                                                                                 agagGCAATGCCCGCAACGTTCTCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGATCCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGatat         <  2:148020/149‑1 (MQ=255)
                                                                                                                                 agagGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGatat         <  1:707335/149‑1 (MQ=255)
                                                                                                                                 agagGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACATCGCCCTGCACGTTCTCAGCTCCGCTGGCACCCAGGATACCAACCAGGCCTGCTttgttgtaat         >  2:411877/1‑145 (MQ=255)
                                                                                                                                     gCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAGCCAATCATGCTTTGTTGATAAcgtg     >  2:57878/1‑147 (MQ=255)
                                                                                                                                       gaTGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGAt   <  1:60506/148‑1 (MQ=255)
                                                                                                                                        aTGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATg  >  2:128327/1‑149 (MQ=255)
                                                                                                                                        aTGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATg  >  1:630337/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATG  >  NC_000913/4455218‑4455502

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