Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,316,485 G→A 9.3% intergenic (+51/‑98) PP_mr15 → / → PP_1150 P20 ncRNA/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,316,4850GA 9.3% 111.0 / 2.5 43intergenic (+51/‑98)PP_mr15/PP_1150P20 ncRNA/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base A (2/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

TTCATCTGGAAGGACTCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTC  >  NC_002947/1316339‑1316630
                                                                                                                                                  |                                                                                                                                                 
ttCATCTGGAAGGACTCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCCGTCCCTTTCAGTGCAGGCCAATCGa                                                                                                                                                  <  2:154317/148‑1 (MQ=255)
     cTGGAAGGACTCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAAtata                                                                                                                                             >  7:38757/1‑148 (MQ=255)
        gAAGGACTCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATAc                                                                                                                                            <  8:38757/146‑1 (MQ=255)
              cTCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATtagc                                                                                                                                   >  5:220892/1‑149 (MQ=255)
              cTCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATtagc                                                                                                                                   >  7:258904/1‑149 (MQ=255)
              cTCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCAAATATACATGATtagc                                                                                                                                   >  3:149942/1‑149 (MQ=255)
              cTCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCAAATATACATGATtagc                                                                                                                                   >  3:144580/1‑149 (MQ=255)
               tCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATtagct                                                                                                                                  >  7:59413/1‑149 (MQ=255)
                cGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATtagcta                                                                                                                                 <  7:120119/149‑1 (MQ=255)
                                      aCCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCttt                                                                                                           <  6:220892/149‑1 (MQ=255)
                                      aCCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCAAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCt                                                                                                             <  4:144580/147‑1 (MQ=255)
                                      aCCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCAAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCt                                                                                                             <  4:149942/147‑1 (MQ=255)
                                             gggTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGATAGCTAATTAATCGTTTCCGCTTTATCCTGc                                                                                                    >  2:271010/1‑149 (MQ=255)
                                              ggTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTTCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCt                                                                                                   <  8:128302/149‑1 (MQ=255)
                                                  cacGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCtcatc                                                                                               <  3:220159/149‑1 (MQ=255)
                                                   acGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCtcatca                                                                                              >  3:300875/1‑149 (MQ=255)
                                                     gAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCAtt                                                                                            >  1:35337/1‑149 (MQ=255)
                                                               ggCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACa                                                                                                                                           >  2:228329/1‑92 (MQ=255)
                                                               ggCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACa                                                                                                                                           <  1:228329/92‑1 (MQ=255)
                                                                gCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCtttt                                                                                 <  8:128957/149‑1 (MQ=255)
                                                                 cATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTCTCCTGCTCATCATTGTCTATCTTTTc                                                                                <  8:258904/149‑1 (MQ=255)
                                                                  aTTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCa                                                                               <  8:68023/149‑1 (MQ=255)
                                                                      aGGGCGACTTCGGTCGCCCA‑TTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCaa                                                                           <  8:173278/148‑1 (MQ=255)
                                                                           gACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACtt                                                                      <  8:59413/149‑1 (MQ=255)
                                                                                   tCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCg                                                               >  2:349424/1‑148 (MQ=255)
                                                                                      cccATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAg                                                           <  1:168800/149‑1 (MQ=255)
                                                                                        cATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGt                                                         >  6:344085/1‑149 (MQ=255)
                                                                                         aTTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTc                                                        >  4:197181/1‑149 (MQ=255)
                                                                                                            tttATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCt                                     >  4:186891/1‑149 (MQ=255)
                                                                                                                  gTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAg                                <  3:186891/148‑1 (MQ=255)
                                                                                                                   tCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAgg                               <  4:305197/148‑1 (MQ=255)
                                                                                                                       aGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGcc                             <  3:197181/146‑1 (MQ=255)
                                                                                                                              ccTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCgcgc                   <  5:344085/149‑1 (MQ=255)
                                                                                                                               cTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGcc                  >  2:53545/1‑149 (MQ=255)
                                                                                                                                     gTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTGCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCCGGCCCTGTTCGCGCCcagcag            >  4:247250/1‑149 (MQ=255)
                                                                                                                                         aGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTc        >  2:283273/1‑149 (MQ=255)
                                                                                                                                           gCCAATCGAATATACATGATTAGCTAGGTAATTAATCGTTTCCGCGTTATCCTGCTCATCATTGTCTATCTTTTAAGCAAGACTTTTCACCGCAAGGTCCGCCTCGTGGCCATCACCTTTCAGGCCCTGTTCTCGCACAGCAGCCTCGc      >  2:284807/1‑149 (MQ=255)
                                                                                                                                           gCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGc      >  3:239200/1‑149 (MQ=255)
                                                                                                                                             aaaTCGAATATAGATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGccc    <  4:300875/148‑1 (MQ=255)
                                                                                                                                              aaTCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCg       <  3:87518/145‑1 (MQ=255)
                                                                                                                                              aaTCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCg       >  4:87518/1‑145 (MQ=255)
                                                                                                                                                tCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTc  >  7:78451/1‑148 (MQ=255)
                                                                                                                                                tCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTCTCTATGTTTTCAGCAAGACTTTTCCACCCAAGCTCCGCCACGTGCCCATAACCTTTCAGGACCTTTTCCCGCCCAGCAGCCTCGCCATc  >  3:2360/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTCATCTGGAAGGACTCGAACAAGCGCAAGCCGCGCTGACCTTCCGGGTACACGAGCACGCAAGGCATTTAGGGCGACTTCGGTCGCCCATTTTTTTTGCCTGCATTTTTTATTGTCTGAGCTGTCCCTTTCAGTGCAGGCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCATTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTC  >  NC_002947/1316339‑1316630

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: ‑