Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,828 Δ1 bp 15.4% intergenic (+28/‑123) PP_0414 → / → rpe polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,8280G.15.4% 70.6 / 9.6 26intergenic (+28/‑123)PP_0414/rpepolyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (8/14);  new base . (2/2);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCC  >  NC_002947/503686‑503977
                                                                                                                                              |                                                                                                                                                      
cccGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTaa                                                                                                                                                 >  1:158195/1‑150 (MQ=255)
          gCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAaga                                                                                                                                      <  1:205202/150‑1 (MQ=255)
                gcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTAAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCgg                                                                                                                                <  2:147220/150‑1 (MQ=255)
                             tGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGTACAggcaagcca                                                                                                                   <  2:304116/150‑8 (MQ=25)
                              gATTGCGTTTTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTGGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTg                                                                                                                  <  2:51301/150‑1 (MQ=255)
                                 tGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                               <  2:653620/150‑1 (MQ=255)
                                          ggTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                      <  2:541331/150‑1 (MQ=255)
                                               tgttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGCCACAGGTTGCGTTGGTCGACACCAATTGATa                                                                                                 <  1:345728/150‑1 (MQ=255)
                                                gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                                <  1:143522/150‑1 (MQ=255)
                                                               ggAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATc                                                                                 <  1:457057/150‑1 (MQ=255)
                                                                     gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAagag                                                                                                                                     >  1:189706/1‑92 (MQ=255)
                                                                     gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAagag                                                                                                                                     <  2:189706/92‑1 (MQ=255)
                                                                                              aGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGt                                                  <  1:304239/150‑1 (MQ=255)
                                                                                                 aCTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                                <  2:736020/101‑1 (MQ=255)
                                                                                                 aCTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                                >  1:736020/1‑101 (MQ=255)
                                                                                                             gCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAgg                                   >  2:228517/1‑150 (MQ=255)
                                                                                                                   tgatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                                          >  2:384869/1‑121 (MQ=255)
                                                                                                                   tgatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                                          <  1:384869/121‑1 (MQ=255)
                                                                                                                        gCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCa                        >  1:681149/1‑150 (MQ=255)
                                                                                                                                    cGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGCTCAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTAt            >  2:655479/1‑150 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                                <  2:232096/48‑1 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                                >  1:232096/1‑48 (MQ=255)
                                                                                                                                      gCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAttt                                            >  1:502283/1‑116 (MQ=38)
                                                                                                                                      gCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAttt                                            <  2:502283/116‑1 (MQ=38)
                                                                                                                                       ccTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGc         <  2:158195/150‑1 (MQ=255)
                                                                                                                                         tctTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccc       >  1:306634/1‑150 (MQ=255)
                                                                                                                                              gcgGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTcc  >  1:433449/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
CCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCC  >  NC_002947/503686‑503977

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