Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,159,187 C→T 14.3% intergenic (+91/‑118) PP_3661 → / → PP_3662 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,159,1870CT14.3% 63.2 / 4.0 28intergenic (+91/‑118)PP_3661/PP_3662membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base T (2/2);  total (14/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 = 7.43e-01

GCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCAT  >  NC_002947/4159042‑4159319
                                                                                                                                                 |                                                                                                                                    
gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc                                                                                                                                  >  2:330951/1‑150 (MQ=255)
  tGGTTGGCGGGGGGCTGGAGATTAACATCGGGGTGCAGAAGCTGTTTTGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                                                                                <  2:387123/150‑1 (MQ=38)
                   gCGATTCCCCTGGGGGTGCAGAAGCTGTTGTGCTGATAGGAGAGCGAGTCGGCTGTTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCTCAAGACGGCAAGGTCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                                                                                               <  2:163103/150‑1 (MQ=255)
                       ttACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCa                                                                                                           <  2:508081/150‑1 (MQ=255)
                                       gAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccc                                                                                           >  1:158675/1‑150 (MQ=255)
                                               tGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGtt                                                                                   <  2:391819/150‑1 (MQ=255)
                                                 gtcTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCAGCTCATACCGGCGGGTGCTGCCCTGAAGACGTCAAGGCCACTGTCGAAACGGATTCATTCTTTAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTga                                                                                 >  1:391819/3‑149 (MQ=255)
                                                              ccGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtct                                                                    <  2:405229/150‑1 (MQ=255)
                                                                 aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTg                                                                                  <  1:76112/133‑1 (MQ=255)
                                                                 aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTg                                                                                  >  2:76112/1‑133 (MQ=255)
                                                                           ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                                                                                               >  1:502473/1‑94 (MQ=255)
                                                                           ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                                                                                               <  2:502473/94‑1 (MQ=255)
                                                                                        ccGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaa                                               <  1:54721/145‑1 (MQ=255)
                                                                                        ccGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaa                                               >  2:54721/1‑145 (MQ=255)
                                                                                              cTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg                                                            <  2:491035/126‑1 (MQ=255)
                                                                                              cTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg                                                            >  1:491035/1‑126 (MQ=255)
                                                                                                    gtgggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCTcac                              >  1:21518/1‑150 (MQ=255)
                                                                                                    gtgggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGcac                              >  1:21566/1‑150 (MQ=255)
                                                                                                      gggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAg                            >  2:376482/1‑150 (MQ=255)
                                                                                                      gggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTCCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAg                            <  1:239984/150‑1 (MQ=255)
                                                                                                            tgcCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGa                      >  1:126715/1‑150 (MQ=255)
                                                                                                            tgcCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGa                      <  1:330951/150‑1 (MQ=255)
                                                                                                             gcCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGAc                     <  1:376482/150‑1 (MQ=255)
                                                                                                                        ggCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTc          >  2:487520/1‑150 (MQ=255)
                                                                                                                         gCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCa         >  1:7860/1‑150 (MQ=255)
                                                                                                                           aaGGCCCCTGTAGGAGCGGGTTTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAAc       <  2:21566/150‑1 (MQ=255)
                                                                                                                           aaGGCCCCTGTAGGAGCGGGTTTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAAc       <  2:21518/150‑1 (MQ=255)
                                                                                                                                cccTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCAt  >  2:117163/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                    
GCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCAT  >  NC_002947/4159042‑4159319

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: