Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,390,631 C→A 9.7% intergenic (+672/‑50) PP_4738 → / → PP_mr56 hypothetical protein/C4 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,390,6310CA 9.7% 148.0 / 3.7 52intergenic (+672/‑50)PP_4738/PP_mr56hypothetical protein/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base A (2/3);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTG  >  NC_002947/5390482‑5390776
                                                                                                                                                     |                                                                                                                                                 
ggCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGc                                                                                                                                                   >  2:422613/1‑150 (MQ=255)
  cGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCtc                                                                                                                                                 <  1:422613/150‑1 (MQ=255)
     cTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCtcatc                                                                                                                                              <  1:164465/150‑1 (MQ=255)
                ccTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCt                                                                                                                                   >  1:1281748/1‑150 (MQ=255)
                 cTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCtcatc                                                                                                                                              >  2:494104/1‑138 (MQ=255)
                 cTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCtcatc                                                                                                                                              <  1:494104/138‑1 (MQ=255)
                 cTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCATGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCgg                                                                                                                                     <  1:1250605/147‑1 (MQ=255)
                 cTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCATGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCgg                                                                                                                                     >  2:1250605/1‑147 (MQ=255)
                  tGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTTCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCttt                                                                                                                                 >  1:299504/1‑150 (MQ=255)
                  tGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCttt                                                                                                                                 >  1:822731/1‑150 (MQ=255)
                  tGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGATGAGCAACGAc                                                                                                                                        <  1:754644/143‑1 (MQ=255)
                  tGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTACGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGATGAGCAACGAc                                                                                                                                        >  2:754644/1‑143 (MQ=255)
                  tGGAATCTGGGCGTTGTGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCttt                                                                                                                                 <  2:968635/150‑1 (MQ=255)
                   ggAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTg                                                                                                                                >  2:471638/1‑150 (MQ=255)
                    gAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGa                                                                                                                               <  1:1144866/150‑1 (MQ=255)
                    gAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCACCATTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGa                                                                                                                               <  1:1397392/150‑1 (MQ=255)
                                    ggTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATac                                                                                                               <  1:471638/150‑1 (MQ=255)
                                           aaaTCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAgcgc                                                                                                        <  1:230036/150‑1 (MQ=255)
                                              tCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCa                                                                                                     <  1:591304/150‑1 (MQ=255)
                                                       aCCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGAc                                                                                                                              >  1:176986/1‑116 (MQ=255)
                                                       aCCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGAc                                                                                                                              <  2:176986/116‑1 (MQ=255)
                                                        ccccTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCa                                                                                           >  2:445044/1‑150 (MQ=255)
                                                           cTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATgg                                                                                        >  1:540415/1‑150 (MQ=255)
                                                             ttCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCg                                                                                      >  2:715530/1‑150 (MQ=255)
                                                              tCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCgg                                                                                     <  1:541521/150‑1 (MQ=255)
                                                                 tGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGGTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTtg                                                                                  <  1:234659/150‑1 (MQ=255)
                                                                      ggCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGATGAGCAACGAc                                                                                                                                        <  1:212525/91‑1 (MQ=255)
                                                                      ggCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGATGAGCAACGAc                                                                                                                                        >  2:212525/1‑91 (MQ=255)
                                                                              ccAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACg                                                                     >  1:128839/1‑150 (MQ=255)
                                                                              ccAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACg                                                                     >  1:544302/1‑150 (MQ=255)
                                                                                                 cTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCt                                                  <  2:540415/150‑1 (MQ=255)
                                                                                                    ttttCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCAt                                                                             <  1:423963/120‑1 (MQ=255)
                                                                                                    ttttCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCAt                                                                             >  2:423963/1‑120 (MQ=255)
                                                                                                    tttgCCGTTGCGGTCATTTTTTCAGCGTTTTAGTCTCGGTGGGTCGTTGCCCATCAACGACCGGCTTTGACCGGCTGCTACCAAACACAGAGCAGTTAGCCACATAGGGCGGGTGTGCATGTGGCACGCTTGCGTGCGCCGGTTTCtgtg                                               <  1:387409/150‑1 (MQ=255)
                                                                                                       tCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTc                                            <  1:445044/150‑1 (MQ=255)
                                                                                                           ttGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCg                                        >  1:597337/1‑150 (MQ=255)
                                                                                                              cGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCACCAACCCCGGGCTTTGACAGGCCGCAAAGATAAACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTc                                     <  1:674517/150‑1 (MQ=255)
                                                                                                                gTCATTTTTTTAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTg                                   >  2:1154400/1‑150 (MQ=255)
                                                                                                                gTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTg                                   >  2:1080712/1‑150 (MQ=255)
                                                                                                                 tCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTAATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGt                                  >  1:871325/1‑150 (MQ=255)
                                                                                                                  cATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTAATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTc                                 <  2:871325/150‑1 (MQ=255)
                                                                                                                   aTTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCa                                >  1:683186/1‑150 (MQ=255)
                                                                                                                    ttttttCAGGGTTTTAGTCTCGGTGGGTCGTTGTTTATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCaa                               <  2:335038/150‑1 (MQ=255)
                                                                                                                    ttttttCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCaa                               <  2:822731/150‑1 (MQ=255)
                                                                                                                        ttCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCa                           <  2:597337/150‑1 (MQ=255)
                                                                                                                            gggTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGca                       >  2:1132397/1‑150 (MQ=255)
                                                                                                                                     tctcGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGcca              >  2:503541/1‑150 (MQ=255)
                                                                                                                                          gccggTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCAtt         >  1:796502/4‑150 (MQ=255)
                                                                                                                                            gggTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTcc       <  2:128839/150‑1 (MQ=255)
                                                                                                                                              gTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGCGCACGCAagcgtgccccat                                                       <  1:594493/100‑6 (MQ=255)
                                                                                                                                              gTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGCGCACGCAagcgtgccccat                                                       >  2:594493/1‑95 (MQ=255)
                                                                                                                                              gTCGTTGATGAGCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCAt                                                                                          <  2:832271/65‑1 (MQ=37)
                                                                                                                                                 gTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTg  >  1:862958/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTG  >  NC_002947/5390482‑5390776

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