Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,711,029 T→G 12.5% intergenic (+289/+25) PP_2373 → / ← PP_2374 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,711,0290TG12.5% 70.8 / 5.7 32intergenic (+289/+25)PP_2373/PP_2374hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/18);  new base G (2/2);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

CCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGAAGTCGCGTGC  >  NC_002947/2710882‑2711173
                                                                                                                                                   |                                                                                                                                                
cccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCCCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGCTGGACACAGGCCGCTGTAGGCGCGGGTTCa                                                                                                                                                <  1:309163/150‑1 (MQ=255)
             cgggccTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGc                                                                                                                                   <  1:328244/150‑1 (MQ=255)
               ggccTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGATGTAGGCGCGGGTTCACCCGCGAAGAGGCCg                                                                                                                                 <  2:308160/150‑1 (MQ=255)
                     ctgGCTAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAg                                                                                                                           <  2:601474/150‑1 (MQ=255)
                     ctgGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAg                                                                                                                           <  1:194287/150‑1 (MQ=255)
                          cAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAg                                                                                                                      >  1:300875/1‑150 (MQ=255)
                                 cACCTGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCa                                                                                                               <  1:190749/150‑1 (MQ=255)
                                 cACCTGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCa                                                                                                               <  1:495785/150‑1 (MQ=255)
                                          ttCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAACGGGTCAGGTGCAggg                                                                                                      >  1:155883/1‑150 (MQ=255)
                                             aCCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCa                                                                                                   >  2:75555/1‑150 (MQ=255)
                                                         gcTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGcc                                                                                       >  1:78734/1‑150 (MQ=255)
                                                                              cccGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCg                                                                                                                    >  1:168979/1‑100 (MQ=255)
                                                                              cccGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCg                                                                                                                    <  2:168979/100‑1 (MQ=255)
                                                                                     cTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGt                                                           <  1:55797/150‑1 (MQ=255)
                                                                                                     gTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAAc                                                       <  1:112856/138‑1 (MQ=255)
                                                                                                     gTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAAc                                                       >  2:112856/1‑138 (MQ=255)
                                                                                                         acacAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAgg                                                             >  2:103887/1‑128 (MQ=255)
                                                                                                         acacAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAgg                                                             <  1:103887/128‑1 (MQ=255)
                                                                                                                    gggTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAgg                                                             >  2:162347/1‑117 (MQ=255)
                                                                                                                    gggTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAgg                                                             <  1:162347/117‑1 (MQ=255)
                                                                                                                         ggACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCtt                       <  1:90382/150‑1 (MQ=255)
                                                                                                                                ggCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCGATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTgg                <  1:75555/150‑1 (MQ=255)
                                                                                                                                 gCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTggg               >  1:451515/1‑150 (MQ=255)
                                                                                                                                     cggtgggCGCGGTGTCACCCGCTAAGAGGCCGGCACAGGTTAGCGGGTAAGGTGCAGGGGCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGaa           <  2:531046/145‑1 (MQ=255)
                                                                                                                                        tcttcgcgGGTGAACCCGCGCAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGaa                                                        <  1:208937/98‑1 (MQ=255)
                                                                                                                                        tcttcgcgGGTGAACCCGCGCAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGaa                                                        >  2:208937/5‑102 (MQ=255)
                                                                                                                                          ttcgcgGGTGAACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAgg                                                             <  2:124231/93‑1 (MQ=38)
                                                                                                                                          ttcgcgGGTGAACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAgg                                                             >  1:124231/3‑95 (MQ=38)
                                                                                                                                          ggCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGAAGTcgc      <  1:184471/150‑1 (MQ=255)
                                                                                                                                          ggCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGAAGTcgc      >  2:602187/1‑150 (MQ=255)
                                                                                                                                           gcgcGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGAAGTcgcg     <  2:20258/150‑1 (MQ=255)
                                                                                                                                              ggggTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGAAGTCGCGTgc  <  2:155883/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGAAGTCGCGTGC  >  NC_002947/2710882‑2711173

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: