Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,390,480 C→A 9.2% intergenic (+521/‑201) 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,4800CA 9.2% 169.0 / 4.0 54intergenic (+521/‑201)PP_4738/PP_mr56hypothetical protein/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base A (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 = 9.67e-01

CACAACTCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTC  >  NC_002947/5390333‑5390618
                                                                                                                                                   |                                                                                                                                          
cacaACTCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCgg                                                                                                                                          <  1:1193551/150‑1 (MQ=255)
cacaACTCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCgg                                                                                                                                          <  2:163379/150‑1 (MQ=255)
  caACTCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGc                                                                                                                                        <  2:1079149/150‑1 (MQ=255)
      tCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCACATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTcc                                                                                                                                    >  2:359334/1‑149 (MQ=255)
      tCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCt                                                                                                                                    >  1:447342/1‑150 (MQ=255)
       cATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCGAGGGCGTCtt                                                                                                                                   <  2:1102602/150‑1 (MQ=255)
        aTCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCttt                                                                                                                                  >  2:1122173/1‑150 (MQ=255)
         tCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCtttt                                                                                                                                 >  1:664416/1‑150 (MQ=255)
          ccAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCCGATAAGAAGACGCCCTCGGGCGTCttttt                                                                                                                                <  1:359334/150‑1 (MQ=255)
          ccAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCttttt                                                                                                                                <  2:13899/150‑1 (MQ=255)
                   ccggCCCGTCGGGCGGTTGGGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGAGGAGTACTGGCGGACGGGTGGTCAAATTCCCCCGGCTCCACCACATTCACAGAAGAACCCGCCCTCGGGCGTCTTTTTTTGTGCCTg                                                                                                                       <  1:28394/147‑1 (MQ=255)
                    cTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTgg                                                                                                                      <  2:664416/150‑1 (MQ=255)
                      gCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGaa                                                                                                                    <  2:313271/150‑1 (MQ=255)
                             gggCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGgc                                                                                                             >  1:912225/1‑150 (MQ=255)
                                    tGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGCCCCCCTCGCGCGCCTTTATTTGTGCCTGGGATTTGGGCGCTGCgg                                                                                                      <  1:228980/150‑1 (MQ=255)
                                    tGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCgg                                                                                                      <  1:754513/150‑1 (MQ=255)
                                         aTTAAATCAGTAGACACGGCTAAGCATCTAGTACCGACAGCGGAGTACTGGCGGACGGGGGGTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGATGACGCCCTCGAGCGTCTTTTTTTCTGCCTGGAATTCGGGCGCCGGGGGGCgg                                                                                                 >  2:457013/1‑149 (MQ=255)
                                                   tAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGAt                                                                                       >  1:64959/1‑150 (MQ=255)
                                                       cacGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTAAAATCCCCCTGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTaa                                                                                   <  1:282102/150‑1 (MQ=255)
                                                        acGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAc                                                                                  >  2:1185041/1‑150 (MQ=255)
                                                        acGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAc                                                                                  <  2:64959/150‑1 (MQ=255)
                                                        acGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAc                                                                                  >  1:736036/1‑150 (MQ=255)
                                                         cGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAcc                                                                                 <  1:251172/150‑1 (MQ=255)
                                                            cTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCt                                                                              <  1:1122173/150‑1 (MQ=255)
                                                               aGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTc                                                                           <  2:912225/150‑1 (MQ=255)
                                                                        taCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGc                                                                  >  2:875214/1‑150 (MQ=255)
                                                                          ccGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGt                                                                >  1:357181/1‑150 (MQ=255)
                                                                            gACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTggg                                                                                                              >  2:295295/1‑102 (MQ=255)
                                                                            gACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTggg                                                                                                              <  1:295295/102‑1 (MQ=255)
                                                                                  ggAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCaaa                                                        <  2:357181/150‑1 (MQ=255)
                                                                                         tGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATActg                                                 >  1:1006362/1‑150 (MQ=255)
                                                                                              gACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATActgcctgc                                            <  1:875214/150‑1 (MQ=255)
                                                                                                    ggTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGt                                      >  1:389108/1‑150 (MQ=255)
                                                                                                     gTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGtt                                     <  2:436101/150‑1 (MQ=255)
                                                                                                      ttCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATGGGATTTAACCCTTTTCTTGTAGGGCGTTTCCCAAAGCTACTGCCCGCGCCTGttt                                    >  2:774701/1‑150 (MQ=255)
                                                                                                            tCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGtt                              >  1:273125/1‑150 (MQ=255)
                                                                                                              cccccGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATActgcctg                                             <  2:656751/133‑1 (MQ=255)
                                                                                                              cccccGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATActgcctg                                             >  1:656751/1‑133 (MQ=255)
                                                                                                                 ccGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGtt                              <  2:171937/145‑1 (MQ=255)
                                                                                                                 ccGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGtt                              >  1:171937/1‑145 (MQ=255)
                                                                                                                        caccaAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCAttttt                  >  2:1130971/1‑150 (MQ=255)
                                                                                                                                 cAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGttt         >  1:1226474/1‑150 (MQ=255)
                                                                                                                                  aaaGAAGAAGACGCCTTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTCTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGtttt        <  1:457013/150‑1 (MQ=255)
                                                                                                                                   aagaagaagACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTa       <  2:273125/150‑1 (MQ=255)
                                                                                                                                    aaaaaaagaCGCCCGAGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCg                           <  2:1006098/124‑1 (MQ=255)
                                                                                                                                    aaaaaaagaCGCCCGAGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCg                           >  1:1006098/6‑129 (MQ=255)
                                                                                                                                     gaagaagaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGt     >  1:98851/1‑150 (MQ=255)
                                                                                                                                      aaaaagaCGCCCGAGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTc                                  <  2:510803/117‑1 (MQ=255)
                                                                                                                                      aaaaagaCGCCCGAGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTc                                  >  1:510803/4‑120 (MQ=255)
                                                                                                                                        gaagaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGtctc  >  1:855539/1‑150 (MQ=255)
                                                                                                                                                ccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATActgc                                                >  2:393090/1‑96 (MQ=255)
                                                                                                                                                ccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATActgc                                                <  1:393090/96‑1 (MQ=255)
                                                                                                                                                ccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGtt                              <  1:21850/114‑1 (MQ=255)
                                                                                                                                                ccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGtt                              >  2:21850/1‑114 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CACAACTCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTC  >  NC_002947/5390333‑5390618

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