Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,338,986 T→C V71V (GTA→GTG)  unknown ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,338,9860TC100.0% 70.1 / NA 26V71V (GTA→GTG) unknownhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/8);  total (18/8)

ATGATACACACTTTCTCCTTTAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATAGGTTCCAACAACAAGAAACTTTACTGAATGCCCATTTTTAAGAGGGAAAATATGAAACTTCCCGTCGGTATAGGCATCAGCAATATAACTAAGTTTTTCTGATAGACACTCTATTTCCTCAATGAATTTGTCAGCTACTTCCTTATTTTGGCTTTGTGTGTAAATATATTCTGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAATGTAGTCTTGTTTCATCGCTTCAATGCCT  >  1/1338901‑1339211
                                                                                     |                                                                                                                                                                                                                                 
aaataacatacaccccttcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTc                                                                              >  1:818686/21‑235 (MQ=255)
     acacacaccccttcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTAt                                                                                                                                       <  1:448188/158‑1 (MQ=255)
                ttcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAata                                                                                   >  1:225443/5‑214 (MQ=255)
                 tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCAttt                                              >  1:604529/4‑250 (MQ=255)
                 tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCAttt                                              >  1:545586/4‑250 (MQ=255)
                 tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGATTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTaaa                                                                                     >  1:79663/4‑211 (MQ=255)
                      aTGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAcc                                        <  1:87745/251‑1 (MQ=255)
                      aTGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAcc                                        <  1:425820/251‑1 (MQ=255)
                        gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCtt                                                                                               >  1:491282/1‑194 (MQ=255)
                        gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCAt                                      >  1:93247/1‑251 (MQ=255)
                        gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCAt                                      >  1:629036/1‑251 (MQ=255)
                                ccTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTa                                                                                                                                        >  1:43191/1‑145 (MQ=255)
                                     tGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGC‑GCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTaatc                        <  1:193324/251‑4 (MQ=255)
                                                cGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCAt               <  1:657619/250‑1 (MQ=255)
                                                      tCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCa        >  1:744739/1‑251 (MQ=255)
                                                      tCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCa        >  1:400564/1‑251 (MQ=255)
                                                       cATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTc                                             >  1:792671/1‑213 (MQ=255)
                                                        ataataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTaatcc                       >  1:39444/1‑230 (MQ=255)
                                                        ataataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAAt      >  1:491142/1‑251 (MQ=255)
                                                         taataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTccgc                                                                           <  1:269537/181‑4 (MQ=255)
                                                         taataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACaacgt                            >  1:592328/1‑225 (MQ=255)
                                                           atatacattCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGc    >  1:623153/8‑250 (MQ=255)
                                                            tatacattCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCt  <  1:665540/245‑1 (MQ=255)
                                                                    cctacaaccaGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTtgtg                                                                                           >  1:767337/4‑154 (MQ=255)
                                                                    cctacaaccaGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAata                                                                                   <  1:488859/159‑1 (MQ=255)
                                                                                 aCTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAc                                         >  1:518623/1‑191 (MQ=255)
                                                                                     |                                                                                                                                                                                                                                 
ATGATACACACTTTCTCCTTTAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATAGGTTCCAACAACAAGAAACTTTACTGAATGCCCATTTTTAAGAGGGAAAATATGAAACTTCCCGTCGGTATAGGCATCAGCAATATAACTAAGTTTTTCTGATAGACACTCTATTTCCTCAATGAATTTGTCAGCTACTTCCTTATTTTGGCTTTGTGTGTAAATATATTCTGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAATGTAGTCTTGTTTCATCGCTTCAATGCCT  >  1/1338901‑1339211

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