Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,647 G→T 23.1% intergenic (+59/‑237) 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,6470GT23.1% 37.7 / 5.1 21intergenic (+59/‑237)PP_2661/PP_2662hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (5/11);  new base T (2/3);  total (7/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACT  >  NC_002947/3048509‑3048796
                                                                                                                                          |                                                                                                                                                     
tgCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCGGCCCTTCGCGGGCTCGCCCGCGAAg                                                                                                                                            <  1:109063/150‑1 (MQ=255)
                   gcagcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTTGCTTCCTGAAGCCAGGAGATTCCGCTCAGAGAATGTAATTGGGGCTGCCTTGCCGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCAGccccaa                                                                                                                         <  2:745941/150‑3 (MQ=37)
                                 gTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTa                                                                                                           <  1:155216/150‑1 (MQ=255)
                                          tgaccGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAccc                                                                                                  <  1:733773/147‑1 (MQ=255)
                                                                       gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcgg                                                                                                                               >  2:651072/1‑92 (MQ=255)
                                                                       gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcgg                                                                                                                               <  1:651072/92‑1 (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                                                         >  2:26205/1‑104 (MQ=38)
                                                                                                                                 ttCGCGGGCTCGCCCGCGAAGGGCTGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGAc                                                         <  1:26205/104‑1 (MQ=39)
                                                                                                                                 ttCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCGGCTCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACCCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACAAAATCACCGGGCCGTACCCTCCCGccc           >  2:395700/1‑150 (MQ=255)
                                                                                                                                          gAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACt  >  2:584491/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
TGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACT  >  NC_002947/3048509‑3048796

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