Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,554 A→C 6.2% intergenic (+593/+21) 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,5540AC 6.2% 1062.1 / 8.1 410intergenic (+593/+21)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (211/153);  new base C (5/25);  total (217/193)
Fisher's exact test for biased strand distribution p-value = 1.20e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 413 total aligned reads displayed.

TCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCAT  >  CP009273/3263488‑3263618
                                                                  |                                                                
     tGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                            <  1:32990946/68‑1 (MQ=35)
      gggTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                           <  1:25615246/68‑1 (MQ=35)
       ggTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                          >  1:8453248/1‑68 (MQ=35)
        gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                         >  1:10256096/1‑68 (MQ=35)
        gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                         >  1:20705168/1‑68 (MQ=35)
        gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                         >  1:16174873/1‑68 (MQ=35)
        gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                         >  1:6628732/1‑68 (MQ=35)
        gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                         >  1:8284218/1‑68 (MQ=35)
        gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                          <  1:28021922/67‑1 (MQ=35)
        gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                          <  1:16538830/67‑1 (MQ=35)
        gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                          <  1:23632950/67‑1 (MQ=35)
        gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                          <  1:30734473/67‑1 (MQ=35)
          ataGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGttttt                                                       >  1:10730764/1‑68 (MQ=25)
          ataGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCCAAGCGGGttttt                                                       >  1:31180916/1‑68 (MQ=25)
            aGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                     <  1:33530490/68‑1 (MQ=35)
            aGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                     <  1:23681298/68‑1 (MQ=35)
             gCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                    >  1:19771093/1‑68 (MQ=33)
             gCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                    >  1:15506450/1‑68 (MQ=33)
             gCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                    >  1:586314/1‑68 (MQ=33)
             gCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                    >  1:19059361/1‑68 (MQ=33)
             gCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                    >  1:7293061/1‑68 (MQ=33)
              cgtcgtGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGt                                                   >  1:17621704/1‑68 (MQ=31)
                    ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCa                                             >  1:30061812/1‑68 (MQ=35)
                    ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCa                                             <  1:4389181/68‑1 (MQ=35)
                    ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAATCCCGCCGAAGCGGGTTTTTACGTAAATCa                                             >  1:25780636/1‑68 (MQ=255)
                     taCATTCATTCATCTTTATGCCCCTCAAAAATCAAAAACCCGCCGAAGCGGGTTTTACGTAAATCAg                                            <  1:6254089/66‑1 (MQ=11)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           >  1:26499182/1‑68 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           >  1:28388662/1‑68 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           >  1:28476545/1‑68 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           >  1:3208833/1‑68 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           >  1:7591374/1‑68 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           >  1:36177370/1‑68 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           >  1:4001229/1‑68 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           >  1:9392762/1‑68 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           >  1:171738/1‑68 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                            <  1:20825156/67‑1 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                            <  1:18084795/67‑1 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                            <  1:1154292/67‑1 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                            <  1:7980419/67‑1 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                            <  1:10828788/67‑1 (MQ=35)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAAACGCCCAAGCGGGTTTTTACGTAAATCAgg                                           >  1:13019315/1‑68 (MQ=11)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGCAGCGGGGTTTTAAGTAAATCAgg                                           >  1:1409681/1‑68 (MQ=11)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAATCGGGTTTTTACGTAAATCAgg                                           >  1:4174834/1‑68 (MQ=255)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGCTTTTTACGTAAATCAgg                                           >  1:22793925/1‑68 (MQ=255)
                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGCTGTTTACGTAAATCAgg                                           >  1:20372196/1‑68 (MQ=11)
                      aCAGTCAATCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                            <  1:11724430/67‑1 (MQ=25)
                       cactcaatcaaCTTTCTGCCACTCAAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                           <  1:19996886/64‑1 (MQ=11)
                         gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTgc                                        >  1:2531937/1‑67 (MQ=25)
                         gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGa                                        >  1:22948597/1‑68 (MQ=35)
                         gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGa                                        <  1:21079887/68‑1 (MQ=35)
                           cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGaa                                       <  1:28284087/67‑1 (MQ=35)
                              tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTg                                   >  1:28815704/1‑68 (MQ=35)
                                atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                  >  1:36987532/1‑67 (MQ=35)
                                atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                  >  1:35448039/1‑67 (MQ=35)
                                atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                  >  1:18738106/1‑67 (MQ=35)
                                     tCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATa                            >  1:13240531/1‑68 (MQ=35)
                                        gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATGAGGTGAAACTGACCGATAAg                          >  1:34223116/1‑67 (MQ=25)
                                        gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAg                          >  1:26186194/1‑67 (MQ=35)
                                        gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAg                          >  1:15087491/1‑67 (MQ=35)
                                        gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAg                          >  1:294768/1‑67 (MQ=35)
                                        gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAg                          >  1:30287682/1‑67 (MQ=35)
                                          cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCg                       >  1:651346/1‑68 (MQ=35)
                                          cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCg                       >  1:10542677/1‑68 (MQ=35)
                                             tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGt                    <  1:21716765/68‑1 (MQ=25)
                                             tCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGt                    >  1:2729713/1‑68 (MQ=40)
                                                 aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGCCCGATAAGCCGGGt                    >  1:21407101/1‑64 (MQ=255)
                                                  taGCAAATACCAGACGAAGCGGGTTTTTACTTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt               <  1:27660880/67‑1 (MQ=17)
                                                  aatcAAAAACCCGCCGTAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGt               <  1:26939068/65‑1 (MQ=17)
                                                  aaGCATAACCCCGCAGAAGCGTGTTTTTACGAAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt               <  1:5576101/68‑1 (MQ=17)
                                                  aaGCATAAACCCGCCGTAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGt               <  1:34654330/68‑1 (MQ=17)
                                                  aaGCACAAACCCGCCGGACCGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGt               <  1:35047151/68‑1 (MQ=14)
                                                  aaGCAAAAACCCGCCGCATCGGGTTATTACGTAAATCAGGTCAAACTGACCGATAAGCCGGGTTCTGt               <  1:14386415/68‑1 (MQ=17)
                                                  aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCATGTGAAACTGACCGATAAGCCGGGTTCTGt               <  1:14989542/68‑1 (MQ=255)
                                                  aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt               <  1:12629878/68‑1 (MQ=255)
                                                  aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt               <  1:10040046/68‑1 (MQ=255)
                                                  aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt               <  1:12158918/68‑1 (MQ=255)
                                                  aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt               <  1:20521706/68‑1 (MQ=255)
                                                   aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTc              >  1:11213238/1‑68 (MQ=255)
                                                     cttaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGt            <  1:36837614/65‑1 (MQ=255)
                                                     ctaaaaCCCGCCTGATCGGTTTATTACGTAAATCATGTTAAACTGACCGATAAGCCGGGTTCTGTCGt            <  1:6941487/66‑1 (MQ=14)
                                                      taaaaCCCGCCGAAGCGGGTTATTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGt            <  1:23013936/66‑1 (MQ=255)
                                                      ataaaCCCGCCGCAGCGGTTTATTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg           <  1:22698722/66‑1 (MQ=18)
                                                         aaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGAc        >  1:32473981/1‑68 (MQ=255)
                                                         aaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGAc        >  1:14687440/1‑68 (MQ=255)
                                                          aCCCGCCGTACCGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGTCGTGGACa       <  1:19314619/68‑1 (MQ=17)
                                                          aCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACa       <  1:25149614/68‑1 (MQ=255)
                                                          aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGGTTCTGTGGTGGaaa       >  1:5176023/1‑66 (MQ=14)
                                                           cccGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAg      >  1:9158635/1‑68 (MQ=255)
                                                            ccgcggATGCGGGTTTTTACGTAAATCAGGAGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt     <  1:29137000/68‑1 (MQ=21)
                                                            ccgccgCATCGGTTTTTTACGTACATCAGGTTAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt     <  1:23478817/68‑1 (MQ=18)
                                                            ccgccgAATCGGTTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt     <  1:23919809/68‑1 (MQ=255)
                                                            ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt     <  1:33795876/68‑1 (MQ=255)
                                                            ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt     <  1:10945356/68‑1 (MQ=255)
                                                             cgccgATGCGGGTTTTTTCGTATATCAGGTGAAACTGAGCGATAAGGCGTCTTCTGTCGTGGACAGtc    <  1:882799/68‑1 (MQ=18)
                                                             cgccgAAGCGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtc    <  1:18475311/68‑1 (MQ=21)
                                                             cgccgAAGCGGGTTTTTACGTAATTCAGGTTAAACTGACCGATAATCCGGGTTCTGTCGTGGACAGtc    <  1:16849087/68‑1 (MQ=21)
                                                             cgccgAAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtc    <  1:24728788/68‑1 (MQ=255)
                                                             cgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtc    <  1:30472758/68‑1 (MQ=255)
                                                             cgacgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAATCCGGGTTCTGTCGTGGACAGtc    <  1:24960006/65‑1 (MQ=255)
                                                              gccgAACCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtca   <  1:15925384/68‑1 (MQ=21)
                                                                  |                                                                
TCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCAT  >  CP009273/3263488‑3263618

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