Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,884,398 G→T 9.1% intergenic (‑731/‑220) 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,3980GT 9.1% 161.2 / 6.5 55intergenic (‑731/‑220)PP_1691/PP_1692hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base T (2/3);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

GCCTCGCCAATGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATCC  >  NC_002947/1884251‑1884542
                                                                                                                                                   |                                                                                                                                                
gcCTCGCCAATGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAgg                                                                                                                                                  >  7:128136/1‑148 (MQ=255)
 ccTCGCCAATGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGa                                                                                                                                                 >  3:126938/1‑148 (MQ=255)
     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                                                                                                                                        >  1:197604/1‑122 (MQ=255)
                                    gAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGTTCCTGTGGAt                                                                                                                                        <  2:197604/122‑1 (MQ=255)
                                     aaaaaaaaCACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGc                                                                                                            >  2:43856/1‑149 (MQ=255)
                                       aaaaaaCACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGt                                                                                                          <  2:44479/149‑1 (MQ=255)
                                        aaaaaCACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTa                                                                                                                    >  1:210502/1‑138 (MQ=255)
                                        aaaaaCACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTa                                                                                                                    <  2:210502/138‑1 (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                                 <  3:142024/129‑1 (MQ=255)
                                                                                                                                    ccGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTc                                 >  4:142024/1‑129 (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                <  7:137221/135‑1 (MQ=255)
                                                                                                                                               acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGaa                >  8:132648/1‑135 (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)
                                                                                                                                                   |                                                                                                                                                
GCCTCGCCAATGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATCC  >  NC_002947/1884251‑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: ‑