Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,951,161 A→T 100% intergenic (‑125/+70) 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,1610AT96.7% 61.7 / ‑6.2 21intergenic (‑125/+70)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (9/11);  minor base G (1/0);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑TTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGG  >  NC_002947/3951021‑3951300
                                                                                                                                              |                                                                                                                                           
aGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGa                                                                                                                                       <  8:103051/147‑1 (MQ=255)
  gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGTGACGc                                                                                                                                    >  3:102115/1‑148 (MQ=255)
                                tttcgcgacgcaaggccgctcccacagggCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGa                                                                                                                                       <  4:254688‑M1/86‑1 (MQ=255)
                                tttcgcgacgcaaggccgctcccacagggCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGa                                                                                                                                       >  3:254688‑M1/30‑115 (MQ=255)
                                                       aGCGGCCTTGCGTCGCTAACGGGGTGCGCAGCACCCCCGGCAGTGGCAGCTCAAGCGCATAAAGAGGGGAGCGCGGCGCGCTC‑‑TTTTCGTGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCa                                                                              <  4:102115/149‑1 (MQ=255)
                                                                   tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCacc                                                                   >  1:14292/1‑148 (MQ=255)
                                                                              cTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTGGCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg                                                         >  1:103559/1‑149 (MQ=255)
                                                                                  gcagcaGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCacc                                                                   >  5:283855/1‑133 (MQ=255)
                                                                                  gcagcaGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCacc                                                                   <  6:283855/133‑1 (MQ=255)
                                                                                         ccccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                             >  4:231800/1‑148 (MQ=255)
                                                                                            cGGCTGTGGCAGCTCATGCCCATAACCCGGGCCGCGCGCCGCGCGC‑‑CTTGCGCTACGAATGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGc                                         <  3:69725/149‑1 (MQ=255)
                                                                                                    gCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcg                                 >  6:212784/1‑149 (MQ=255)
                                                                                                     cAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcgc                                <  3:231800/149‑1 (MQ=255)
                                                                                                        cTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTgg                             <  2:103559/149‑1 (MQ=255)
                                                                                                                cAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTggg                            <  3:277760/142‑1 (MQ=255)
                                                                                                                cAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTggg                            >  4:277760/1‑142 (MQ=255)
                                                                                                                 aGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGAAACTCACCGACCGGCCAGCCGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgcc                    >  8:291312/1‑149 (MQ=255)
                                                                                                                    gcgaaaggAGCGCGCAGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCt                 >  6:43446/7‑149 (MQ=255)
                                                                                                                         gggAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCgg            <  5:212784/149‑1 (MQ=255)
                                                                                                                          ggAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGt           <  2:14292/149‑1 (MQ=255)
                                                                                                                                   cgcgcTC‑‑CTTGCGCTACTCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCggg  <  1:78230/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
AGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑TTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGG  >  NC_002947/3951021‑3951300

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