Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,104,110 T→A 23.6% intergenic (+32/‑82) PP_5353 → / → PP_5354 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,104,1100TA23.6% 60.2 / 19.1 34intergenic (+32/‑82)PP_5353/PP_5354hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (11/15);  new base A (4/4);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

GAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCT  >  NC_002947/6103967‑6104253
                                                                                                                                               |                                                                                                                                               
gaagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                                                                                                           >  1:101792/1‑150 (MQ=255)
gaagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                                                                                                           >  1:594623/1‑150 (MQ=255)
 aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCggg                                                                                                                                          <  2:427420/150‑1 (MQ=255)
  agaaTTCGCGCTTCGACAGAGGTGTAGTATTTAATAAAATCTTGCATATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCgggg                                                                                                                                         <  2:223697/150‑1 (MQ=255)
           gcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctt                                                                                                                                >  1:222310/1‑150 (MQ=255)
             ttCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctttt                                                                                                                              <  2:75154/150‑1 (MQ=255)
                  ccAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGTCGAGAAGCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCaa                                                                                                                         >  2:87516/1‑150 (MQ=255)
                     gggCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGcc                                                                                                                           <  1:132673/145‑1 (MQ=255)
                     gggCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGcc                                                                                                                           >  2:132673/1‑145 (MQ=255)
                         gAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGc                                                                                                                  <  1:551522/150‑1 (MQ=255)
                                           tGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgcgc                                                                                                >  1:30852/1‑150 (MQ=255)
                                            ggCACATGGTTCGCGCCCAAGGCGAGACACAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCa                                                                                               <  2:101792/150‑1 (MQ=255)
                                                caTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGTTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGt                                                                                           >  1:497628/1‑150 (MQ=255)
                                                 aTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt                                                                                          <  2:594623/150‑1 (MQ=255)
                                                                ggCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCaga                                                                           <  1:89407/150‑1 (MQ=255)
                                                                         cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAg                                                                                 <  1:341368/135‑1 (MQ=255)
                                                                         cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAg                                                                                 >  2:341368/1‑135 (MQ=255)
                                                                          aGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAAc                                                                 <  1:609730/150‑1 (MQ=255)
                                                                              ttGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTg                                                             >  1:448672/1‑150 (MQ=255)
                                                                              ttGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGCGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTg                                                             <  2:135493/150‑1 (MQ=255)
                                                                                         ccGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGAt                                                  <  2:448672/150‑1 (MQ=255)
                                                                                                   tCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTc                                        >  1:250632/1‑150 (MQ=255)
                                                                                                        aCAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGa                                   <  1:74313/150‑1 (MQ=255)
                                                                                                           gCTAATTGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGc                                <  2:497628/150‑1 (MQ=255)
                                                                                                                             aaCACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgc                                                                                                  <  2:22485/66‑1 (MQ=255)
                                                                                                                             aaCACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgc                                                                                                  >  1:22485/1‑66 (MQ=255)
                                                                                                                                 aaaaaCCCCGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGccc          >  2:419516/4‑150 (MQ=255)
                                                                                                                                      ccccGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGa     <  2:250632/150‑1 (MQ=255)
                                                                                                                                      ccccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGaa                             >  2:221352/1‑126 (MQ=255)
                                                                                                                                      ccccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGaa                             <  1:221352/126‑1 (MQ=255)
                                                                                                                                       cccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGaa                                              >  2:1464/1‑108 (MQ=255)
                                                                                                                                       cccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGaa                                              <  1:1464/108‑1 (MQ=255)
                                                                                                                                       cccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAg    >  2:570049/1‑150 (MQ=255)
                                                                                                                                         cGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCt  <  1:87516/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCT  >  NC_002947/6103967‑6104253

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