Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,556 G→C 5.9% intergenic (+595/+19) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,263,5560GC 5.9% 1096.5 / 8.0 407intergenic (+595/+19)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (202/157);  new base C (3/24);  total (206/201)
Fisher's exact test for biased strand distribution p-value = 3.57e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 413 total aligned reads displayed.

TTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCAT  >  CP009273/3263490‑3263622
                                                                  |                                                                  
ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGc                                                                   <  1:3852784/68‑1 (MQ=35)
 tttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGGCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                                                  >  1:19449065/1‑68 (MQ=255)
    gggTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                               <  1:9673075/68‑1 (MQ=35)
    gggTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                               <  1:13117911/68‑1 (MQ=35)
      gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                             >  1:10425674/1‑68 (MQ=35)
      gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                             >  1:8277261/1‑68 (MQ=35)
      gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                             >  1:27378698/1‑68 (MQ=35)
      gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                              <  1:2614371/67‑1 (MQ=35)
      gTATAGCGTCGTGGACAGTCATTAATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                             >  1:9325045/1‑68 (MQ=25)
        ataGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                           >  1:1788070/1‑68 (MQ=35)
        ataGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                           >  1:14569663/1‑68 (MQ=35)
        ataGCGTCGTGGACAGTCATTCATCATTCTGCCCCTCCAAAAGCAAAAACCCGCCCAAGCGGGttttt                                                           >  1:17617948/1‑68 (MQ=11)
          aGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:8620676/68‑1 (MQ=35)
           gCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                        >  1:4279349/1‑68 (MQ=33)
           gCGTCGTGGACAGTCATTAATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                        >  1:26544061/1‑68 (MQ=255)
            cgtcgtGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGt                                                       >  1:30039195/1‑68 (MQ=31)
                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCTGGTTATTACGTAAATCa                                                 >  1:4009208/1‑68 (MQ=255)
                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCa                                                 >  1:20161470/1‑68 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                               >  1:23459505/1‑68 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                               >  1:21865549/1‑68 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                               >  1:16710805/1‑68 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                               >  1:10553144/1‑68 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                               >  1:12725620/1‑68 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:22349665/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:31402429/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:22071466/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:22858423/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:9187057/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:2458706/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:28455307/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:20771888/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:20459827/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                                <  1:12886101/67‑1 (MQ=35)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                               >  1:22600131/1‑68 (MQ=25)
                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCAGCCGAAGCGGGGTTTTACGTAAATAAgg                                               >  1:1654736/1‑68 (MQ=11)
                    aCAGTCATTCATCTTTCTGCCACTCCAAAAGAAAAAACCCGCCGAACCGGGTTTTTACGTAAATCAgg                                               >  1:25724472/1‑68 (MQ=11)
                    aCAGTCATTCATCTTTCTGACCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                               >  1:14486103/1‑68 (MQ=25)
                      aGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTg                                             >  1:10708836/1‑68 (MQ=25)
                       gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGa                                            >  1:5127536/1‑68 (MQ=35)
                       gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGa                                            >  1:4572592/1‑68 (MQ=35)
                       gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGa                                            >  1:2440927/1‑68 (MQ=35)
                         cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                              <  1:14667612/64‑1 (MQ=35)
                         cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                              <  1:18294476/64‑1 (MQ=35)
                         cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGaa                                           <  1:17282755/67‑1 (MQ=35)
                            tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTg                                       >  1:9050192/1‑68 (MQ=35)
                              atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                      >  1:27617756/1‑67 (MQ=35)
                              atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                      >  1:29438750/1‑67 (MQ=35)
                              atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                      >  1:10098377/1‑67 (MQ=35)
                              atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                      >  1:23300427/1‑67 (MQ=35)
                              atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                      >  1:1157026/1‑67 (MQ=35)
                               tCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAcc                                    >  1:5866307/1‑68 (MQ=35)
                                   tCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATa                                <  1:5451380/68‑1 (MQ=35)
                                      gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAg                              >  1:30351905/1‑67 (MQ=35)
                                      gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAg                              >  1:7103932/1‑67 (MQ=35)
                                      gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                             <  1:8169924/68‑1 (MQ=35)
                                      gCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAATGAGGTGAAACTGACCGATAAg                              >  1:21610035/1‑67 (MQ=11)
                                       ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGcc                            >  1:12188675/1‑68 (MQ=35)
                                        cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCg                           >  1:18775042/1‑68 (MQ=35)
                                         ccTCCAAAAGCACAAACCCGCCGAACCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                          <  1:3270368/68‑1 (MQ=255)
                                         ccTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                          <  1:31737192/68‑1 (MQ=255)
                                         ccTCAAAAATCTAAAACCCGCCGTACCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                          <  1:23830739/68‑1 (MQ=11)
                                         acTCCAAAAGAAAAAACCCGCCGTATCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                          <  1:3727059/67‑1 (MQ=11)
                                          cTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCggg                         <  1:3541466/68‑1 (MQ=39)
                                           tCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGt                        <  1:24315493/68‑1 (MQ=40)
                                             cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTc                      <  1:6624095/68‑1 (MQ=255)
                                              aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCt                     <  1:10994428/68‑1 (MQ=255)
                                                aatcAAAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGt                   <  1:2705528/65‑1 (MQ=17)
                                                aaGCATAAACCCGCCGTACCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                   <  1:11939162/68‑1 (MQ=17)
                                                aaGCACAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGt                   <  1:30553744/68‑1 (MQ=255)
                                                aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGt                        >  1:13408470/1‑63 (MQ=40)
                                                aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                   <  1:2400511/68‑1 (MQ=255)
                                                aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                   <  1:11299460/68‑1 (MQ=255)
                                                aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                   <  1:14109998/68‑1 (MQ=255)
                                                aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                   <  1:14915190/68‑1 (MQ=255)
                                                aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                   <  1:19070798/68‑1 (MQ=255)
                                                aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                   <  1:12421849/68‑1 (MQ=255)
                                                aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                   <  1:31045748/68‑1 (MQ=255)
                                                aaGCAAAAACCCGCCGAACCGGTTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                   <  1:28847047/68‑1 (MQ=255)
                                                  gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCg                 >  1:16356005/1‑68 (MQ=255)
                                                    aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGt                <  1:28093357/67‑1 (MQ=255)
                                                      aaaCCCGCCGTACCGGGTTTTTACGTAAATCAGGTCAAACTGACCGATAAGCCGGGTTCTGTCGTGGa             <  1:15509570/68‑1 (MQ=18)
                                                        aCCCGCCGAAGCGGTTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACa           <  1:21192393/68‑1 (MQ=255)
                                                         cccGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAg          >  1:24750937/1‑68 (MQ=255)
                                                          ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt         <  1:6268899/68‑1 (MQ=255)
                                                          ccgccgAACCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt         <  1:24984491/68‑1 (MQ=21)
                                                           cgccgAAGCGTGTTTTTACGTAAAACCGGTTAAACTGACAGATAAGCCGTGTTCTGTCGTGGACAGtc        <  1:7747869/68‑1 (MQ=17)
                                                           cgccgAAGCGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtc        <  1:7523900/68‑1 (MQ=21)
                                                           cgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtc        <  1:5621403/68‑1 (MQ=255)
                                                           cgccgAACCGGTTTATTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtc        <  1:15304648/68‑1 (MQ=21)
                                                           cgccgAACCGGGTTATTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt         <  1:22169819/67‑1 (MQ=255)
                                                            gccgCAGCGGTTTTTTACGTAAATCAGGTCAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtca       <  1:13545782/68‑1 (MQ=21)
                                                            gccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtca       <  1:25572064/68‑1 (MQ=255)
                                                            gccgAAGCGGGTTATTACGTACATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtca       <  1:15040068/68‑1 (MQ=255)
                                                              cgAAGCGTTTTTTTACGTAAATCCCGTGAACCAGACCGATAAGCCGGGTTCTGTCGTGGACAGtcatt     <  1:31472238/68‑1 (MQ=25)
                                                               gAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattc    >  1:26771172/1‑68 (MQ=255)
                                                                 cgCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATACGCCGGGTTCTGTCGTGGACAGtcattcat  <  1:25449683/67‑1 (MQ=255)
                                                                 aGCGTGTTTTTACGTAAATCTGGTCAAACTGACAGATAATCCGGGTTCTGTCGTGGACAGtcattca   <  1:13824765/67‑1 (MQ=25)
                                                                 aGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattca   <  1:17056436/67‑1 (MQ=255)
                                                                  gCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattca   >  1:30692000/1‑66 (MQ=255)
                                                                  gCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattca   >  1:15920505/1‑66 (MQ=255)
                                                                  |                                                                  
TTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCAT  >  CP009273/3263490‑3263622

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