Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,856,165 C→T 14.8% D12D (GAC→GAT)  PP_5133 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,856,1650CT14.8% 59.5 / 4.1 27D12D (GAC→GAT) PP_5133membrane protein
Reads supporting (aligned to +/- strand):  ref base C (12/11);  new base T (2/2);  total (14/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 = 6.29e-01

GCCGATCACGCGGTTCTCGGCATAAGCCGCGATCAGGCTGAGGGGGAGTGATGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTGGCCACGCGCATCGAAACACGTGCGCTGTTCCTGGCCGAGCGCGAT  >  NC_002947/5856016‑5856312
                                                                                                                                                     |                                                                                                                                                   
gcggATCACGCGGTTCTCGGCATAAGCCGCGATCAGGCTGAGGGGGAGTGATGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGAc                                                                                                                                                     >  1:563207/4‑150 (MQ=255)
                tcGGCATAAGCCGCGATCAGGCTGAGGGGGAGTGATGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCg                                                                                                                                     >  1:412538/1‑150 (MQ=255)
                                tCAGGCTGAGGGGGAGTGATGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTgcg                                                                                                                     <  2:412538/150‑1 (MQ=255)
                                    gCTGAGGGGGAGTGATGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAGTGGAAGACCTTGTGACTGCACCCACGCCACTGGATAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCgaa                                                                                                                 >  2:546987/1‑150 (MQ=255)
                                          ggggAGTGATGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACAt                                                                                                           <  1:450471/150‑1 (MQ=255)
                                          gaggACTGCTGTATTCATGCCTGCGACCATAACACCTCACTAGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACAt                                                                                                           <  2:204697/148‑1 (MQ=255)
                                              aGTGATGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCg                                                                                                       >  2:86527/1‑150 (MQ=255)
                                                tGATGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGc                                                                                                     >  1:362311/1‑150 (MQ=255)
                                                  aTGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGATAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGccc                                                                                                   <  1:546987/150‑1 (MQ=255)
                                                  aTGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGc                                                                                                         >  2:539689/1‑144 (MQ=255)
                                                  aTGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGc                                                                                                         <  1:539689/144‑1 (MQ=255)
                                                           tGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGa                                                                                          >  2:274780/1‑150 (MQ=255)
                                                              cAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCa                                                                                       <  1:65840/150‑1 (MQ=255)
                                                                gCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAgg                                                                                     >  1:504248/1‑150 (MQ=255)
                                                                          cACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCg                                                                           >  2:439119/1‑150 (MQ=255)
                                                                                gCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAAGGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCAcgc                                                                     <  1:67888/150‑1 (MQ=255)
                                                                                        gTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGc                                                             >  1:105975/1‑150 (MQ=255)
                                                                                                     cTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCt                                                >  2:352642/1‑150 (MQ=255)
                                                                                                      ttGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTg                                               <  1:80430/150‑1 (MQ=255)
                                                                                                       tGACTTGTGCAATGGAAGACCTTTTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTgg                                              <  2:391644/150‑1 (MQ=255)
                                                                                                       tGACCTGTCCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTAACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCACTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTgg                                              <  1:165468/150‑1 (MQ=255)
                                                                                                                         aCCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTGGCCACGCGCATCGAAacac                            <  2:508440/150‑1 (MQ=255)
                                                                                                                           cTTGTGACTGCACCCACGCCACTGGATAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGAc                                                                                         >  2:30/1‑87 (MQ=255)
                                                                                                                           cTTGTGACTGCACCCACGCCACTGGATAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGAc                                                                                         <  1:30/87‑1 (MQ=255)
                                                                                                                           cTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTGGCCACGCGCATCGAAACACGt                          >  2:377725/1‑150 (MQ=255)
                                                                                                                                gACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTGGCCACGCGCATCGAAACACGTGCGCt                     >  1:222217/1‑150 (MQ=255)
                                                                                                                                   tGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTGGCCACGCGCATCGAAACACGTGCGCTGtt                  >  2:343686/1‑150 (MQ=255)
                                                                                                                                                   gACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTGGCCACGCGCATCGAAACACGTGCGCTGTTCCTGGCCGAGCGCGat  <  1:86527/150‑1 (MQ=255)
                                                                                                                                                     cAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAg                                                          <  2:496228/92‑1 (MQ=255)
                                                                                                                                                     cAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAg                                                          >  1:496228/1‑92 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCGATCACGCGGTTCTCGGCATAAGCCGCGATCAGGCTGAGGGGGAGTGATGTATTCATGCCAGCGAGGATAGCACCGGGCGTGTGGGTTATGCTTCTGCCTTGACCTGTGCAATGGAAGACCTTGTGACTGCACCCACGCCACTGGACAAGCGCTGGCTCACCGAAGCGGTACGCCTGCGCGAAGAACATGCCGGCCCGCTGGAAGACCAGGAAGCCAACCGCCACGCCCGCCAGCTAGGTGGCGACCTGGCCACGCGCATCGAAACACGTGCGCTGTTCCTGGCCGAGCGCGAT  >  NC_002947/5856016‑5856312

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