Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 910,204 G→A V5I (GTA→ATA)  unknown → phage‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1910,2040GA100.0% 80.5 / NA 32V5I (GTA→ATA) unknownphage‑related protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/14);  total (18/14)

GGTATTTCCGGCAAGTCACATACTCACCGGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAGTAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACCGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATCTGGAACATGGCTTGCCGTGGCTTGAACAAATGGGGCGCGGTGTGAATCTGGCGGATTGGGAAATTAAAATAAAGCGGCATGTACTGGGAACGGAAGGAGTGTCAAAAA  >  1/910130‑910444
                                                                          |                                                                                                                                                                                                                                                
gaaaaatacgttaaGTCACATACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCt                                                                                                                                                                              <  1:757406/131‑1 (MQ=255)
    aatacgttaaGTCACATACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATc                                                              <  1:313048/243‑1 (MQ=255)
    aagacgttaaaccacaTACACACCCGGGCGACAGCGGCGGAACAACAGGAACACCCCAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGCCGTTTGGGCAAGGTTGTTCGAACTACGCAACTGAATCTGAAGCTATCGCCCAGAATGTACAACCCCGCCTTTGGTCGTTTGCCAATGGCTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATc                                                              <  1:674304/239‑1 (MQ=255)
              gTCACATACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTg                                                                                                                                                                                                              <  1:189321/97‑1 (MQ=255)
              gTCACATACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCGATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTg                                                    >  1:576381/1‑251 (MQ=255)
              gTCACATACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGtt                                                                                                                      >  1:468856/1‑185 (MQ=255)
              gTCACATACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGagc                                                                                   >  1:801992/1‑218 (MQ=255)
              gTCACATACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTg                                                    >  1:168851/1‑251 (MQ=255)
              gTCACATACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTg                                                    >  1:217588/1‑251 (MQ=255)
              gTCACATACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTg                                                    >  1:795527/1‑251 (MQ=255)
                cacaTACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAAt                                                                                                                                                                                <  1:667333/125‑1 (MQ=255)
                cacaTACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTAc                                                                                                                                                        <  1:139432/149‑1 (MQ=255)
                cacaTACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTggg                                                  >  1:139580/1‑251 (MQ=255)
                  caTACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATc                                                              >  1:199358/1‑237 (MQ=255)
                  caTACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGaa                                                >  1:294071/1‑251 (MQ=255)
                  caTACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGaa                                                >  1:123428/1‑251 (MQ=255)
                   aTACTCACCCGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGAtttagaac                                                                                                        <  1:575264/194‑8 (MQ=255)
                             gggCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAAttaa                                     >  1:603440/1‑248 (MQ=255)
                              ggCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAAttaaa                                    >  1:868683/1‑247 (MQ=255)
                                         cggAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgt                          <  1:391640/250‑3 (MQ=255)
                                           gAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGAt                                                                                                                                                                                                                        >  1:4864/1‑58 (MQ=255)
                                              caacaGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAAt                                                               <  1:754933/208‑1 (MQ=255)
                                               aacaGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagac                   <  1:528246/251‑10 (MQ=255)
                                               aacaGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagac                   <  1:538397/251‑10 (MQ=255)
                                                acaGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggaga                    >  1:301100/1‑241 (MQ=255)
                                                 caGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggaga                    >  1:771209/1‑240 (MQ=255)
                                                    gAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTg                                                    >  1:785865/1‑213 (MQ=255)
                                                    gAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagactgaag              >  1:880160/1‑237 (MQ=255)
                                                     aaCACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagactgaagg             <  1:484478/251‑16 (MQ=255)
                                                     aaCACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagactgaagg             <  1:381610/251‑16 (MQ=255)
                                                         ccacAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACa                                                                                                                                                       <  1:870998/109‑1 (MQ=255)
                                                                gACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCaaaaa  >  1:400948/1‑251 (MQ=255)
                                                                          |                                                                                                                                                                                                                                                
GGTATTTCCGGCAAGTCACATACTCACCGGGGCGACAGCGGCGGAACAACAGGAACACCACAATGACAACAAGAGTAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACCGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATCTGGAACATGGCTTGCCGTGGCTTGAACAAATGGGGCGCGGTGTGAATCTGGCGGATTGGGAAATTAAAATAAAGCGGCATGTACTGGGAACGGAAGGAGTGTCAAAAA  >  1/910130‑910444

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