Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,854,734 C→G 7.9% intergenic (+56/‑255) PP_1658 → / → PP_1659 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,854,7340CG 7.9% 150.0 / 2.4 49intergenic (+56/‑255)PP_1658/PP_1659hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/17);  new base G (2/2);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTC  >  NC_002947/1854591‑1854877
                                                                                                                                                |                                                                                                                                               
tATTCTGCCGCTCAACCCCCATGTCGAGGAAACGGACGCACCCGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCgg                                                                                                                                            <  4:54871/149‑1 (MQ=255)
       ctgctCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCAGAGCGGCCCtg                                                                                                                                       <  6:14337/147‑2 (MQ=255)
       ctgctCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCAGAGCGGCCCt                                                                                                                                        <  3:290863/146‑1 (MQ=255)
         gctCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTcgc                                                                                                                                   >  7:92064/1‑149 (MQ=255)
              aCAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGacac                                                                                                                              >  5:27424/1‑149 (MQ=255)
                   cATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACAcgcgcg                                                                                                                         >  2:53494/1‑149 (MQ=255)
                         gAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTtctc                                                                                                                    <  8:92064/148‑1 (MQ=255)
                         gAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTtctc                                                                                                                    >  1:86208/1‑148 (MQ=255)
                          ccgAAACGGACGAAGACGCCCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTa                                                                                                                  <  8:235247/147‑1 (MQ=255)
                            gAAACGGACGAAGACCACCCGGACGCGGCTTACAACATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTaga                                                                                                                <  7:71524/149‑1 (MQ=255)
                               aCGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTagaaga                                                                                                             <  7:36311/149‑1 (MQ=255)
                                    cGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTg                                                                                                        <  3:238133/149‑1 (MQ=255)
                                    cGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTg                                                                                                        >  3:216026/1‑149 (MQ=255)
                                     gAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCGCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTgg                                                                                                       <  2:188503/149‑1 (MQ=255)
                                         acgacCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCtt                                                                                                   >  2:68067/1‑149 (MQ=255)
                                          cgacCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTg                                                                                                  >  1:271646/1‑149 (MQ=255)
                                              ccGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTc                                                                                              <  6:27424/149‑1 (MQ=255)
                                              ccGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTc                                                                                              >  2:255685/1‑149 (MQ=255)
                                                           aCAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTtc                                                                                  >  1:230335/1‑149 (MQ=255)
                                                                 tCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTag                                                                                                                 >  4:13809/1‑111 (MQ=255)
                                                                 tCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTag                                                                                                                 <  3:13809/111‑1 (MQ=255)
                                                                    cTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCAcc                                                                        >  7:58985/1‑149 (MQ=255)
                                                                      gatgatCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCg                                                                       >  1:237013/1‑148 (MQ=255)
                                                                                  tgGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGtt                                                          >  7:224861/1‑149 (MQ=255)
                                                                                   ggTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGttt                                                         >  2:216649/1‑149 (MQ=255)
                                                                                      aaGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCAc                                                                         <  1:137053/130‑1 (MQ=255)
                                                                                      aaGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCAc                                                                         >  2:137053/1‑130 (MQ=255)
                                                                                          aCGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGtt                                                  <  1:255685/149‑1 (MQ=255)
                                                                                                   ccGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCg                                         >  2:105213/1‑149 (MQ=255)
                                                                                                    cGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGt                                        >  1:268978/1‑149 (MQ=255)
                                                                                                     gCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTc                                       <  2:237013/149‑1 (MQ=255)
                                                                                                       cTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCtt                                     >  1:306235/1‑149 (MQ=255)
                                                                                                              gTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGCCACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTAGCGAAAGGGTTCGTCGCTCCCCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTcc                               >  4:200658/1‑149 (MQ=255)
                                                                                                               ttAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTcc                               >  6:244531/1‑148 (MQ=255)
                                                                                                                 aaTGGGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCG‑AAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTcc                           >  5:14780/1‑149 (MQ=255)
                                                                                                                      ggTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTtc                       <  8:224861/149‑1 (MQ=255)
                                                                                                                          cattcaAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCTGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTtcgtc                    <  1:68067/148‑1 (MQ=255)
                                                                                                                          cattcaAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTtcgtc                    <  2:230335/148‑1 (MQ=255)
                                                                                                                            ttcaAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGttt                                                         <  2:18217/109‑1 (MQ=255)
                                                                                                                            ttcaAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGttt                                                         >  1:18217/1‑109 (MQ=255)
                                                                                                                              caAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCTCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACa               >  8:130909/1‑149 (MQ=255)
                                                                                                                               aaGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACAt              >  8:16116/1‑149 (MQ=255)
                                                                                                                                  gcgaaaGGGCCGCAGAGCGGCCCTTTCGCGACACCCGCGATCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTc           >  8:75177/6‑149 (MQ=255)
                                                                                                                                   tCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCG‑AAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCtt         <  6:14780/149‑1 (MQ=255)
                                                                                                                                     aaaGGGCCGCAGAGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTGCCTCCGTTCGTCGGACATTTCt          >  8:194466/3‑147 (MQ=255)
                                                                                                                                       aGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAgg      >  2:24137/1‑149 (MQ=255)
                                                                                                                                       aGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAgg      >  2:213052/1‑149 (MQ=255)
                                                                                                                                          gCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGtt   >  2:21481/1‑149 (MQ=255)
                                                                                                                                           ccGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTc  >  2:204244/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGC‑GTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTC  >  NC_002947/1854591‑1854877

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: ‑