Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,560,141 C→G 7.0% C20W (TGC→TGG)  PP_2243 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,560,1410CG 7.0% 181.5 / 5.3 57C20W (TGC→TGG) PP_2243hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base G (2/2);  total (31/26)
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.08e-01

TGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGTTCAGGCGTGCAACTGCTGA  >  NC_002947/2559997‑2560283
                                                                                                                                                |                                                                                                                                              
tgcGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGt                                                                                                                                           <  2:193085/150‑1 (MQ=255)
     tGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAg                                                                                                                                      <  1:868855/150‑1 (MQ=255)
      gAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAga                                                                                                                                     >  1:440870/1‑149 (MQ=255)
                ttCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGaaa                                                                                                                                       >  1:131729/1‑138 (MQ=255)
                ttCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGaaa                                                                                                                                       <  2:131729/138‑1 (MQ=255)
                  cTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGAc                                                                                                                         >  1:224945/1‑150 (MQ=255)
                     aTCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCt                                                                                                                      <  1:357039/150‑1 (MQ=255)
                      tCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTg                                                                                                                     <  2:328990/150‑1 (MQ=255)
                      tCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTg                                                                                                                     >  1:1038968/1‑150 (MQ=255)
                       cccAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTTTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAGATGCCTGGTGAAAGCCATCCGCCTGACCCTGa                                                                                                                    >  1:1597892/1‑150 (MQ=255)
                        ccAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGAt                                                                                                                   >  2:425387/1‑150 (MQ=255)
                         cAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATc                                                                                                                  >  2:1307107/1‑150 (MQ=255)
                          aaGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCt                                                                                                                 >  2:1422933/1‑150 (MQ=255)
                                    cTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCt                                                                                                       <  2:885463/150‑1 (MQ=255)
                                        ggCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGAcc                                                                                                   <  1:608677/150‑1 (MQ=255)
                                            cAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAgg                                                                                               >  2:413528/1‑150 (MQ=255)
                                              aTCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCATCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGcc                                                                                             >  1:1271011/1‑150 (MQ=255)
                                              aTCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGcc                                                                                             >  1:913142/1‑150 (MQ=255)
                                                     ttCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCagaaga                                                                                      >  1:1063103/1‑150 (MQ=255)
                                                        aaTCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAGAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACg                                                                                   <  2:1315043/150‑1 (MQ=255)
                                                          tCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCagaa                                                                                        >  2:22530/1‑143 (MQ=255)
                                                          tCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCagaa                                                                                        <  1:22530/143‑1 (MQ=255)
                                                          tCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACggg                                                                                 >  2:237617/1‑150 (MQ=255)
                                                                  cgAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCa                                                                         <  1:425387/150‑1 (MQ=255)
                                                                              gacATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCAc                                                                        >  2:655989/1‑139 (MQ=255)
                                                                              gacATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCAc                                                                        <  1:655989/139‑1 (MQ=255)
                                                                                   tCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCAt                                                                                                            <  1:1647247/98‑1 (MQ=255)
                                                                                   tCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCAt                                                                                                            >  2:1647247/1‑98 (MQ=255)
                                                                                     ttGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATc                                                      <  1:413528/150‑1 (MQ=255)
                                                                                       gCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCt                                                    >  2:361018/1‑150 (MQ=255)
                                                                                            cAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGc                                               <  1:1307107/150‑1 (MQ=255)
                                                                                             aaaCGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCt                                              >  2:251088/1‑150 (MQ=255)
                                                                                                    tcGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAgaagaa                                                                                     <  1:33531/104‑1 (MQ=255)
                                                                                                    tcGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAgaagaa                                                                                     >  2:33531/1‑104 (MQ=255)
                                                                                                          cGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCAtcagtca                                 <  2:1271011/150‑1 (MQ=255)
                                                                                                          cGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCAtcagtca                                 <  2:1038968/150‑1 (MQ=255)
                                                                                                               tAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGc                            >  2:1604951/1‑150 (MQ=255)
                                                                                                               tAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGc                            >  1:415285/1‑150 (MQ=255)
                                                                                                                 aaCCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCg                          >  1:1520717/1‑150 (MQ=255)
                                                                                                                      tGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGt                     >  1:1414787/1‑150 (MQ=255)
                                                                                                                         aCCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACTCTGTTCTTCTTTTCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGTTCa                  >  2:1480031/1‑150 (MQ=255)
                                                                                                                           tttagcacGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGTTCAgg                >  2:1556663/8‑150 (MQ=255)
                                                                                                                            ccAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGTAGAACGGTCGCCTGCACCATACGCATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGTTCAGGc               <  2:581858/150‑1 (MQ=255)
                                                                                                                             cAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATAc                                                                   >  2:453760/1‑97 (MQ=255)
                                                                                                                             cAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATAc                                                                   <  1:453760/97‑1 (MQ=255)
                                                                                                                                cccGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATc                                                      >  2:1223415/1‑107 (MQ=255)
                                                                                                                                cccGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATc                                                      <  1:1223415/107‑1 (MQ=255)
                                                                                                                                cccGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGTTCAGGCGTGc           <  2:1520717/150‑1 (MQ=255)
                                                                                                                                 ccGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGc                                           >  1:612747/1‑117 (MQ=255)
                                                                                                                                 ccGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGc                                           <  2:612747/117‑1 (MQ=255)
                                                                                                                                     tGCGAGAAATGGCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGgc                                                                                >  2:309604/1‑76 (MQ=255)
                                                                                                                                     tGCGAGAAATGGCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGgc                                                                                >  2:309782/1‑76 (MQ=255)
                                                                                                                                     tGCGAGAAATGGCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGgc                                                                                <  1:309604/76‑1 (MQ=255)
                                                                                                                                     tGCGAGAAATGGCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGgc                                                                                <  1:309782/76‑1 (MQ=255)
                                                                                                                                       cGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGTTCAGGCGTGCAActgct    <  2:1414787/150‑1 (MQ=255)
                                                                                                                                        gagaAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGTTCAGGCGTGCAActgctg   >  1:711068/1‑150 (MQ=255)
                                                                                                                                         agaAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGTTCAGGCGTGCAACTGCTGa  <  1:237617/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCCTGACCCAGGCCCAGAAGAACGGGCGCCTGCACCATACGGATGACGGCATCCTGCCGCTAGCCCATCAGTCAGACGCCGACCGTTCAGGCGTGCAACTGCTGA  >  NC_002947/2559997‑2560283

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