Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,870,339 G→A 37.7% G134S (GGC→AGC)  iap → aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,870,3390GA37.7% 37.4 / 37.0 50G134S (GGC→AGC) iapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (19/12);  new base A (12/7);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGACCAGCGGCGAAGAGGAAGGGAAATTAGG  >  CP009273/2870202‑2870481
                                                                                                                                         |                                                                                                                                              
ttCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACg                                                                                                                                     <  1:543878/148‑1 (MQ=255)
             gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATgg                                                                                                                        <  1:597739/149‑1 (MQ=255)
                   cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATGATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATa                                                                                                                  >  2:795531/1‑149 (MQ=255)
                       aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAAcgc                                                                                                              <  2:225005/149‑1 (MQ=255)
                         gggAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCACGCTGAGCGATGCTGACGCCGATGCCGATCTCAGCGGGCTGAGGATACAAGGAAGGGATGAAAAcgccg                                                                                                            >  1:333321/1‑149 (MQ=255)
                         gggAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGACGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTAAAAGGAATGGATGATAAagccg                                                                                                            >  1:2733/1‑149 (MQ=255)
                               tACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCCCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGAAGGtt                                                                                                      >  1:273511/1‑149 (MQ=255)
                                   ggtaTTGCCGCTCACGAAGGCAAATCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAgg                                                                                                  <  2:291659/146‑1 (MQ=255)
                                    tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATGATCATTATGGCGCATCTGGATACTTGCGCCCCGCTGAGCGATGCTGACGCAGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgt                                                                                                 >  2:562593/1‑149 (MQ=255)
                                          ccGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                                                                           >  2:780472/1‑149 (MQ=255)
                                          ccGCTCATGAAGGCAAAGCGCAGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTTGACGGGCTGACGTTACAAGGACTGGATGATAACGCCGCAGGTTTAGGTGTCATg                                                                                           >  1:480393/1‑149 (MQ=255)
                                                 tGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAAt                                                                                    <  1:665836/149‑1 (MQ=255)
                                                       cAAGGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAg                                                                              <  2:415293/149‑1 (MQ=255)
                                                       cAAAGCGCCGCAGCAGATCAACATTATGGCGTATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAg                                                                              <  2:847989/149‑1 (MQ=255)
                                                          aGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAAc                                                                           >  2:53894/1‑149 (MQ=255)
                                                          aGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCACGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAGTTGGCAGAAc                                                                           >  2:855329/1‑149 (MQ=255)
                                                              ccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCAAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCt                                                                       >  1:333167/1‑149 (MQ=255)
                                                               cgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGGTTAGGTGTCATGATGGACTTGGCAGAACGccgg                                                                      >  1:352442/1‑147 (MQ=255)
                                                                 cagcagATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaa                                                                    >  1:513950/1‑149 (MQ=255)
                                                                 cagcagATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaa                                                                    >  1:256613/1‑149 (MQ=255)
                                                                 cagcagATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCGGaa                                                                    >  1:852309/1‑149 (MQ=255)
                                                                 cagcagATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGCTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaa                                                                    >  1:492880/1‑149 (MQ=255)
                                                                       atcatcATTATGGCGCAGCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGACGATAACGCCGAAGGATTAAGTGTAAGGCTGGACGTGGGCGAAGGCCAGaaaaaggg                                                              >  1:677886/1‑146 (MQ=255)
                                                                        tcataaTTATCGCGCATCTGGATACTTCCGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACg                                                             <  1:783324/149‑1 (MQ=255)
                                                                        ccctcatTTTGGCGCACCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGGTACAAGGAATGGATGATAACGCCGCAGGTGTAGGTGACCGGCTGGAACTGGCAGAACGCCTGAAAAATACg                                                             <  1:405192/146‑1 (MQ=255)
                                                                          atcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGcc                                                           >  2:294712/1‑149 (MQ=255)
                                                                              ttATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGCCGTGACAAGGAATGGATGATAACGCCGGAGGTTTAGGGTTCGTGCTGGAATTGGCAGCACGACTGAAAAATACGCCTaag                                                       >  1:240303/1‑147 (MQ=255)
                                                                                 gggCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATACCGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAGAATACGCCTACCGAg                                                    <  1:484545/148‑1 (MQ=255)
                                                                                   gcgcATCTGGATACTTACGCCCCGCTGAGTGATGCTGCCGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGACGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTa                                                  >  2:564233/1‑149 (MQ=255)
                                                                                   gcgcATCTGGATACTTACGCCCCGCTGAGGGATGCTGACGCGGATGCCAATCTCAGCGGGCTGACGTTACAAGGGATGGATGATAACGCCGCAGGTTTAGGTGTGATGATGGAATTGGGAGAGGGTCTGAAAAGTGTGCCTACTGAGTa                                                  >  1:731980/1‑149 (MQ=255)
                                                                                   gcgcATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACCAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTCCCGAgaa                                                  >  1:353342/1‑147 (MQ=255)
                                                                                   gcgcATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAACATACGCCTACCGAgaa                                                  >  1:470955/1‑147 (MQ=255)
                                                                                     gcATCTGGATAGTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTAATGCTGGCATTGGCAGAACGCCCGAAAAAGACGCCTACCGAGTagt                                                >  1:362126/1‑147 (MQ=255)
                                                                                     gcATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCGGAAAAATACGCCTACCGAGTATg                                                >  1:718984/1‑149 (MQ=255)
                                                                                       aTCTGGATACTTACGCCCCGCTGGGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGt                                              <  1:860456/149‑1 (MQ=255)
                                                                                        tCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTa                                             >  2:700552/1‑149 (MQ=255)
                                                                                             gtACTTACGCCCCGCTGAGCGATGCTGACGCCGCTGCCAATCGCAGCGGGCTGACGTTACAAGGAATGGATGCTAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCGACCGAGTATGGTATTCGa                                        <  1:121790/148‑1 (MQ=255)
                                                                                             aTACTTACGCCCCGCTGAGCGATGCGGACGCCGATGCCAATCGCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGa                                        <  2:78795/149‑1 (MQ=255)
                                                                                                cTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGCATGGGATTCGAttt                                     >  1:855228/1‑149 (MQ=255)
                                                                                                 ttACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGGATGGTATTCGATTtg                                    >  2:370336/1‑149 (MQ=255)
                                                                                                   tcGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTtgtg                                  <  2:39477/148‑1 (MQ=255)
                                                                                                       cccGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGa                              >  1:721726/1‑149 (MQ=255)
                                                                                                       cccGCGGAGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCTTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGa                              <  1:186112/149‑1 (MQ=255)
                                                                                                              aGCGATGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGGAGGTTTAGGTTTTATGCGGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGACCAgcgg                       >  1:357267/1‑149 (MQ=255)
                                                                                                                  aggcgggcgtcgcggCCAATCGCGGCGGGCTGAGGTTACAAGGAATGGATGGTAACGCCGCAGGTTTAGGTGTCAGGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGACCAGCGGCGaa                   <  1:759775/135‑1 (MQ=255)
                                                                                                                  aTGCTGACGCCGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGACCAGCGGCGaa                   <  2:505479/149‑1 (MQ=255)
                                                                                                                            cGATGCCAATCTCAGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGACCAGCGGCGAAGAGGAAGGGa         >  2:332517/1‑149 (MQ=255)
                                                                                                                              aTGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGACCAGCGGCGAAGAGGAAGGGaaa       <  1:562593/149‑1 (MQ=255)
                                                                                                                                gCCAATCTCGGCGGGCTGACGTTACAAGGAATGGAAGTAAATGCGGCGGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACGGAGGATGGTATTCGATTTGTGGCGACCAGCGGCGAAGAGGAACGGgaagt     >  1:365493/1‑147 (MQ=255)
                                                                                                                                   aaTCTCGGCGGGCTGACGTTACAAGGATTGGATGATAACGCCGCAGGTTTAGGTGTCAGGCTGGAATTGGCAGACCGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGACCAGCGGCGAAGAGGAAGGGAAATTAgg  <  2:767382/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                              
ATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGTATTCGATTTGTGGCGACCAGCGGCGAAGAGGAAGGGAAATTAGG  >  CP009273/2870202‑2870481

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