Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,884,399 A→T 9.2% intergenic (‑732/‑219) 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,3990AT 9.2% 153.1 / 5.2 54intergenic (‑732/‑219)PP_1691/PP_1692hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base T (2/3);  total (26/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 = 8.22e-01

CCTCGCCAATGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATCC  >  NC_002947/1884252‑1884542
                                                                                                                                                   |                                                                                                                                               
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                                                                                                                                        <  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                                                         <  8:134382/108‑1 (MQ=255)
                                                                                                                                ccACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCg                                                         >  7:134382/1‑108 (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                          >  5:18071/2‑130 (MQ=255)
                                                                                                                                         aaTCCACAGTTCCTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGaaaa                          <  6:18071/129‑1 (MQ=255)
                                                                                                                                          aTCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCAc                                            >  7:34152/1‑111 (MQ=255)
                                                                                                                                          aTCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCAc                                            <  8:34152/111‑1 (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                <  7:132648/135‑1 (MQ=255)
                                                                                                                                              acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGaa                >  8:132648/1‑135 (MQ=255)
                                                                                                                                              acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGaa                >  8:137221/1‑135 (MQ=255)
                                                                                                                                              acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATcc  >  7:309062/1‑149 (MQ=255)
                                                                                                                                              acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATcc  >  5:10136/1‑149 (MQ=255)
                                                                                                                                              acaGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATcc  >  4:250569/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCTCGCCAATGCCTTGACTGCGTGGTTAAGTCAATGAAAAAAAACACTATTTTTAGCTGGTGAAAAAATCGTCAGTTTGACCAAAGCCCCCATTTCACGGGGGTTTGAGGGGTATGAGAGGGGGTTGTCCACCGAGTTATCCACAGGAACTGTGGATTGTCCCGAGCGCTTGCTCTAGCACGGGCGTGCCAGCAATGGGATGGACTGTCCGACAATCGCCATGCTCCAACAGACCGGTTGGTCGTCCACCTGAAAACGTCAAGAAAAGATCATTGAATTTTTTTTCAATCC  >  NC_002947/1884252‑1884542

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