Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,951,159 T→C 100% intergenic (‑123/+72) PP_3486 ← / ← PP_3487 cytochrome c/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,951,1590TC100.0% 66.1 / NA 22intergenic (‑123/+72)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/7);  total (15/7)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGG  >  NC_002947/3951024‑3951300
                                                                                                                                       |                                                                                                                                                
tCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTTCGTCGCGAAACGGCTTCGCAGCAGCCCCCGCAGTGGCCGCGCAAGCGCAGAAACAGGGGGGCCCTGCGCGCTCCGTTCGCGACGc                                                                                                                                       >  7:98639/1‑147 (MQ=255)
                                             cTGTGGGAGCGGCCTTTAGTCGCGAAAGGGCTTCGCAGCAGCCCCGGCAGTGGGAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGGTCCTTTAGCGACACCAGGCCCCTCCCACAAAAACGGAAGCCACAGAG‑‑‑CGCAAATACCCaag                                                                                     >  2:281028/1‑147 (MQ=255)
                                                         ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCACGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATa                                                                         >  3:183013/1‑149 (MQ=255)
                                                             gCGTCGCTAATGGGCTGCGCAGCCGCCCCGGCAGTGGCAGCTCAAGCGCAGAACGAGGGTCGCGCTGCGCGCTCCTTTCGCTACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACt                                                                       <  7:262231/147‑1 (MQ=255)
                                                                tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCaccg                                                                  >  3:319908/1‑149 (MQ=255)
                                                                  gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCaccgac                                                                >  7:76154/1‑149 (MQ=255)
                                                                      acgggCTGCGCACCACCCCCGGCTGTGGCAGCTCAAGCCCACAACGCGGGCCGCGCTGCGCGCTCCTTTCGCGACGAAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCaccgaccg                                                              <  1:55814/145‑1 (MQ=255)
                                                                          gCTGCGCAGCAGCCCGGGCTGTGGCAGCGCAAGCGCATAAAGCGGGTCGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc                                                        <  5:97770/149‑1 (MQ=38)
                                                                                  gcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                                >  8:111308/1‑149 (MQ=255)
                                                                                      ccccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                             >  7:70100/1‑148 (MQ=255)
                                                                                                  cAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCAGGGGGGGAGGGATTTCAGGAAGCGGCCCATGATACTAACCGACCCGCCAGCGGCAGGGTCAAAAGCTTCATg                                   >  1:307951/1‑149 (MQ=19)
                                                                                                      tCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTggg                            >  4:210913/1‑149 (MQ=255)
                                                                                                          gcgcAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAAt                        >  4:2402/1‑149 (MQ=255)
                                                                                                              aGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGa                                                                           >  2:260402/1‑94 (MQ=255)
                                                                                                              aGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGa                                                                           <  1:260402/94‑1 (MQ=255)
                                                                                                              aGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgcc                    <  3:210913/149‑1 (MQ=255)
                                                                                                              aGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgcc                    >  2:69184/1‑149 (MQ=255)
                                                                                                                  cgaaaggAGCGCGCAGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTg                >  2:61522/6‑149 (MQ=255)
                                                                                                                   gagGGTAGCGCTGCGCGCTCCTTTCGCTACGAAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTgg               <  1:281028/149‑1 (MQ=255)
                                                                                                                       ggAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGt           <  3:2402/149‑1 (MQ=255)
                                                                                                                        gAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTg          >  4:289278/1‑149 (MQ=255)
                                                                                                                                cgcgcTCCTTTCGCGACGCAAGGCCGCTCCGACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCggg  >  5:119626/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                
TCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAG‑‑‑GCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGG  >  NC_002947/3951024‑3951300

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