Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 858,223 T→C L182P (CTG→CCG)  unknown → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1858,2230TC100.0% 43.5 / NA 21L182P (CTG→CCG) unknownhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/3);  total (18/3)

TCATCTTCTCAAAACTAAAATAATCAAAGTAATCGAAACATTTGACGGTTATAAGGATTTAAACTGCTTAGCCAGTTGGCTATTATTTAGCGGTCAGCAAGTCGCAACACTTAGTGATCTTTATCATAAAGATTTTTGGAATTGTATTAGAAAATCAGATTATCCAAATGCCCATGGTTATTTGGATGGTATAGAGATTACCAATGAATCATTCCACACGTTACTTCTGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCAATAATTGGCAGGCATTTGATGGTGCACAAAGAATTA  >  1/857996‑858465
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
tcatcTTCTCAAAACTAAAATAATCAAAGTAATCGAAACATTTGACGGTTATAAGGATTTAAACTGCTTAGCCAGTTGGCTATTATTTAGCGGTCAGCAAGTCGCAACACTTAGTGATCTTTATCATAAAGATTTTTGGAATTGTATTAGAAAATCAGATTATCCAAATGCCCATGGTTATTTGGATGGTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAggg                                                                                                                                                                                                                             >  1:258403/1‑251 (MQ=255)
           aaaCTAAAATAATCAAAGTAATCGAAACATTTGACGGTTATAAGGATTTAAACTGCTTAGCCAGTTGGCTATTATTTAGCGGTCAGCAAGTCGCAACACTTAGTGATCTTTATCATAAAGATTTTTGGAATTGTATTAGAAAATCAGATTATCCAAATGCCCATGGTTATTTGGATGGTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTtgtg                                                                                                                                                                                                                  <  1:396270/251‑1 (MQ=255)
           aaaCTAAAATAATCAAAGTAATCGAAACATTTGACGGTTATAAGGATTTAAACTGCTTAGCCAGTTGGCTATTATTTAGCGGTCAGCAAGTCGCAACACTTAGTGATCTTTATCATAAAGATTTTTGGAATTGTATTAGAAAATCAGATTATCCAAATGCCCATGGTTATTTGGATGGTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTtgtg                                                                                                                                                                                                                  >  1:551408/1‑251 (MQ=255)
                    taatCAAAGTAATCGAAACATTTGACGGTTATAAGGATTTAAACTGCTTAGCCAGTTGGCTATTATTTAGCGGTCAGCAAGTCGCAACACTTAGTGATCTTTATCATAAAGATTTTTGGAATTGTATTAGAAAATCAGATTATCCAAATGCCCATGGTTATTTGGATGGTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTtg                                                                                                                                                                                                                    <  1:785725/240‑1 (MQ=255)
                                        tttGACGGTTATAAGGATTTAAACTGCTTAGCCAGTTGGCTATTATTTAGCGGTCAGCAAGTCGCAACACTTAGTGATCTTTATCATAAAGATTTTTGGAATTGTATTAGAAAATCAGATTATCCAAATGCCCATGGTTATTTGGATGGTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCa                                                                                                                                                                                                                                         <  1:547277/199‑1 (MQ=14)
                                                                                                                                                                              tGGTTATTTGGATGGTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAgg                                               >  1:362020/1‑251 (MQ=11)
                                                                                                                                                                                            gTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAg                                                >  1:209066/1‑236 (MQ=11)
                                                                                                                                                                                            gTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCaataat                                 >  1:240976/1‑251 (MQ=14)
                                                                                                                                                                                            gTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCaataat                                 >  1:38697/1‑251 (MQ=14)
                                                                                                                                                                                            gTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCaataat                                 >  1:174486/1‑251 (MQ=14)
                                                                                                                                                                                            gTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCaataat                                 >  1:860379/1‑251 (MQ=14)
                                                                                                                                                                                            gTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCaataat                                 >  1:512813/1‑251 (MQ=14)
                                                                                                                                                                                            gTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCaataat                                 >  1:542332/1‑251 (MQ=14)
                                                                                                                                                                                            gTATAGAGATTACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCaataat                                 >  1:750506/1‑251 (MQ=14)
                                                                                                                                                                                                     ttACCAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCAATAATTggcaggca                        >  1:870993/1‑251 (MQ=18)
                                                                                                                                                                                                        ccAATGAATCATTCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCAATAATTGGCAGGCAttt                     >  1:394526/1‑251 (MQ=18)
                                                                                                                                                                                                                   ttCCACACGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCaataat                                 >  1:826012/1‑228 (MQ=17)
                                                                                                                                                                                                                          cGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCAATAATTGGCAGGCATTTGATGGTGCACAAAGAAtt   >  1:749754/1‑251 (MQ=37)
                                                                                                                                                                                                                          cGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCAATAATTGGCAGGCATTTGATGGTGCACAAAGAAtt   >  1:85443/1‑251 (MQ=37)
                                                                                                                                                                                                                          cGTTACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCAATAATTGGCAGGCATTTGATGGTGCACAAAGAAtt   >  1:406892/1‑251 (MQ=37)
                                                                                                                                                                                                                            ttACTTCCGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCAATAATTGGCAGGCATTTGATGGTGCACAAAGAATTa  >  1:527464/1‑250 (MQ=37)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
TCATCTTCTCAAAACTAAAATAATCAAAGTAATCGAAACATTTGACGGTTATAAGGATTTAAACTGCTTAGCCAGTTGGCTATTATTTAGCGGTCAGCAAGTCGCAACACTTAGTGATCTTTATCATAAAGATTTTTGGAATTGTATTAGAAAATCAGATTATCCAAATGCCCATGGTTATTTGGATGGTATAGAGATTACCAATGAATCATTCCACACGTTACTTCTGCGTTTTATCAATGACGAAAGGGCGGTATTTGTGTTAGACCCGCCATATTTGTGTACGAAACAAGAGAGCTATAAGCAAGCGCATTATTTCGATCTTATCGACTTCCTACGCTTAATTAACATTACACGACCACCGTATATCTTCTTCTCATCGACAAAAAGCGAGTTCGTGCGGTTTATTGAGTATATGTGTGAGGATAAAGTCAATAATTGGCAGGCATTTGATGGTGCACAAAGAATTA  >  1/857996‑858465

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