Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,963,552 C→T G226G (GGG→GGA)  gpp ← guanosine‑5'‑triphosphate,3'‑diphosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,963,5520CT100.0% 65.8 / NA 24G226G (GGG→GGA) gppguanosine‑5'‑triphosphate,3'‑diphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/11);  total (13/11)

GCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTCCAGCCGTGATACCGTAATTCATCGGCAACCGGACGTAACACTTCGCGTGCCGCTTTTTCTGCAGCATCAAAATTTTCCTGCCCCAGATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTT  >  NC_000913/3963361‑3963771
                                                                                                                                                                                               |                                                                                                                                                                                                                           
gcACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGcacan                                                                                                                                                                  <  1:308563/251‑2 (MQ=255)
          ccAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTCCAGCCg                                                                                                                                                        >  1:383173/1‑251 (MQ=255)
              gnTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGAGCCGGAAGCACCAACGCACACTTTCCAGCCGTGAt                                                                                                                                                    >  1:396078/1‑251 (MQ=255)
                                                       ggATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTCCAGCCGTGATACCGTAATTCATCGGCAACCGGACGTAACACTTCGCGTGcc                                                                                                           >  1:423415/1‑251 (MQ=255)
                                                        gATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGACGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTCCAGCCGTGATACCGTAATTCATCGGCAACCGGACGTAaca                                                                                                                      >  1:37415/1‑239 (MQ=255)
                                                           cGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGc                                                                                                                                                                                                                     <  1:209143/141‑1 (MQ=255)
                                                                    ccNNGGGANCCCTAACGCACGTTCCAGCGTCAGNCCGTCAATCTCCAGTTCCTCCAGCCGACCGCNANGAAGGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTNATCCATTCCCNGTGNNATCANGNTNTCCTGGAACGCCTGCACGGTGCCGGAAGCNCCAACGCACACTTTCCAGCCGTGNNACNNNAANNCATCg                                                                                                                                      <  2:383173/210‑1 (MQ=255)
                                                                         ggAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGAAACTTTTCCAGGGTACTGCGTTCATCCATTCCCTGTGc                                                                                                                                                                                                                     >  1:110374/1‑127 (MQ=255)
                                                                         ggAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCNCCAGTTCCTNNAGNCGACCNCNNNGNNNGGNTCGCTGTTTCAATTGCTGCNACNNTTCCAGNGTAATGCGTTCATCCATTCCCTGTGc                                                                                                                                                                                                                     <  2:110374/127‑1 (MQ=255)
                                                                                 naCGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTc                                                                                                                                                              >  1:126042/2‑174 (MQ=255)
                                                                                   cGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTTTAACGCCTGCACGGTGCCGTAAGCACCAACGCACACGTTCCAGCCGTGATACCTTAATTCATCGGCAACCGGACGGAACACTTCGCGTGCCGCTTTTTCTGCAGCATCAAAATTTTc                                                                                 <  1:292392/249‑1 (MQ=255)
                                                                                     cACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTCCAGCCGTGATACCGTAATTCATCGGCAAcc                                                                                                                                >  1:141699/1‑200 (MQ=255)
                                                                                        gTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCAGCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGa                                                                                                                                                                                 <  1:217290/148‑1 (MQ=255)
                                                                                         ttCCAGCGTCAGCACGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATAGCTCGCTTTTTCCCTTGCTGCACCTCTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACg                                                                                                                                                                                         >  1:189920/1‑139 (MQ=255)
                                                                                                  cAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTCCAGCCGTGATACCGTAAt                                                                                                                                            >  1:468594/1‑175 (MQ=255)
                                                                                                              ctcCAGTTCCGCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGGAACGCCTGCACGGTGCCGGAAGCACCAACGc                                                                                                                                                                      <  1:196611/137‑1 (MQ=255)
                                                                                                                       ctccAACCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTCCAGCCGTGATACCGTAATTCATCGGCAACCGGACGTAACACTTCGCGTGCCGCTTTTTCTGCAGCATCAAAATTTTcc                                                                                >  1:116234/1‑214 (MQ=255)
                                                                                                                                      atgaatgGNTNNNTNNTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCNTGTGCCATCATGATTNCCNGTAANGCCNGCACGGTGCCGGAAGCACCAACg                                                                                                                                                                       >  2:362870/1‑112 (MQ=255)
                                                                                                                                      atgaaggGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACg                                                                                                                                                                       <  1:362870/112‑1 (MQ=255)
                                                                                                                                             gCNCGCNTNNTNNATTNCTGCAACTTTTCCAGGGTAATNNGTTCATCCANTCCCTGTTCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAANGCACACGTTNNAGNCGTGANANNNTNNNTCNTCGGCAACCGGACGTAACANTTNNCGTGCCNCTTTTTCTGCAGCATCAAAATTTTcc                                                                                <  2:116234/192‑1 (MQ=255)
                                                                                                                                                      ttCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTCCAGCCGTGATACCGTAATTCATCGGCAACCGGACGTAACACTTCGCGTGCCGCTTTTTCTGCAGCATc                                                                                          <  1:151341/173‑1 (MQ=255)
                                                                                                                                                             gctgcAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCCACGCACACTTTCCAGCCGTGATACCGTAATTCATCGGCAACCGGACGTAACACTTCGCGTGCCGCTTTTTCTGCAGCATCAAAATTTTCCTGCCCCAGATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAggtg     >  1:198417/1‑251 (MQ=255)
                                                                                                                                                              ctgcAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCGGTAACGCCTGCTCGTTGTCGGAACCACCAACGCACACTTGCCAGCCGTGAGTCCGGAAGTCATCGGTAACCGGACGTAACACTTCGCGTg                                                                                                             <  1:181143/146‑1 (MQ=255)
                                                                                                                                                                gcAACTTTTCCAGGGTAATGCGTTCATCCATTCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGCTGCCGGAAACACCAACGCACACTTTCCAGCCGTGATACCGTAATTCATCGGCAACCGGACGTAACACTTCGCGTGCCGCTTTTTCTGCAGCATCAAACTTTTCCTGCCCCAGATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGtt  >  1:227925/1‑251 (MQ=255)
                                                                                                                                                                                               |                                                                                                                                                                                                                           
GCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAGCACCAACGCACACTTTCCAGCCGTGATACCGTAATTCATCGGCAACCGGACGTAACACTTCGCGTGCCGCTTTTTCTGCAGCATCAAAATTTTCCTGCCCCAGATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTT  >  NC_000913/3963361‑3963771

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