Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,662 C→T 18.1% intergenic (+74/‑222) PP_2661 → / → PP_2662 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,048,6620CT18.1% 43.9 / 3.3 22intergenic (+74/‑222)PP_2661/PP_2662hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/11);  new base T (2/2);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGCCACCCCACGCC  >  NC_002947/3048528‑3048810
                                                                                                                                      |                                                                                                                                                    
gcagcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTTGCTTCCTGAAGCCAGGAGATTCCGCTCAGAGAATGTAATTGGGGCTGCCTTGCCGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCAGccccaa                                                                                                                                       <  2:745941/150‑3 (MQ=37)
              gTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTa                                                                                                                         <  1:155216/150‑1 (MQ=255)
                       tgaccGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAccc                                                                                                                <  1:733773/147‑1 (MQ=255)
                                                    gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcgg                                                                                                                                             <  1:651072/92‑1 (MQ=255)
                                                    gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcgg                                                                                                                                             >  2:651072/1‑92 (MQ=255)
                                                     cTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAg                                                                                                                                                 <  2:950/87‑1 (MQ=255)
                                                     cTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAg                                                                                                                                                 >  1:950/1‑87 (MQ=255)
                                                         gtgagtccAGGAGATTCCGCTCAGATAAAGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTCATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGaa                                                                              <  2:252861/144‑1 (MQ=255)
                                                         acgggtccAGGACAGTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGaa                                                                              <  2:442262/144‑1 (MQ=255)
                                                               ccAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCACCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGa                                                                        <  1:649743/150‑1 (MQ=255)
                                                                 aGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAc                                                                                                                  >  1:515181/1‑106 (MQ=255)
                                                                 aGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAc                                                                                                                  <  2:515181/106‑1 (MQ=255)
                                                                               aGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCgg                                                        <  1:629820/150‑1 (MQ=255)
                                                                                     tGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCaa                                                  >  2:334119/1‑150 (MQ=255)
                                                                                          ttGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCAc                                             <  2:440419/150‑1 (MQ=255)
                                                                                               gctgctTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGc                                        <  1:334119/150‑1 (MQ=255)
                                                                                                             cTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCAGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGcc                          >  1:337774/1‑150 (MQ=255)
                                                                                                              ttCGCGGGCTCGCCCGCGAAGGGCTGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGAc                                                                       <  1:26205/104‑1 (MQ=39)
                                                                                                              ttCGCGGGCTCGCCCGCGAAGGGCTGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGAc                                                                       >  2:26205/1‑104 (MQ=38)
                                                                                                              ttCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCGGCTCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACCCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACAAAATCACCGGGCCGTACCCTCCCGccc                         >  2:395700/1‑150 (MQ=255)
                                                                                                                       gAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACt                >  2:584491/1‑150 (MQ=255)
                                                                                                                                   ggCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGCCACCCCACg    >  2:96665/1‑150 (MQ=255)
                                                                                                                                     ccGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGCCACCCCACGcc  >  1:693384/1‑150 (MQ=255)
                                                                                                                                      |                                                                                                                                                    
GCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGCCACCCCACGCC  >  NC_002947/3048528‑3048810

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