Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,691,783 C→A 16.0% intergenic (+50/‑121) PP_1487 → / → PP_1488 hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,691,7830CA16.0% 122.5 / 14.7 50intergenic (+50/‑121)PP_1487/PP_1488hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base A (4/4);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCC  >  NC_002947/1691641‑1691925
                                                                                                                                              |                                                                                                                                              
cGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTgcg                                                                                                                                          <  4:74967/149‑1 (MQ=255)
  tGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTgcggc                                                                                                                                        >  2:316505/1‑149 (MQ=255)
         cTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGt                                                                                                                                 >  2:236299/1‑149 (MQ=255)
         cTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTcgat                                                                                                                                 <  7:244741/149‑4 (MQ=39)
           ggACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCACGGGCGGCGCTcgatct                                                                                                                               >  7:250069/1‑144 (MQ=255)
            gACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTcc                                                                                                                                   >  4:46576/1‑144 (MQ=255)
            gACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTcc                                                                                                                                   <  3:46576/144‑1 (MQ=255)
            gACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGtct                                                                                                                               >  6:215802/1‑148 (MQ=255)
            gACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCACGGGCGGCGCTCGAtctc                                                                                                                              <  8:250069/149‑1 (MQ=255)
             aCGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTca                                                                                                                             <  3:265698/149‑1 (MQ=255)
                aGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAg                                                                                                                          <  7:29910/149‑1 (MQ=38)
                  cTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGgc                                                                                                                        <  3:141753/149‑1 (MQ=255)
                       ccGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCaa                                                                                                                    >  7:154111/1‑148 (MQ=255)
                       ccGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCaa                                                                                                                    >  1:113964/1‑148 (MQ=255)
                          ccAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAAtt                                                                                                                <  5:215802/149‑1 (MQ=255)
                           cAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAAtt                                                                                                                >  3:73506/1‑148 (MQ=255)
                           cAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTc                                                                                                               <  4:80664/149‑1 (MQ=255)
                                  ggCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTgacga                                                                                                        <  1:316505/149‑1 (MQ=255)
                                  ggCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTgacga                                                                                                        >  1:296959/1‑149 (MQ=255)
                                        aCCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATca                                                                                                  <  5:286623/149‑1 (MQ=255)
                                        aCCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATca                                                                                                  <  1:236299/149‑1 (MQ=255)
                                        aCCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTca                                                                                                                             >  5:160651/1‑122 (MQ=37)
                                        aCCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTca                                                                                                                             <  6:160651/122‑1 (MQ=37)
                                              gCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCaa                                                                                            <  2:113964/149‑1 (MQ=255)
                                                ttGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATa                                                                                          <  8:154111/149‑1 (MQ=255)
                                                   cgcCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCACTACCg                                                                                       >  8:142420/1‑149 (MQ=255)
                                                         gAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGAc                                                                                 >  6:96163/1‑149 (MQ=255)
                                                           aCTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACct                                                                               >  3:114164/1‑149 (MQ=255)
                                                             tGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACctct                                                                             >  7:140885/1‑149 (MQ=255)
                                                               gAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTgg                                                                           <  2:234965/149‑1 (MQ=255)
                                                                      aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCaa                                                                                            <  4:282213/125‑1 (MQ=255)
                                                                      aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCaa                                                                                            >  3:282213/1‑125 (MQ=255)
                                                                            ggTTGGCTGCAGTGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGc                                                              >  4:6372/1‑149 (MQ=255)
                                                                            ggTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGc                                                              <  4:73506/149‑1 (MQ=255)
                                                                                      aGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTccc                                                    <  1:190563/149‑1 (MQ=255)
                                                                                          tAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGAGCCGCCCGTGCGGCGCTCCGTCTCACAGGATCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTAGGCGATTTTTCCCGCTc                                                >  1:199565/1‑149 (MQ=255)
                                                                                              ggCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTtg                                            >  3:161545/1‑149 (MQ=255)
                                                                                                 ggTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCACGCTCGTtgtg                                          >  8:62503/1‑148 (MQ=255)
                                                                                                                 tATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTGTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAGGACCGCTGGGCTATACTTCGGCGAGTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTg                           >  1:222585/1‑147 (MQ=255)
                                                                                                                     gAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgt                      <  3:6372/148‑1 (MQ=255)
                                                                                                                        cgcCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAt                        <  4:266234/143‑1 (MQ=255)
                                                                                                                        cgcCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAt                        >  3:266234/1‑143 (MQ=255)
                                                                                                                           cTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCGTGACGACAATCACATCATTACCGTAAGAGCTCTGGGGTATACTTCGGCGATTTTTCCCGCTGGTTGTGAATTGCTGAGGGCTGATTGTgaagaag               >  2:21358/1‑148 (MQ=38)
                                                                                                                             gtgAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTg             >  6:159291/1‑149 (MQ=255)
                                                                                                                              tgAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTg             >  2:303244/1‑148 (MQ=38)
                                                                                                                                agacggAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACTACAATCACATCAATACCGTAAGACCTATGGGCTATACTTCGGCGATTTTTCCAGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGaa          >  7:150130/6‑148 (MQ=38)
                                                                                                                                agaTCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGa           <  8:140885/148‑1 (MQ=255)
                                                                                                                                agaTCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTATGAATTCCTGAGGGCTGAATGTGAAGAACTGGAc          >  1:36690/1‑149 (MQ=255)
                                                                                                                                     gAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTgc     <  4:114164/149‑1 (MQ=255)
                                                                                                                                     gAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTgc     <  7:62503/149‑1 (MQ=255)
                                                                                                                                         gccgccCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTACTGAGGGCTGAATGTGAAGAACTGTACATTgcgac  >  2:35928/1‑146 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCC  >  NC_002947/1691641‑1691925

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