Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,273,080 C→T 12.8% intergenic (+45/+246) 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,0800CT12.8% 115.8 / 7.9 51intergenic (+45/+246)PP_2003/PP_2004hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (15/28);  new base T (2/4);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01

ACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGCAAAATGAGCAA  >  NC_002947/2272935‑2273227
                                                                                                                                                 |                                                                                                                                                    
aCAGCCAGTGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGGCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCg                                                                                                                                                    <  3:145038/148‑1 (MQ=38)
aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCg                                                                                                                                                    <  8:233076/148‑1 (MQ=39)
aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCgg                                                                                                                                                   <  2:272920/149‑1 (MQ=255)
  aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCAGAACCGCCCCGATAGCCCCCCGCCCGCGAAAGGGCCggat                                                                                                                                                 >  2:116784/1‑147 (MQ=255)
   gCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCAACCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCATTCGCCGGCAGCCCGCCGAAGGGGCCGTTc                                                                                                                                                >  7:151722/1‑149 (MQ=255)
    ccAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGtt                                                                                                                                                 >  5:59184/1‑147 (MQ=255)
        ggTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCAGAAACGCCCCATTAGCCGGCACCCCGCGAAAGGGCCGGGTCCgggg                                                                                                                                           >  8:157553/1‑147 (MQ=255)
         gTTCGCAGTCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTg                                                                                                                                          >  6:89568/1‑149 (MQ=255)
         gTTCGCAGCCTGTGGCCCCGCCGGCCTAAGCTGCACCCATTCAGGTGCATCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGTCGGTTCAGGcga                                                                                                                                          <  8:151722/149‑3 (MQ=255)
              cAGCCTGTGGCCCCGCCGGCCGACGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAAc                                                                                                                                     <  5:83795/149‑1 (MQ=255)
              cAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGCCCCTTTCGCCGGCCGCCCGCGACAGGGCCGGTTCAGGCTGTaaa                                                                                                                                      >  3:100626/1‑148 (MQ=255)
              cAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAAc                                                                                                                                     <  5:266286/149‑1 (MQ=255)
                gCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACAt                                                                                                                                   <  8:65634/149‑1 (MQ=255)
                gCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCGGAACCCACCCGTTCGCCCGCAGCCCGGGAAAGGGGGCGGTCAGGGGGTAAACag                                                                                                                                   >  5:184588/1‑148 (MQ=255)
                  cTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATcc                                                                                                                                 <  4:64099/149‑1 (MQ=255)
                     tgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATcccc                                                                                                                               <  4:218449/148‑1 (MQ=255)
                     tgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATccc                                                                                                                                <  7:283611/147‑1 (MQ=255)
                       gCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAAACGCACCATTAGCCGGCCGCCCCCGAAAGGGCCGGTTGAGGCTGTAAACATACCCAGt                                                                                                                            >  6:19144/1‑149 (MQ=255)
                       gCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCACTCCCTTTTATGGGTAGCCTGAACCGACCCTTTAGCCGGCACCCCGCCAAAGGGCCGCTTCAGGCTGTAAACATCCCCTGt                                                                                                                            >  6:300373/1‑149 (MQ=255)
                            gccggccgAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAg                                                                                                                       <  5:89568/149‑1 (MQ=255)
                                            cccATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCATTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACt                                                                                                       >  4:50046/1‑149 (MQ=255)
                                              cATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACttt                                                                                                     <  5:19144/149‑1 (MQ=255)
                                              cATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCCACCCTTTCCCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATACCCTGTTAAAGACGCCCCCCGATAACttt                                                                                                     >  2:203605/1‑149 (MQ=255)
                                                   aGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGaa                                                                                                <  5:56961/149‑1 (MQ=255)
                                                        cagcagGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCt                                                                                           <  1:203605/149‑1 (MQ=255)
                                                        cagcagGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCt                                                                                           <  1:116784/149‑1 (MQ=255)
                                                         agcagGATGCCCAGCCTGCCGAGGCCGCCGACCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTg                                                                                          <  4:276845/149‑1 (MQ=255)
                                                         agcagGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGCCCCTTTCGCCGGCCGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTg                                                                                          >  5:102858/1‑149 (MQ=255)
                                                         agcagGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCCACCCTTTAGCCCGCACCCCGCGAAAGGGCCGGTTCACGCTGTAAACATCCCCAGTTACAGACGCCCCCCGATAAATTTCTGAACACCTg                                                                                          >  4:106637/1‑149 (MQ=255)
                                                                  cccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGttt                                                                                  <  1:244628/148‑1 (MQ=255)
                                                                   ccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAc                                                                                <  7:157553/149‑1 (MQ=255)
                                                                             gAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGCCCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAg                                                                      <  4:269928/149‑1 (MQ=255)
                                                                                 ccgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTg                                                                  <  3:50046/149‑1 (MQ=255)
                                                                                  cgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGa                                                                                                            <  5:254225/106‑1 (MQ=255)
                                                                                         ccctccctTTGATGGGTCGCCAGAACCCAACCATTACCCGGCACCCCGCGAAACGGCCGGGTAAGGCAGTTAAAAAACCCAGTTTAAGAAGCCCCCCGAAAACTTTAGTAACACCAGCAACCTTTAACACACATAAACATGtaaaaaaa                                                          >  2:174130/1‑142 (MQ=255)
                                                                                           ctccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATccc                                                        <  7:39131/149‑1 (MQ=255)
                                                                                             ccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCtt                                                      <  4:159312/149‑1 (MQ=255)
                                                                                                   gATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATcc                                                 <  4:100626/148‑1 (MQ=255)
                                                                                                         tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGc                                          >  3:80288/1‑149 (MQ=255)
                                                                                                         tCGCCTGAACCGACCCTTTCGCAGGCAGCCCGCGAAAGGGCGGGTTCAGGCTGTAATCATCCCCTGTTAAAGACGCCCCGCGATAACTTTCTGTACACCTGCTACGTTTACCACTCATAAGCTTGTGATCATACCTTCATCCCACCTGc                                          <  6:234171/149‑1 (MQ=255)
                                                                                                         gcGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGc                                          <  5:62997/148‑1 (MQ=255)
                                                                                                                 aCCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCt                                  >  8:26444/1‑148 (MQ=38)
                                                                                                                    gACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAg                              <  6:59184/149‑1 (MQ=255)
                                                                                                                     gcccTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAc                                                                                <  2:255613/98‑1 (MQ=37)
                                                                                                                     gcccTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAc                                                                                >  1:255613/2‑99 (MQ=37)
                                                                                                                           tCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTACATCCCACCTGCTTAGGCTGAAGGCGAt                         <  1:174130/147‑1 (MQ=255)
                                                                                                                                   aGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAg               <  3:106637/149‑1 (MQ=255)
                                                                                                                                            aaGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCGCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGCAAAATGa      <  6:123313/149‑1 (MQ=255)
                                                                                                                                            aaGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGCAAAATGa      <  5:300373/149‑1 (MQ=255)
                                                                                                                                             aGGGCCGGTTAAGGCTGTAAACATCCACTGTTACAGACGCCCCGCGATAACTTTCTCAACACCTCCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCAACCTGCTTAGGCTGAAGGCGATACCTACAAAGGCAAAATGat     >  4:287746/1‑148 (MQ=255)
                                                                                                                                                gCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCACCAGCATAAGCATCAAGCCGAACCAACAAAAGCaaaaataacaa  >  6:199052/1‑141 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGCAAAATGAGCAA  >  NC_002947/2272935‑2273227

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