Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 363,780 C→T A842A (GCG→GCA lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913363,7800CT100.0% 84.3 / NA 25A842A (GCG→GCAlacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/13);  total (12/13)

GGCCAATCCGCGCCGGATGCGGTGTATCGCTCGCCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCACGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGCCTGGTAATGGCCCGCCGCCTTCCAGCGTTCGAC  >  NC_000913/363653‑363884
                                                                                                                               |                                                                                                        
ggCCAATCCGCGCCGGATGCGGTGTATCGCTCGCCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTg                                                                                                     <  1:428433/133‑1 (MQ=255)
  ccAATCCGCGCCGGATGCGGTGTATCGCTCGCCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAg                                                                                                   <  1:124969/133‑1 (MQ=255)
  ccAATCCGCGCCGGATGCGGTGTATCGCTCGCCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAg                                                                                                   <  2:411610/133‑1 (MQ=255)
  ccAATCCGCGCCGGATGCGGTGTATCGCCCGCAACTTCAACTTCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAg                                                                                                   <  2:159172/133‑1 (MQ=255)
      tCCGCGCCGGATGCGGTGTATCGCTCGCCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGt                                                                                               >  2:23820/1‑133 (MQ=255)
           gcCGGATGCGGTGTATCGCTCGCCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTaa                                                                                          >  2:409272/1‑133 (MQ=255)
             cGGATGCGGTGTATCGCTCGCCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATc                                                                                        <  1:23820/133‑1 (MQ=255)
                     gtgtATCGCTCGCCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGc                                                                                <  2:391674/133‑1 (MQ=255)
                             ctcgcCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCaa                                                                        >  2:123773/1‑133 (MQ=255)
                                gcCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGATCGGTCGTTATCGGCGCCGCATCAGCAAgtg                                                                     >  1:188323/1‑133 (MQ=255)
                                gcCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCTCCGCATAAGCAAgtg                                                                     >  1:185039/1‑133 (MQ=255)
                                  cACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTa                                                                   >  1:482925/1‑133 (MQ=255)
                                           aTCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTg                                                          <  2:482925/133‑1 (MQ=255)
                                           aTCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTg                                                          <  1:426161/133‑1 (MQ=255)
                                               aCGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACt                                                      >  2:370759/1‑133 (MQ=255)
                                                             ttGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACgctgct                                        >  2:389021/1‑133 (MQ=255)
                                                              tGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCtt                                       <  1:378588/133‑1 (MQ=255)
                                                                  cACTCCCACCAATCCCGTAGGTCTTCCGGCAGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGc                                   <  1:192992/133‑1 (MQ=255)
                                                                   aCTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGcc                                  <  1:123773/133‑1 (MQ=255)
                                                                          tCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGCCTGGTAAt                           <  1:442082/133‑1 (MQ=255)
                                                                                   tAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGCCTGGTAATGGCccgccg                  >  1:41448/1‑133 (MQ=255)
                                                                                   tAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGCCTGGTAATGGCccgccg                  >  2:448174/1‑133 (MQ=255)
                                                                                      gTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGCCTGGTAATGGCCCGCCGCCt               >  1:406288/1‑133 (MQ=255)
                                                                                        tttCCGGCTGATTAATAAGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGCCTGGTAATGGCCCGCCGCCTTc             >  1:86634/1‑133 (MQ=255)
                                                                                                   taaataaGGTTTTCCCCTGATGCTGCCATGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGCCTGGTAATGGCCCGCCGCCTTCCAGCGTTCGAc  <  1:448174/133‑1 (MQ=255)
                                                                                                                               |                                                                                                        
GGCCAATCCGCGCCGGATGCGGTGTATCGCTCGCCACTTCAACATCAACGGTAATCGCCATTTGACCACTACCATCAATCCGGTAGGTTTTCCGGCTGATAAATAAGGTTTTCCCCTGATGCTGCCACGCGTGAGCGGTCGTAATCAGCACCGCATCAGCAAGTGTATCTGCCGTGCACTGCAACAACGCTGCTTCGGCCTGGTAATGGCCCGCCGCCTTCCAGCGTTCGAC  >  NC_000913/363653‑363884

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: