Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,273,066 A→T 15.9% intergenic (+31/+260) 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,0660AT15.9% 100.8 / 12.9 50intergenic (+31/+260)PP_2003/PP_2004hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (16/23);  new base T (4/4);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

GTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGA  >  NC_002947/2272928‑2273203
                                                                                                                                          |                                                                                                                                          
gtggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGaaa                                                                                                                                  <  1:160853/149‑1 (MQ=255)
   gATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCGGGCTGCCGGCGAAAggg                                                                                                                               <  3:116620/149‑1 (MQ=255)
   gATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAgg                                                                                                                                <  2:339705/148‑1 (MQ=255)
       aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCg                                                                                                                            <  6:248963/148‑1 (MQ=39)
       aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCg                                                                                                                            <  6:301691/148‑1 (MQ=255)
       aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTAGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCg                                                                                                                            <  3:67483/148‑1 (MQ=38)
        cAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGt                                                                                                                          <  7:147925/149‑1 (MQ=255)
         aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGtt                                                                                                                         <  7:129729/149‑1 (MQ=255)
           ccAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCa                                                                                                                       <  7:1066/149‑1 (MQ=255)
                gTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTg                                                                                                                  >  8:59202/1‑149 (MQ=255)
                     cAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAAc                                                                                                             <  5:128421/149‑1 (MQ=255)
                           gtgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCCACCCATTCGCCGGCCGCCACCGAAATGGCCAGTTAAGGCTGAAAACAAccc                                                                                                        >  6:147971/1‑148 (MQ=255)
                                  cgccgGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACa                                                                                                <  2:236058/149‑1 (MQ=255)
                                  cgccgGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCGCCCTTTGATGGGTCGCCTGAACCCCCCCTTTCGCCGGCCGCCCGCGAAAGGGCCCGGTCACGCTGTAAACCTccccgggta                                                                                                  >  5:152831/1‑142 (MQ=255)
                                     cggccgAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGAc                                                                                             <  7:239678/149‑1 (MQ=255)
                                                  aCCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAAc                                                                               >  5:111292/1‑149 (MQ=255)
                                                   cccATTCCGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTAGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCAGTAAACATCCCCTGTTACAGACGCCCCGCGAAAACt                                                                              >  2:84496/1‑149 (MQ=255)
                                                      aTTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTc                                                                           <  7:29679/149‑1 (MQ=255)
                                                      aTTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCAGAACCCACCCTTTCGCCGGCAGCCCGCCAAAGGGCCGGTTCAGGCTGTAAAAAACCCCAGTTAAAGACGCCCCGCCAATACttta                                                                           >  4:227614/1‑148 (MQ=255)
                                                       ttCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCt                                                                          <  3:227614/149‑1 (MQ=255)
                                                              gcagcaGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCGTTCGCCGGCAGCCCGCGAAAGGGGCGGGTCAGGCTGTTAACATACCCGGTTTCAGACGCCCCGCGATAACTTTTTGCAAAcc                                                                   >  7:121902/1‑149 (MQ=255)
                                                                    ggATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTAc                                                             <  7:59202/149‑1 (MQ=255)
                                                                      aTGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTAGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATACCCTGTTACAGACGCCCCGCGATAACTTTCTg                                                                         >  8:172404/1‑135 (MQ=255)
                                                                          ccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCAGAACCGACCCATTAGCCGGCCGCCCGCGAAACGGCCGGTTAAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCAGAACCCCTGGAACGTTTaa                                                       >  6:138164/1‑148 (MQ=255)
                                                                           cAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGCAAGGGCCGGTTCAGGCTGGAAACATCCCCTGTTACAGCCGCCCCGCGATAACTTTGTGAACACCTGCTACGTTTAc                                                       >  5:264650/1‑148 (MQ=255)
                                                                                        ccgccgccCCCTCCATTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCTAAAGGGCCGGTTCAGGCTGAAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCtt                                          <  5:147971/148‑1 (MQ=255)
                                                                                         cgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGt                                        >  5:72372/1‑149 (MQ=255)
                                                                                          gccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGt                                        >  6:168447/1‑148 (MQ=255)
                                                                                          gccgccCCCTCCCTTTGATGGGTCGCCGGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtc                                       <  6:152831/149‑1 (MQ=255)
                                                                                              cccccTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtcatca                                   <  5:168447/149‑1 (MQ=255)
                                                                                                ccctccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGCCGCCCCGCGATAACGTTCTGCACACCTGGTACGTTTACCACTGAGAAGCTTGGCataatc                                 >  7:143160/1‑149 (MQ=255)
                                                                                                ccctccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGCATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTcatcat                                  <  1:84496/149‑1 (MQ=255)
                                                                                                 cctccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCGGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATcc                                >  7:163675/1‑149 (MQ=255)
                                                                                                  ctccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCAt                                                <  2:234147/132‑1 (MQ=255)
                                                                                                           aTGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCcg                                                                                       <  3:79745/85‑1 (MQ=255)
                                                                                                              ggTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCTGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCt                   <  8:163675/149‑1 (MQ=255)
                                                                                                               gTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTg                  <  6:72372/149‑1 (MQ=255)
                                                                                                                tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCt                              <  8:137647/136‑1 (MQ=255)
                                                                                                                tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCt                              >  7:137647/1‑136 (MQ=255)
                                                                                                                   ccTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCt                >  4:164544/1‑147 (MQ=255)
                                                                                                                    cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACCCCTGCTACGTTTACCAGTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAg             >  3:275780/1‑149 (MQ=255)
                                                                                                                             cccTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtc                                       <  7:121816/114‑1 (MQ=255)
                                                                                                                             cccTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtc                                       >  8:121816/1‑114 (MQ=255)
                                                                                                                              ccTTTCGCGGGCTGCCGGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACCCCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGc    >  6:247900/1‑148 (MQ=255)
                                                                                                                              ccTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGc                                            >  6:30245/1‑108 (MQ=255)
                                                                                                                              ccTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGc                                            <  5:30245/108‑1 (MQ=255)
                                                                                                                               cTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGCCAGTAAACATCCCATGTTACAGACGCCACGCGATAACTTTCTGAACACCTGCTACGTTTACCAATCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCg   >  6:53286/1‑148 (MQ=255)
                                                                                                                               cTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGa  >  3:80901/1‑149 (MQ=255)
                                                                                                                                tttCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAg             <  3:228163/137‑1 (MQ=255)
                                                                                                                                tttCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAg             >  4:228163/1‑137 (MQ=255)
                                                                                                                                          |                                                                                                                                          
GTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGA  >  NC_002947/2272928‑2273203

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: