Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,273,070 C→G 12.3% intergenic (+35/+256) PP_2003 → / ← PP_2004 hypothetical protein/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,273,0700CG12.3% 128.8 / 10.1 49intergenic (+35/+256)PP_2003/PP_2004hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (19/24);  new base G (4/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

CGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGG  >  NC_002947/2272923‑2273215
                                                                                                                                                   |                                                                                                                                                   
cGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCcg                                                                                                                                                    >  6:291509/1‑149 (MQ=255)
         aTGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGc                                                                                                                                           >  4:107416/1‑149 (MQ=255)
         aTGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGc                                                                                                                                           <  2:306738/149‑1 (MQ=255)
            aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCg                                                                                                                                         <  1:373370/148‑1 (MQ=255)
              aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAgggt                                                                                                                                           <  1:159892/144‑2 (MQ=255)
              aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAgggt                                                                                                                                           >  2:159892/1‑143 (MQ=255)
              aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGATCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGtt                                                                                                                                      >  5:369812/1‑149 (MQ=255)
                ccAGGGTTCGCAGCCTGTGGCCCCGCCGGGCGAAGCTGCACCCATTCTGGTGCAGCAGGATGCCCAGCCTGGCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACGCTTTCGGCGGGAGGCCGGGTAAGGGCCGGttga                                                                                                                                    <  8:26395/149‑3 (MQ=255)
                ccAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGtt                                                                                                                                      >  2:254113/1‑147 (MQ=255)
                ccAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑a                                                                                                                                  >  6:184555/1‑149 (MQ=255)
                   gggTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCGGAAACCCCCCATTGGCCGGCCGCCCGCGCAAGGGGCGGGTG‑‑Aggg                                                                                                                               >  1:375123/1‑148 (MQ=255)
                      cgcGCAGCCTGTGGCCCCGCCGGCCGAAGCGGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGGCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGt                                                                                                                            <  6:190371/147‑1 (MQ=255)
                        cGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTa                                                                                                                           <  1:14858/148‑1 (MQ=255)
                         gCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCGGGCTGCCGGCGAAAgg                                                                                                                                             >  5:265975/1‑131 (MQ=255)
                         gCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCGGGCTGCCGGCGAAAgg                                                                                                                                             <  6:265975/131‑1 (MQ=255)
                             ccTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCCGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTACACATc                                                                                                                     >  3:103167/1‑149 (MQ=255)
                              cTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATcc                                                                                                                    >  2:106750/1‑149 (MQ=255)
                               tgtgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAAACCACCCATTGGCCCGCAGCCCGCCAAATGGGCCGTTA‑‑AGGCTGTaaaaaaacc                                                                                                                   >  6:232259/1‑143 (MQ=255)
                                 tgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCt                                                                                                                 <  3:107416/149‑1 (MQ=255)
                                      ccgccgGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTAc                                                                                                            <  5:184555/149‑1 (MQ=255)
                                       cgccgGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACa                                                                                                           <  6:65494/149‑1 (MQ=255)
                                       cgccgGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACa                                                                                                           <  5:291509/149‑1 (MQ=255)
                                                     gCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAAACCCCCCATTCGCCGGCCGCCCCCGAAAGGGCCCGGTA‑‑AGGCTGTAAACCTGCCCGGTTACAGACGCCCCGCGAt                                                                                              >  3:28493/1‑148 (MQ=255)
                                                            ttCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTAGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGAAGCCCCGCGCTTACTTTCt                                                                                      >  6:202791/1‑149 (MQ=255)
                                                              cAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTg                                                                                     <  7:38508/148‑1 (MQ=255)
                                                              cAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTg                                                                                     <  1:106750/148‑1 (MQ=255)
                                                                       cagGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCt                                                                           <  4:301894/149‑1 (MQ=255)
                                                                               ccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATccc                                                                                                                   <  2:273973/101‑1 (MQ=255)
                                                                               ccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATccc                                                                                                                   >  1:273973/1‑101 (MQ=255)
                                                                                 aGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCa                                                                 <  2:375123/149‑1 (MQ=255)
                                                                                    cggccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTc                                                              <  3:233335/149‑1 (MQ=255)
                                                                                          aGGCCGCCGCCCCCTCCCTTTGATGGATCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGc                                                        <  6:369812/149‑1 (MQ=255)
                                                                                            gccgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCtt                                                      <  4:28493/149‑1 (MQ=255)
                                                                                                          cctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCt                                          >  1:298336/1‑149 (MQ=255)
                                                                                                           ctTTGATGGGTCGCCTGAACCGACCCTTTAGCCGGCATCCCGCGAAACGGCCGGTTA‑‑AGGCTGTAAACATACCCTGTTACAGACGCCCCCCGATAAATTTATTAACACCAGCTAAGTTTACCACTCATAACCATGTCATCATACCtta                                        >  4:87660/1‑147 (MQ=255)
                                                                                                             ttGATGGGTCGCCTGAAACGCCACTTTCGCCGGCAGCCCGCCAAAGGGCCCGTTC‑‑AGGCTGTAAACAACCCCTGTTACAGACGCCCCCCCAAAACATTCTGAACACCTGCTACGTTTAACAAACAAAATCTTGTCATCATACCTTCAta                                     >  2:122801/1‑148 (MQ=255)
                                                                                                                            aaCCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGc                       <  2:264501/148‑1 (MQ=255)
                                                                                                                              ccGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGAAAACATCCCCTGTTACAGACGCCCCGCGATATCTTTCTGAACACCTGCTACGTTTAACAGTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTg                     >  4:216672/1‑148 (MQ=255)
                                                                                                                               cGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGaa                   <  5:232259/149‑1 (MQ=255)
                                                                                                                                ggcccTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAg                  >  7:247531/3‑149 (MQ=255)
                                                                                                                                gACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAg                  >  4:144873/1‑149 (MQ=255)
                                                                                                                                    cTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGa              >  5:357662/1‑149 (MQ=255)
                                                                                                                                     tttCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAg                  <  4:161518/144‑1 (MQ=255)
                                                                                                                                     tttCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAg                  >  3:161518/1‑144 (MQ=255)
                                                                                                                                     tttCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGa              >  2:117978/1‑148 (MQ=255)
                                                                                                                                          ccGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTa        <  3:144873/149‑1 (MQ=255)
                                                                                                                                            ggCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCtt                                                      >  2:194705/1‑101 (MQ=255)
                                                                                                                                            ggCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCtt                                                      <  1:194705/101‑1 (MQ=255)
                                                                                                                                               aGCCCGCGAAATGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAg   <  3:87660/149‑1 (MQ=255)
                                                                                                                                                gCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAgg  <  2:298336/149‑1 (MQ=255)
                                                                                                                                                  ccGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAc                                                                   <  7:220063/82‑1 (MQ=255)
                                                                                                                                                  ccGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAc                                                                   >  8:220063/1‑82 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTC‑‑AGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGG  >  NC_002947/2272923‑2273215

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