Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,333,763 +G coding (5/615 nt) serB2 → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,333,7601.G100.0% 79.1 / NA 22M1R (ATG→AGG) †serB2phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (12/10);  total (12/10)

CGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTG  >  AM260479/3333645‑3333872
                                                                                                                    |                                                                                                                
cGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATg                                                                                                                 <  1:426120/117‑1 (MQ=255)
cGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATg                                                                                                                 >  2:426120/1‑117 (MQ=255)
                 cccgcccgACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTcgac                                                                                >  2:733356/1‑134 (MQ=255)
                                                gccgccCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCgtcgtc                                                                                   <  2:809728/100‑1 (MQ=255)
                                                gccgccCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCgtcgtc                                                                                   >  1:809728/1‑100 (MQ=255)
                                                                        tGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCg                            <  1:733356/131‑1 (MQ=255)
                                                                                  cacaCCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCt                                                                                                      >  1:804992/1‑47 (MQ=39)
                                                                                  cacaCCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCt                                                                                                      <  2:804992/47‑1 (MQ=39)
                                                                                    cacCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGaa                                                                          <  1:732680/73‑1 (MQ=255)
                                                                                    cacCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGaa                                                                          >  2:732680/1‑73 (MQ=255)
                                                                                             cccGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGt                                                 >  1:438393/1‑89 (MQ=255)
                                                                                             cccGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGt                                                 <  2:438393/89‑1 (MQ=255)
                                                                                             cccGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCg                            >  1:706292/1‑110 (MQ=255)
                                                                                             cccGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCg                            <  2:706292/110‑1 (MQ=255)
                                                                                               cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCag                                                              >  1:469510/1‑74 (MQ=255)
                                                                                               cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCag                                                              <  2:469510/74‑1 (MQ=255)
                                                                                               cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAg                               >  1:968952/1‑105 (MQ=255)
                                                                                               cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAg                               <  2:968952/105‑1 (MQ=255)
                                                                                               cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTg  >  1:1008780/1‑134 (MQ=255)
                                                                                               cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTg  >  2:888567/1‑134 (MQ=255)
                                                                                                     acAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAg                                                  <  1:991928/80‑1 (MQ=255)
                                                                                                     acAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAg                                                  >  2:991928/1‑80 (MQ=255)
                                                                                                               tGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCg                                                                  >  2:109154/1‑54 (MQ=39)
                                                                                                               tGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCg                                                                  <  1:109154/54‑1 (MQ=39)
                                                                                                                    |                                                                                                                
CGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTG  >  AM260479/3333645‑3333872

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