Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,895,519 G→A 15.2% intergenic (‑75/‑72) PP_4303 ← / → PP_4304 hypothetical protein/VIC family cation transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,895,5190GA15.2% 72.6 / 6.6 33intergenic (‑75/‑72)PP_4303/PP_4304hypothetical protein/VIC family cation transporter
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base A (2/3);  total (15/18)
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.63e-01

TTCTGCAAAAAGCCGTAAGCGTTATATGACAGTGTTCGCCGAGCCAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTC  >  NC_002947/4895371‑4895668
                                                                                                                                                    |                                                                                                                                                      
ttCTGCAAAAAGCCGTAAGCGTTATATGACAGTGTTCGCCGAGCCAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGc                                                                                                                                                       <  1:4158/150‑1 (MQ=255)
     cAAAAAGCCGTAAGCGTTATATGACAGTGTTCGCCGAGCCAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTACCAGtt                                                                                                                                                  >  1:387548/1‑150 (MQ=255)
         aaGCCGTAAGCGTTATATGACAGTGTTCGCCGAGCCAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGc                                                                                                                                             <  2:432954/150‑1 (MQ=255)
                    gTTATATGACAGTGTTCGCCGAGCCAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGt                                                                                                                                  <  1:507319/150‑1 (MQ=255)
                                tgtTCGCCGAGCCAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGAc                                                                                                                      <  2:469836/150‑1 (MQ=255)
                                         aGCCAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTc                                                                                                                                                 >  2:559275/1‑115 (MQ=255)
                                         aGCCAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTc                                                                                                                                                 <  1:559275/115‑1 (MQ=255)
                                           ccAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGAc                                                                                                                      >  1:402353/1‑139 (MQ=255)
                                           ccAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGAc                                                                                                                      <  2:402353/139‑1 (MQ=255)
                                                             ttACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTACCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGa                                                                                         >  1:73896/1‑150 (MQ=255)
                                                                cccTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATca                                                                                      >  2:553901/1‑150 (MQ=255)
                                                                      tCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGAGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGcc                                                                                >  2:447408/1‑150 (MQ=255)
                                                                      tCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGAGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGcc                                                                                >  2:447397/1‑150 (MQ=255)
                                                                        aTGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGAc                                                                                                                      >  2:391079/1‑110 (MQ=255)
                                                                        aTGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGAc                                                                                                                      <  1:391079/110‑1 (MQ=255)
                                                                         tGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATg                                                                             >  2:602572/1‑150 (MQ=255)
                                                                                 tcACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAAcccc                                                                     <  1:553901/150‑1 (MQ=255)
                                                                                   acacacTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACGTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGc                                                                   >  1:596296/1‑150 (MQ=255)
                                                                                   acacacTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGc                                                                   >  1:332602/1‑150 (MQ=255)
                                                                                     acacTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGcaa                                                                 >  1:43461/1‑150 (MQ=255)
                                                                                     acacTTGCCCGCGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAACAGTCGAGCTCAAGCGCATTTTGCCGGTTCCAGGTCCGCACACCGACATTTTAGGCTTGTTGCTTATGCACCTGGATCACAAGCCATGCAAACCCCGcaa                                                                 <  1:447397/150‑1 (MQ=255)
                                                                                     acacTTCCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAACCTGACTTCAGACGTGTGCTCTTCTGTTGCAGGTCCGTACGCTGCCCTTTTAGGCTTGTTCCTGACGCAGCCGGATCACAAGCCATGCAAACCCCGcaa                                                                 <  1:447408/150‑1 (MQ=255)
                                                                                      cacTTGCCCGTGCCTTACGCTATTCTGCGCCCACCACTCAGATACCTAAACCTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGcaat                                                                <  2:320118/150‑1 (MQ=255)
                                                                                         ttGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTACCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGcaatcaa                                                             <  2:387548/150‑1 (MQ=255)
                                                                                                        gCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGATCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATCAGCCAGCGGcctt                                              >  2:438401/1‑147 (MQ=255)
                                                                                                        gCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGctct                                              <  2:261501/150‑1 (MQ=255)
                                                                                                          tATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTc                                            <  2:43461/150‑1 (MQ=255)
                                                                                                          tATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTACCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTc                                            <  2:73896/150‑1 (MQ=255)
                                                                                                            ttCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTCGt                                          >  1:300921/1‑150 (MQ=255)
                                                                                                                gcgcCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTCGTAATc                                      <  2:332602/150‑1 (MQ=255)
                                                                                                                  gcCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTCGTCATCGt                                    >  1:418030/1‑150 (MQ=255)
                                                                                                                                         cTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCg             <  1:602572/150‑1 (MQ=255)
                                                                                                                                                   tGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCtt   <  2:300921/150‑1 (MQ=255)
                                                                                                                                                    gCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTc  <  1:364313/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TTCTGCAAAAAGCCGTAAGCGTTATATGACAGTGTTCGCCGAGCCAAGCCCCGTCGTTCAGTTACCCTTTTCATGACTGTCTCACACACTTGCCCGTGCCTGACGCTATTCTGCGCCCACCACTCAGATACCTAAAACTGGCATTTTTGCCAGTTCAGCCTGTTGCAGGTCCGGACACTGACCTTTTAGGCTTGTTGCTTACGCACCCGGATCACAAGCCATGCAAACCCCGCAATCAATACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTC  >  NC_002947/4895371‑4895668

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