Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 363,362 C→T G982S (GGC→AGC)  lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913363,3620CT100.0% 52.3 / NA 20G982S (GGC→AGC) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/9);  total (11/9)

AAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGC  >  NC_000913/363119‑363577
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                       
aaaaGTAAAAGAAAAAGAATAAACCGAACATCCAAAATTGTGGGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAAGACGGGGAGACATGGCCTGCCTGGTTATGATTATTTTTGACCCCAGTCCAACTGGCAATGGTAGCAACCGTGGCTCAGCTGGAATCCCGCCGATACTGCCGGGCGCCAGGAGTCGTCGCCACCAATCCGCATATGGAAACCGTCGATATTCAGCCATGTGCTTTCTTcc                                                                                                                                                                                                                  <  1:283324/251‑1 (MQ=255)
                                                                                                      ccGGTTATTATTATTTTTGACACCAGACCAACTGGGAATGGTAGCGACCCGCGCTCAGCTGGGATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAAACGTAGCTATTCCGCCCTGTGCTTTCTTACGCGt                                                                                                                                                                                                              >  1:324925/1‑153 (MQ=255)
                                                                                                      ccGGTTANTNTTATTTTTGACACCAGACCAACGGGTAATGGTAGCGACCGGCGCTCAGNTGGAATTCCGNCGANNCTGANGNGNTCCAGGAGTCGTCGCCACCAATCCCCATANGGAAACCGTCGATATTCAGCCATNNTCNNNCTNNCGCGt                                                                                                                                                                                                              <  2:324925/153‑1 (MQ=255)
                                                                                                                          nacCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGCAACCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGACGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAAc                                                                                                                                                >  1:53245/2‑195 (MQ=255)
                                                                                                                           acCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTc                                                                                                                                        >  1:20495/1‑202 (MQ=255)
                                                                                                                               gACCAGCTGGTAATGGTAGCGCCCGTCGCGCAGCTGGAATTCCGCCGATTCTGACGGGCTCCAGGAGTCGTCGCCTCCAATCCCCATATGGTACCCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGGTGCTGTTGCCTGTAGCGTCGGCTGTTGAGCTGGCAGCCGGCGTTCCACTGGTTTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGaa                                                                                   <  1:167197/251‑2 (MQ=255)
                                                                                                                                                aGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTcgn                                                                                                                                      <  1:91327/183‑2 (MQ=255)
                                                                                                                                                 gCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGCAACCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGCACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGAGTCCCGCAGCGCAGACCGTTTTCGCTCGGGa                                                                    >  1:294050/1‑248 (MQ=255)
                                                                                                                                                  cGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACTGTCGATATTCAGCCATGTGATTTCTTCCGTGTGCAGCAGATGGAGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCCTAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGAc                                                                >  1:291047/1‑251 (MQ=255)
                                                                                                                                                             cAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTACTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGtttn                                                                                                                                                                                     <  1:280144/123‑2 (MQ=255)
                                                                                                                                                             cAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTACTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTc                                                                                                                                                                                     >  2:280144/1‑123 (MQ=255)
                                                                                                                                                                  ggAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATCTGGAACCCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTgcagca                                                                                                                                                                                                        >  1:75915/1‑99 (MQ=255)
                                                                                                                                                                                                 gtcgtcGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTcgc                                                                                                                                      <  1:26744/134‑1 (MQ=255)
                                                                                                                                                                                                            caATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGntgtn                                                                                                                       >  2:275635/1‑137 (MQ=255)
                                                                                                                                                                                                            caATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGcn                                                                                                                   <  1:275635/142‑2 (MQ=255)
                                                                                                                                                                                                                ccccATATGGAAACCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCcgcg                                                                                                                                  >  1:383964/1‑123 (MQ=255)
                                                                                                                                                                                                                ccccATAGGGAATCCGTCTATATTCAGCCATGTGCTTTCTTCCACGAGCAGCAGAGGGCAATGGCGTGATTCCAGCAGTTGCTGTTGACTTTAGCGGCTGATTTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAAGACGCACGTCCCGCATCGCAGACCGTTACCGCTCGGCAAGACGTACGTGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGc  <  1:61842/251‑1 (MQ=255)
                                                                                                                                                                                                                 cccATATGGAAACCGTCGATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATAAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTcc                                                                                               >  1:209453/1‑157 (MQ=255)
                                                                                                                                                                                                                   cATATGGAAACCGTCGATATTCTGCCAGGTGCTTTCTTCCGCGGGCAGCATATGTCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGa                                                                                                                                                  <  1:224771/104‑1 (MQ=255)
                                                                                                                                                                                                                                  gATATTCAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTACATAAGTTGCTGTTGACTGTCGCGACTGATGTTGAACTGGAAGTCGCCCCGACACTGGTGTGGGCCATaattca                                                                                                           >  1:330081/1‑128 (MQ=255)
                                                                                                                                                                                                                                        cAGCCATGTGCTTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCTCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAAc                 <  1:383508/212‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                       
AAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGC  >  NC_000913/363119‑363577

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