Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,970,230 A→G L252P (CTG→CCG)  tap ← methyl‑accepting protein IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,970,2300AG100.0% 163.5 / NA 46L252P (CTG→CCG) tapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/26);  total (20/26)

ATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTGGTTGAACAATCATCTTGCG  >  NC_000913/1970096‑1970354
                                                                                                                                      |                                                                                                                            
aTGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCCACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACt                                                                                                                          <  1:916538/139‑1 (MQ=255)
             ttCGGTACGACTTGAGAGATCGTTATTACCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCg                                                                                                             <  2:166199/139‑1 (MQ=255)
                  tACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACACAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGAt                                                                                                        <  2:131162/139‑1 (MQ=255)
                   aCGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGTGATACCAATGTCCATCTCCTGTCTTCCCTTACGCACGTCAGTTACGGTACAACGGTAAGCCTGCTGGATGGTCTTCGGACTGGTAAAAATTGCGGTGAtc                                                                                                       <  2:680866/139‑1 (MQ=255)
                                 cGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGAcc                                                                                            >  1:888230/1‑136 (MQ=255)
                                 cGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGAcc                                                                                            <  2:888230/136‑1 (MQ=255)
                                 cGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCata                                                                                         <  1:631051/139‑1 (MQ=255)
                                     aTTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTAc                                                                                                <  2:378376/128‑1 (MQ=255)
                                     aTTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTAc                                                                                                >  1:378376/1‑128 (MQ=255)
                                             cGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGAt                                                                                                        <  1:204834/112‑1 (MQ=255)
                                             cGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGAt                                                                                                        >  2:204834/1‑112 (MQ=255)
                                              gACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGa                                                                            <  2:287763/139‑1 (MQ=255)
                                                            tACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTAcc                                                              <  1:386757/139‑1 (MQ=255)
                                                            tACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTAcc                                                              <  1:568427/139‑1 (MQ=255)
                                                             aCCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTAccc                                                             >  2:509709/1‑139 (MQ=255)
                                                              ccAATGTGCATCTCCTGGCTTCTCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTa                                                                                                 >  1:893144/1‑102 (MQ=255)
                                                              ccAATGTGCATCTCCTGGCTTCTCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTa                                                                                                 <  2:893144/102‑1 (MQ=255)
                                                                 aTGTGCACCTCCTGGCTTCCCATACGCCCATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCgcag                                                         <  2:390974/139‑1 (MQ=255)
                                                                   gtgCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCgcagca                                                       >  1:862978/1‑139 (MQ=255)
                                                                      cATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGAcgc                                                                         >  2:811606/1‑118 (MQ=255)
                                                                      cATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGAcgc                                                                         <  1:811606/118‑1 (MQ=255)
                                                                         ctcCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTg                                                 >  1:236605/1‑139 (MQ=255)
                                                                           ccTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTc                                               >  1:674330/1‑139 (MQ=255)
                                                                              ggCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCaaa                                            >  1:619010/1‑139 (MQ=255)
                                                                                cTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAAt                                          >  2:269312/1‑139 (MQ=255)
                                                                                    ccTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGAcgc                                                                         >  2:243592/1‑104 (MQ=255)
                                                                                    ccTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGAcgc                                                                         <  1:243592/104‑1 (MQ=255)
                                                                                              aTCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTa                                     >  2:893803/1‑130 (MQ=255)
                                                                                              aTCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTa                                     <  1:893803/130‑1 (MQ=255)
                                                                                                cACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGc                          <  1:803420/139‑1 (MQ=255)
                                                                                                 aCTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGTTGTCAAAATGGCTACCGATAATGGcc                         <  1:309970/139‑1 (MQ=255)
                                                                                                 aCTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGcc                         <  1:509709/139‑1 (MQ=255)
                                                                                                         tCCCACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTGGTTg                 <  1:361152/139‑1 (MQ=255)
                                                                                                          cccACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACg                                                                                               <  2:750438/60‑1 (MQ=255)
                                                                                                          cccACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCGCATTACg                                                                                               >  1:750438/1‑60 (MQ=255)
                                                                                                           ccACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTa                                                                                                 <  1:1113724/57‑1 (MQ=255)
                                                                                                           ccACGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTa                                                                                                 >  2:1113724/1‑57 (MQ=255)
                                                                                                             aCGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTg                     >  2:44354/1‑131 (MQ=255)
                                                                                                             aCGCAAAGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTg                     <  1:44354/131‑1 (MQ=255)
                                                                                                                   aGCCTGCTGCATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTGGTTGAACAatcatc       >  1:311231/1‑139 (MQ=255)
                                                                                                                       tgctgcATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTGGTTGAACAATCATCTTgc   <  1:614275/139‑1 (MQ=255)
                                                                                                                        gctgcATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTGGTTGAACAATCATCTTgc   <  1:951236/138‑1 (MQ=255)
                                                                                                                        gctgcATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTGGTTGAACAATCATCTTgc   >  2:951236/1‑138 (MQ=255)
                                                                                                                        gctgcATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCACAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTGGTTGAACAATCATCTTgcg  >  1:1082253/1‑139 (MQ=255)
                                                                                                                         ctgcATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGa                                                                                              <  2:301830/46‑1 (MQ=255)
                                                                                                                         ctgcATGGTCTTCGGACTGGCAAAAATGGCGGTGATCTCATTACGa                                                                                              >  1:301830/1‑46 (MQ=255)
                                                                                                                                      |                                                                                                                            
ATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAATGGCCAGTGGTTGAACAATCATCTTGCG  >  NC_000913/1970096‑1970354

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