Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,602,066 T→A 11.1% intergenic (+79/‑88) PP_1404 → / → PP_5469 hypothetical protein/SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,602,0660TA11.1% 114.7 / 6.0 45intergenic (+79/‑88)PP_1404/PP_5469hypothetical protein/SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base A (3/2);  total (23/22)
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.62e-01

CTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTTCTCGA  >  NC_002947/1601919‑1602209
                                                                                                                                                   |                                                                                                                                               
ctTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCtt                                                                                                                                                >  5:167465/1‑149 (MQ=255)
     ccAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCg                                                                                                                                           >  7:113711/1‑149 (MQ=255)
      cAGTACGTGCCGCCAACCCATGCACGCAAGGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACCGGCGACAGGGTGTCTGGCACCGGCTTTGCCgg                                                                                                                                          >  5:25257/1‑149 (MQ=255)
         cACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGACACCGGCAAAGCCGgtgc                                                                                                                                       <  6:204542/149‑2 (MQ=255)
             tGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                   >  3:119977/1‑149 (MQ=255)
             tGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                   <  6:52196/149‑1 (MQ=255)
                cgccAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTCGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCggg                                                                                                                                <  5:328396/149‑1 (MQ=255)
                   cAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCggg                                                                                                                                <  1:12628/146‑1 (MQ=255)
                   cAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCggg                                                                                                                                >  2:12628/1‑146 (MQ=255)
                       ccATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGccc                                                                                                                         >  3:307943/1‑149 (MQ=255)
                                 aaCGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAg                                                                                                               <  3:249555/149‑1 (MQ=255)
                                          cccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTg                                                                                                       <  6:321895/148‑1 (MQ=255)
                                           ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTggg                                                                                                     >  6:275465/1‑149 (MQ=255)
                                           ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTggg                                                                                                     >  6:270289/1‑149 (MQ=255)
                                           ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCAAAGCCGgtgc                                                                                                                                       >  7:108908/1‑114 (MQ=255)
                                           ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCAAAGCCGgtgc                                                                                                                                       <  8:108908/115‑2 (MQ=255)
                                            cAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTgggg                                                                                                    <  6:181005/149‑1 (MQ=255)
                                             aaCCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTgggg                                                                                                    >  2:227611/1‑148 (MQ=255)
                                               ccGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCAc                                                                                                 <  4:307943/149‑1 (MQ=255)
                                                     ggTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTa                                                                                             <  4:119977/147‑1 (MQ=255)
                                                         aTTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAg                                                                                       <  6:205536/149‑1 (MQ=255)
                                                             tcGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGt                                                                                   >  8:292316/1‑149 (MQ=255)
                                                                 cTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCg                                                                               >  8:159487/1‑149 (MQ=255)
                                                                    aaGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCtt                                                                            <  8:246836/149‑1 (MQ=255)
                                                                    aaGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTACCGCtt                                                                            >  8:167193/1‑149 (MQ=255)
                                                                        aGCGTCCAGCCTTGCTGGATTTCCTGTGGACACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTccc                                                                         >  1:27022/1‑148 (MQ=255)
                                                                          cGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAg                                                                      <  7:349892/149‑1 (MQ=255)
                                                                           gTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAgg                                                                     <  8:216487/149‑1 (MQ=255)
                                                                            tCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGc                                                                    >  6:14380/1‑149 (MQ=255)
                                                                               aGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCGTTGTAAACAGGGGACAGGGTGTCTGGAACCGGCTTTGCCGGTGTTCGCGGGAAAGCCCGATCCCAAAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGcccc                                                                 <  6:167465/149‑1 (MQ=38)
                                                                                 ccTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCt                                                               >  6:66341/1‑149 (MQ=255)
                                                                                   ttGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCttt                                                             >  2:177209/1‑149 (MQ=255)
                                                                                      cTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGGCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCttttt                                                           <  3:212999/148‑1 (MQ=255)
                                                                                      cTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCttttt                                                           <  7:159487/148‑1 (MQ=255)
                                                                                         gATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATg                                                       >  4:80162/1‑149 (MQ=255)
                                                                                         gATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATg                                                       >  1:26721/1‑149 (MQ=255)
                                                                                                 gtgGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAg                                               >  4:186388/1‑149 (MQ=255)
                                                                                                    ggCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTAc                                            <  8:113711/149‑1 (MQ=255)
                                                                                                        caGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCt                                        <  6:311636/149‑1 (MQ=255)
                                                                                                                      aaaCAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGcc                          >  1:23152/1‑149 (MQ=255)
                                                                                                                               aCAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATgcgcg                 <  1:227611/149‑1 (MQ=255)
                                                                                                                                         gaacACCGGCAAAGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGttt       >  4:261559/4‑149 (MQ=255)
                                                                                                                                          aacACCGGCAAAGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTGCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGACGGCTGAGCGCGCTCAGTTtc      >  8:121138/3‑149 (MQ=38)
                                                                                                                                           gCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTtc      <  3:186388/148‑1 (MQ=255)
                                                                                                                                            cACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATACGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTtct     <  6:25257/148‑1 (MQ=255)
                                                                                                                                              ccGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTTCTCGa  >  4:162343/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTTCTCGA  >  NC_002947/1601919‑1602209

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