Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,963,709 A→G V3A (GTA→GCA) 
*157Q (TAA→CAA) 
ppdB ←
ppdA ←
conserved protein PpdB
conserved protein PpdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,963,7090AG100.0% 99.4 / NA 27V3A (GTA→GCA) 
*157Q (TAA→CAA) 
ppdB
ppdA
conserved protein PpdB
conserved protein PpdA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/13);  total (14/13)

CTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTACAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATGCGGTATTGCGCAGGCCAAAG  >  NC_000913/2963595‑2963839
                                                                                                                  |                                                                                                                                  
cTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCtgct                                                                                                           >  2:54916/1‑137 (MQ=255)
cTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGc                                                                                                           >  2:268245/1‑140 (MQ=255)
cTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGc                                                                                                           >  1:78742/1‑140 (MQ=255)
           ttCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCACCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCg                                                                                                <  1:124133/140‑1 (MQ=255)
                                    cGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAaccacc                                                                       <  1:268245/140‑1 (MQ=255)
                                              ccAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACt                                                             <  2:78742/140‑1 (MQ=255)
                                                        tACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGa                                                   >  1:34657/1‑140 (MQ=255)
                                                                tGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTg                                                    >  2:72395/1‑131 (MQ=255)
                                                                tGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTg                                                    <  1:72395/131‑1 (MQ=255)
                                                                  aTCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAg                                         >  1:283795/1‑140 (MQ=255)
                                                                  aTCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAg                                         >  1:315208/1‑140 (MQ=255)
                                                                  aTCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTACTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAg                                         >  1:260022/1‑140 (MQ=255)
                                                                    cGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCg                                       >  1:289065/1‑140 (MQ=255)
                                                                       cATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAat                                    >  2:110387/1‑140 (MQ=255)
                                                                         tAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAatat                                  >  1:230637/1‑140 (MQ=255)
                                                                             aaTCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGccc                              <  1:110387/140‑1 (MQ=255)
                                                                               tCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGc                           <  1:54916/140‑1 (MQ=255)
                                                                                cAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTACTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGACCGGc                           <  2:260022/140‑1 (MQ=255)
                                                                                 aaCACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGcc                               >  2:31109/1‑135 (MQ=255)
                                                                                 aaCACTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGcc                               <  1:31109/135‑1 (MQ=255)
                                                                                     cTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATg                      <  2:230637/140‑1 (MQ=255)
                                                                                     cTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATg                      <  2:283795/140‑1 (MQ=255)
                                                                                     cTTCCAGCAGAGAAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATg                      <  2:315208/140‑1 (MQ=255)
                                                                                               agaAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATGCGGTATTgcg            <  2:34657/140‑1 (MQ=255)
                                                                                               agaAAAACCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATGCGGTATTgcg            <  1:197641/140‑1 (MQ=255)
                                                                                                 aaaaaCCTTGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATGCGGTATTGCGCa          >  2:253668/1‑140 (MQ=255)
                                                                                                         tGCTCTTTTGCAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATGCGGTATTGCGCAGGCCaaag  >  1:293185/1‑140 (MQ=255)
                                                                                                                  |                                                                                                                                  
CTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTACAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATGCGGTATTGCGCAGGCCAAAG  >  NC_000913/2963595‑2963839

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