Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 721,950 A→G intergenic (+240/‑94) unknown → / → unknown Protoporphyrinogen IX oxidase, oxygen‑independent, HemG/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1721,9500AG100.0% 62.7 / NA 22intergenic (+240/‑94)unknown/unknownProtoporphyrinogen IX oxidase, oxygen‑independent, HemG/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/5);  total (17/5)

GTAAGTTGTATTTTTTCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAAGATTTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATAATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCTT  >  1/721711‑722095
                                                                                                                                                                                                                                               |                                                                                                                                                 
gtaagtTGTATTTTTTCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGAtt                                                                                                                                        >  1:819090/1‑251 (MQ=255)
       gTATTTTTTCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGAt                                                                                                                                         >  1:433639/1‑243 (MQ=255)
            ttttCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAAc                                                                                                                                    >  1:772768/1‑243 (MQ=255)
             tttCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGt                                                                                                                           >  1:259269/1‑251 (MQ=255)
                                     gTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAg                                                                                                                                           <  1:644177/211‑1 (MQ=255)
                                               gCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAgg                                                                                                                            >  1:536855/1‑216 (MQ=255)
                                                    aTCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGAt                                                                                                                                         >  1:831553/1‑198 (MQ=255)
                                                    aTCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGatt                                                                                    >  1:633800/1‑251 (MQ=255)
                                                    aTCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGatt                                                                                    >  1:642263/1‑251 (MQ=255)
                                                       cTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTGCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         <  1:450090/243‑1 (MQ=255)
                                                            aTAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTg                                                                                                                          >  1:516909/1‑205 (MQ=255)
                                                                  aTCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTGCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTa                                                                      <  1:760200/251‑1 (MQ=255)
                                                                                cATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAg                                                         >  1:760756/1‑250 (MQ=255)
                                                                                        cTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         <  1:39879/210‑1 (MQ=255)
                                                                                              gggCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAag                                          >  1:575543/1‑251 (MQ=255)
                                                                                                gCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         >  1:708114/1‑202 (MQ=255)
                                                                                                                        aaaTATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAgg                                                                                                                            >  1:835565/1‑143 (MQ=255)
                                                                                                                            atTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTgg            >  1:735344/1‑251 (MQ=255)
                                                                                                                                      cGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCtt  >  1:235485/1‑251 (MQ=255)
                                                                                                                                           aataatATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         >  1:341244/1‑159 (MQ=255)
                                                                                                                                                 atcttCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         <  1:335128/150‑1 (MQ=255)
                                                                                                                                                                             tAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTttatttatatta                                                                                                 >  1:370777/1‑117 (MQ=21)
                                                                                                                                                                                                                                               |                                                                                                                                                 
GTAAGTTGTATTTTTTCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAAGATTTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATAATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCTT  >  1/721711‑722095

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