Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,884,400 A→C 10.9% intergenic (‑733/‑218) PP_1691 ← / → PP_1692 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,884,4000AC10.9% 136.8 / 8.0 54intergenic (‑733/‑218)PP_1691/PP_1692hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/24);  new base C (2/4);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GCCAATGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATCC  >  NC_002947/1884256‑1884542
                                                                                                                                                |                                                                                                                                              
gCCAATGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACt                                                                                                                                              <  7:181354/147‑1 (MQ=255)
   aaTGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGTTCCTGTGGa                                                                                                                                         <  3:161321/149‑1 (MQ=255)
     tggcTTGACTGCGTGGTTAAATCAATGAAAAAAAACACTATTTTTTGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGAtt                                                                                                                                       <  5:66380/146‑1 (MQ=255)
        cTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTc                                                                                                                                    <  4:126938/149‑1 (MQ=255)
                cGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAgcgc                                                                                                                            >  6:203759/1‑149 (MQ=255)
                        aGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGctct                                                                                                                     <  7:206753/148‑1 (MQ=255)
                            aaTGAAAAAACACTCTATTTTTAGCTGGTGAAAAAATCTTCAGTTTGACCAAAGCCCCCATTTCAGGTGGGTTTGAGGGGTATGAGAGTGGGTTGTCAACCGATTTATCAACAGGAACTGTGGATTGTCCCGATCGCTTGCTCTAGCa                                                                                                                 <  1:175713/148‑1 (MQ=255)
                             aTGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACg                                                                                                               >  3:269756/1‑149 (MQ=255)
                               gAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGTTCCTGTGGAt                                                                                                                                        <  2:197604/122‑1 (MQ=255)
                               gAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGTTCCTGTGGAt                                                                                                                                        >  1:197604/1‑122 (MQ=255)
                                aaaaaaaaCACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGc                                                                                                            >  2:43856/1‑149 (MQ=255)
                                  aaaaaaCACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGt                                                                                                          <  2:44479/149‑1 (MQ=255)
                                   aaaaaCACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTa                                                                                                                    <  2:210502/138‑1 (MQ=255)
                                   aaaaaCACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTa                                                                                                                    >  1:210502/1‑138 (MQ=255)
                                       acacTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAg                                                                                                     <  1:147712/149‑1 (MQ=255)
                                          cTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCaa                                                                                                  >  3:268455/1‑149 (MQ=255)
                                                 tAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGAt                                                                                            <  5:203759/148‑1 (MQ=255)
                                                         gAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGt                                                                                     >  4:31686/1‑147 (MQ=255)
                                                         gAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTcc                                                                                   >  7:259174/1‑149 (MQ=255)
                                                               aTCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGATCTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAAt                                                                             >  7:65278/1‑149 (MQ=255)
                                                                            ccAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCaa                                                                <  1:43856/149‑1 (MQ=255)
                                                                              aaaGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGAt                                                                                            <  6:195605/119‑1 (MQ=255)
                                                                              aaaGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGAt                                                                                            >  5:195605/1‑119 (MQ=255)
                                                                                    cccATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCgg                                                        >  3:15480/1‑149 (MQ=255)
                                                                                     cgagctCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGt                                                       <  1:58164/144‑1 (MQ=255)
                                                                                           cACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGgtc                                                  <  4:268455/148‑1 (MQ=255)
                                                                                            aCGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGt                                                       >  2:58164/1‑142 (MQ=255)
                                                                                                gggTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCAc                                            <  4:105207/149‑1 (MQ=255)
                                                                                                    ttGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTg                                                     >  3:111448/1‑136 (MQ=255)
                                                                                                         gggTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACg                                   >  3:143166/1‑149 (MQ=255)
                                                                                                              tGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTc                                 >  4:247572/1‑146 (MQ=255)
                                                                                                               gagaGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGAGCCAGCAATGGGATGGACTGTCAGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAg                              >  6:177327/1‑148 (MQ=255)
                                                                                                                       gTTGTCCACCGAGTTATCCACAGGACCTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCa                     <  4:269756/149‑1 (MQ=255)
                                                                                                                          gTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCAtt                   <  3:247572/148‑1 (MQ=255)
                                                                                                                            ccACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCg                                                         >  7:134382/1‑108 (MQ=255)
                                                                                                                            ccACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCg                                                         <  8:134382/108‑1 (MQ=255)
                                                                                                                             cACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAAt               >  1:73047/1‑149 (MQ=255)
                                                                                                                               ccGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTc                                 >  4:142024/1‑129 (MQ=255)
                                                                                                                               ccGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTc                                 <  3:142024/129‑1 (MQ=255)
                                                                                                                                    ttATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAAtttttttt        <  3:31686/149‑1 (MQ=255)
                                                                                                                                     aaTCCACAGTTCCTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGaaaa                          <  6:18071/129‑1 (MQ=255)
                                                                                                                                     aaTCCACAGTTCCTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGaaaa                          >  5:18071/2‑130 (MQ=255)
                                                                                                                                      aTCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCAc                                            <  8:34152/111‑1 (MQ=255)
                                                                                                                                      aTCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCAc                                            >  7:34152/1‑111 (MQ=255)
                                                                                                                                      aTCCACACGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATTGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCa      <  8:128136/149‑1 (MQ=255)
                                                                                                                                       tCCACAGGAACTGTTGATTGTCCCGAGCGCTTGGTCTAGCTCGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCa      <  8:259174/148‑1 (MQ=255)
                                                                                                                                       tCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCaa     <  3:232666/149‑1 (MQ=255)
                                                                                                                                          acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGaa                >  8:132648/1‑135 (MQ=255)
                                                                                                                                          acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGaa                <  7:137221/135‑1 (MQ=255)
                                                                                                                                          acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGaa                >  8:137221/1‑135 (MQ=255)
                                                                                                                                          acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGaa                <  7:132648/135‑1 (MQ=255)
                                                                                                                                          acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATcc  >  7:309062/1‑149 (MQ=255)
                                                                                                                                          acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATcc  >  5:10136/1‑149 (MQ=255)
                                                                                                                                          acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATcc  >  4:250569/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GCCAATGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATCC  >  NC_002947/1884256‑1884542

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