Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,886,731 T→C I163M (ATA→ATG)  cas3 ← CRISPR‑associated endonuclease/helicase Cas3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,886,7310TC100.0% 71.1 / NA 22I163M (ATA→ATG) cas3CRISPR‑associated endonuclease/helicase Cas3
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/9);  total (13/9)

CCTCATTAAACAGAAAGGTATTCGTTGTAGTCCAGGAGCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAAT  >  NC_000913/2886487‑2886974
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
ccTCATTAAACAGCAAGGTATTCGTTGTAGTCCAGGAGCCTAACCAGTCAGCAAGCGAGCAAAAACCTACTAACAGTGAAGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTCCAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGTTCCTGGGAATGTAACATaaaa                                                                                                                                                                                                                                               <  2:77247/251‑1 (MQ=255)
   cATTAAACAGAAAGGTATTCGTTGTAGTCCAGGAGCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAAc                                                                                                                                                                                                                                              <  1:50109/249‑1 (MQ=255)
     ttAAACAGAAAGGTATTCGTTGTAGTCCAGGAGCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGAGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGCTCCTGGGAATGTAACATAAAACCatg                                                                                                                                                                                                                                          >  1:126771/1‑251 (MQ=255)
                   aTTCGTTGTAGTCCAGGAGCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGTGAATGTAACATAAAACCATGATGTCCTGTAACgg                                                                                                                                                                                                                            <  2:162836/251‑1 (MQ=255)
                                    aGCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGCTGTCCTGTAACGGCCTCTACCCATGGAAAc                                                                                                                                                                                                           >  1:343996/1‑251 (MQ=255)
                                         aaCCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAg                                                                                                                                                                                                        <  1:72866/249‑1 (MQ=255)
                                                             aaaaCCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCgg                                                                                                                                                                                  >  1:285927/1‑251 (MQ=255)
                                                                    gCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCACGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCAACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACAATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCTTCac                                                                                                                                                                            >  2:227550/1‑249 (MQ=255)
                                                                            tgatgaACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTg                                                                                                                                                                    >  1:350067/1‑250 (MQ=255)
                                                                                                               aCCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATcc                                                                                                                                >  1:293380/1‑251 (MQ=255)
                                                                                                                      cccctacccgtcgctcttccgatctACAGATATCCACCCCCCACGGGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTaa                                                                                                                         <  2:307764/227‑1 (MQ=255)
                                                                                                                                          ccAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGCAAGTGAATCCTGGTTAACCCAATACAGACCGGCTGc                                                                                                      >  2:156151/1‑250 (MQ=255)
                                                                                                                                              tACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACt                                                                                                                                                                                                    >  2:207772/1‑152 (MQ=255)
                                                                                                                                              tACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACt                                                                                                                                                                                                    <  1:207772/152‑1 (MQ=255)
                                                                                                                                               aCAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTaa                                                                                                                         >  1:307764/1‑226 (MQ=255)
                                                                                                                                                  gATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGAtt                                                                                             >  1:229977/1‑251 (MQ=255)
                                                                                                                                                        cACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAAc                                                                                                                                                                                                                              <  1:2437/116‑1 (MQ=255)
                                                                                                                                                        cACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAAc                                                                                                                                                                                                                              >  2:2437/1‑116 (MQ=255)
                                                                                                                                                                                    cagcaaaaagttccAGAGAAGCCGGCATTTCCCAACGCTACGTAC‑CTGGATCCTGGGAATGTAACATAAAACCATGATGGCCGGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGTGGCATCAAAAAAACTGAAAAAATCCCCGAGAGAGTTCTCAGAAAGGGAAGCCTGGTTAAACCAATACGGACCGGCTGCACCATGATGAAATCGACGGCACATTTGTGTCGATGGACCAt                                                             <  2:68441/237‑1 (MQ=255)
                                                                                                                                                                                          aaGATGCAATAGAATCGGGCATTTCCCAACGCTA‑TTGTACCTGATCCTGGGAATGTAACATAAAACCATGGTGTCCTGTAACGGCCTCTACCCATGGAAACCATGACTCATAAGGATGAGGAGCGTCAGTAAAAAAACTGAATAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACTAATACAGACCGGCTGCAACATGATTAAATTTACGGCACATTTGTGTTCATGGACCATTAAGTGAt                                                     <  2:466085/251‑1 (MQ=255)
                                                                                                                                                                                                                 tCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGAGTTGTAGGATTTAATTTc                                >  1:265284/1‑249 (MQ=255)
                                                                                                                                                                                                                                              aaTGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAAt  >  2:232983/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
CCTCATTAAACAGAAAGGTATTCGTTGTAGTCCAGGAGCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTAT‑CTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAAT  >  NC_000913/2886487‑2886974

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