Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,787,907 Δ1 bp 100% coding (728/1392 nt) PP_3350 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,787,9070T.100.0% 82.6 / NA 19coding (728/1392 nt)PP_3350hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (10/9);  total (10/9)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCGGTTTTCAGCGGGCGCACAGCGAATGCATCCAGCTTGCGTCCGTCCTTGAGGTTGAG  >  NC_002947/3787789‑3788049
                                                                                                                        |                                                                                                                                              
gtccagtGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCg                                                                                                                   >  5:123426/1‑147 (MQ=255)
         cTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCg                                                                                                          >  2:101281/1‑147 (MQ=255)
                cggcCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTg                                                                                                 <  5:7987/149‑1 (MQ=255)
                                 gCCGAAAGTGTAGCGCTGCTCGTCACCCGGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTtgctgc                                                                                 >  1:20865/1‑149 (MQ=255)
                                    gAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCAt                                                                              >  1:94207/1‑149 (MQ=255)
                                                      gTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCAcc                                                           >  4:251692/1‑149 (MQ=255)
                                                      gTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCAcc                                                           >  6:278043/1‑149 (MQ=255)
                                                       tCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCg                                                          <  6:123426/149‑1 (MQ=255)
                                                         aCCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCGGt                                                        <  2:20865/149‑1 (MQ=255)
                                                            cGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCGGttt                                                      >  6:335949/1‑148 (MQ=255)
                                                              gACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCGGTTTTc                                                    <  5:278043/149‑1 (MQ=255)
                                                                   gCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGAtgttgtt                                                                                         >  7:240503/1‑107 (MQ=255)
                                                                   gCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGAtgttgtt                                                                                         <  8:240503/107‑1 (MQ=255)
                                                                     cGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCGGTTTTCAGCGGgc                                             <  2:94207/149‑1 (MQ=255)
                                                                                              ccAGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCGGTTTTCAGCGGGCGCACAGCGAATGCATCCAGCTTg                      <  5:335949/148‑1 (MQ=255)
                                                                                                aGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGAtgttg                                                                                           <  3:114940/77‑1 (MQ=255)
                                                                                                aGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGAtgttg                                                                                           >  4:114940/1‑77 (MQ=255)
                                                                                                aGGCCAAAGGCTAGCAGGTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCGGTTTTCAGCGGGCGCACAGCGAATGCATCCAGCTTGcgt                   <  1:101281/149‑1 (MQ=255)
                                                                                                                 gTCCATGTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCGGTTTTCAGCGGGCGCACAGCGAATGCATCCAGCTTGCGTCCGTCCttgaggttgag  >  6:163144/1‑149 (MQ=255)
                                                                                                                        |                                                                                                                                              
GTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTTCCACGACGCCGACAGCGGCAAGGTGCCGTAGAGGCCGTAGAATTTGATGTTGTTGTTGCTGCCATCGTCCCAACTGTCTTCACCGGTTTTCAGCGGGCGCACAGCGAATGCATCCAGCTTGCGTCCGTCCTTGAGGTTGAG  >  NC_002947/3787789‑3788049

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