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% 68.9 / NA 19M1R (ATG→AGG) †serB2phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (11/8);  total (11/8)

CCAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTG  >  AM260479/3333633‑3333872
                                                                                                                                |                                                                                                                
ccAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATggggcc                                                                                                            >  1:51486/1‑131 (MQ=255)
       ggCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCt                                                                                                      >  2:747391/1‑134 (MQ=255)
              gcgcCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCAACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGc                                                                                               >  1:142517/1‑134 (MQ=255)
                                                tAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCag                                                              >  1:170993/1‑133 (MQ=255)
                                                            gccgccCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAg                                                  <  2:51486/133‑1 (MQ=255)
                                                                cccATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTa                                             <  1:747391/134‑1 (MQ=255)
                                                                                      aCCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCAcc                       >  2:692099/1‑134 (MQ=255)
                                                                                       ccTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCa                         >  1:847867/1‑131 (MQ=255)
                                                                                       ccTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCa                         <  2:847867/131‑1 (MQ=255)
                                                                                              cacaCCCTGATCCCGACCGACAGCGACCATGAATgg                                                                                                                >  2:915807/1‑36 (MQ=255)
                                                                                              cacaCCCTGATCCCGACCGACAGCGACCATGAATgg                                                                                                                <  1:915807/36‑1 (MQ=255)
                                                                                               acacCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATc              <  2:170993/134‑1 (MQ=255)
                                                                                                cacCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACaaa                                <  1:864554/115‑1 (MQ=255)
                                                                                                cacCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACaaa                                >  2:864554/1‑115 (MQ=255)
                                                                                                cacCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATc              >  1:1018918/1‑133 (MQ=255)
                                                                                                cacCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATc              >  1:867749/1‑133 (MQ=255)
                                                                                                cacCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATc              <  2:1018918/133‑1 (MQ=255)
                                                                                                           cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTa                                             <  2:94043/91‑1 (MQ=255)
                                                                                                           cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTa                                             >  1:94043/1‑91 (MQ=255)
                                                                                                           cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTg  >  1:228328/1‑134 (MQ=255)
                                                                                                                 acAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAg                                                  <  2:230852/80‑1 (MQ=255)
                                                                                                                 acAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAg                                                  >  1:230852/1‑80 (MQ=255)
                                                                                                                                |                                                                                                                
CCAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTG  >  AM260479/3333633‑3333872

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