Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,466 G→T 17.7% intergenic (+505/+109) 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,4660GT17.7% 1054.0 / 42.7 438intergenic (+505/+109)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (227/109);  new base T (39/40);  total (273/165)
Fisher's exact test for biased strand distribution p-value = 3.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 440 total aligned reads displayed.

AGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCC  >  CP009273/3263399‑3263532
                                                                   |                                                                  
aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACACGCCGAAGCGGGTTTTTACGTAAACCAGGTg                                                                    >  1:32558967/1‑68 (MQ=25)
 gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGa                                                                   >  1:8769626/1‑68 (MQ=35)
 gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGa                                                                   >  1:13546092/1‑68 (MQ=35)
     ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACt                                                               <  1:37180110/68‑1 (MQ=35)
      tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTg                                                              >  1:30276102/1‑68 (MQ=35)
       catATTTCTGCCACTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGa                                                             <  1:6940714/68‑1 (MQ=11)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                             >  1:30722774/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                             >  1:39172167/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                             >  1:18699003/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                             >  1:4755773/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                             >  1:35117956/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                             >  1:34489556/1‑67 (MQ=35)
         tcagtttcacctggtttacgtAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                           <  1:15891377‑M1/47‑1 (MQ=255)
         tCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                           >  1:24971653/1‑68 (MQ=35)
         tCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                           >  1:22359666/1‑68 (MQ=35)
         tCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                           >  1:35589682/1‑68 (MQ=35)
            ttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAt                                                        >  1:7341988/1‑68 (MQ=35)
                gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTTCGTAAACCAGGTGAAACTGACCGATAAg                                                     >  1:16714799/1‑67 (MQ=25)
                gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                                                     >  1:19742564/1‑67 (MQ=35)
                gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                                                     >  1:16463995/1‑67 (MQ=35)
                gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                                                     >  1:15795569/1‑67 (MQ=35)
                gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                                                     >  1:29108414/1‑67 (MQ=35)
                gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                    <  1:14648785/68‑1 (MQ=35)
                gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                    <  1:2329502/68‑1 (MQ=35)
                 ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                   >  1:39405277/1‑68 (MQ=35)
                  cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                   >  1:1010002/1‑67 (MQ=35)
                  cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                   >  1:11617178/1‑67 (MQ=35)
                        aaaaTCATAAACCCGCCTGATCGGTTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGctttct                                            <  1:28471424/68‑1 (MQ=14)
                        aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                            <  1:17140604/68‑1 (MQ=255)
                         aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttctt                                           >  1:30123574/1‑68 (MQ=255)
                         aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGcgttctt                                           >  1:18880148/1‑68 (MQ=255)
                          acgaaaaaaccccctccggcgggtttttACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                          <  1:27702725‑M2/40‑1 (MQ=255)
                          aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGgtttctgt                                          >  1:38673839/1‑66 (MQ=9)
                          aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                          <  1:28004331/68‑1 (MQ=40)
                           aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                         >  1:31651023/1‑68 (MQ=40)
                           aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                         >  1:3125321/1‑68 (MQ=40)
                           aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCtttt                                         >  1:22183306/1‑68 (MQ=255)
                           aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGACCGATAAGCCGCGTTCtttt                                         >  1:36901697/1‑68 (MQ=255)
                            gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                        >  1:17756316/1‑68 (MQ=39)
                            gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                        >  1:37719823/1‑68 (MQ=39)
                            gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                        >  1:14892737/1‑68 (MQ=39)
                             ctcAAACCCGCCGTAACGGTTTTTTACGTACACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                       >  1:23878050‑M1/4‑68 (MQ=255)
                             catAAACCCGCCGTATCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                                       >  1:21312577‑M1/4‑68 (MQ=255)
                             cagAAACCCGCCGTACCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                                       >  1:5725759‑M1/4‑68 (MQ=255)
                             cAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                       >  1:4958956/1‑68 (MQ=37)
                              AAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg                                      >  1:20630587‑M1/1‑68 (MQ=255)
                              AAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg                                      >  1:36713299‑M1/1‑68 (MQ=255)
                              AAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                        >  1:4173455‑M1/1‑66 (MQ=255)
                               tAAACCCGCCGTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:28841590‑M1/2‑68 (MQ=255)
                               tAAACCCGCCGGAACGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:971509‑M1/2‑68 (MQ=255)
                               tAAACCCGCCGCATCGGTTTTTTACGTAAACCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:34887395‑M1/2‑68 (MQ=255)
                               tAAACCCGCCGCATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:18543309‑M1/2‑68 (MQ=255)
                               tAAACCCGCCGCATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:26897709‑M1/2‑68 (MQ=255)
                               tAAACCCGCCGAACCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:26260854‑M1/2‑68 (MQ=255)
                               cAAACCCGCCGTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:26387527‑M1/2‑68 (MQ=255)
                               cAAACCCGCCGTATCGGGTTTTTACGTAAACCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:38210255‑M1/2‑68 (MQ=255)
                               AAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:6282653‑M1/1‑68 (MQ=255)
                               AAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:11952948‑M1/1‑68 (MQ=255)
                               AAAACCCGCCGAACCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                     >  1:10423937‑M1/1‑68 (MQ=255)
                                 AACCCGCCGCAGCGGGTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                   >  1:2026735‑M1/1‑68 (MQ=255)
                                  ACCCGCCGAATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                  >  1:11290440‑M1/1‑68 (MQ=255)
                                   CCCGCCGCACCGGTTTTTTACGTAAAGCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                 >  1:29550620‑M1/1‑68 (MQ=255)
                                    ccgccgCATCGGGTTTTTACGTACATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGc                                <  1:12910723/68‑1 (MQ=11)
                                    CCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGc                                >  1:31903766‑M1/1‑68 (MQ=255)
                                      gccgCATCGGGTTTTTACGTACAGCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                              <  1:30563230/68‑1 (MQ=11)
                                      GCCGAATCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGt                              >  1:21473731‑M1/1‑68 (MQ=255)
                                      GCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGt                              >  1:17267341‑M1/1‑68 (MQ=255)
                                      GCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGt                              >  1:35969315‑M1/1‑68 (MQ=255)
                                      GCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGt                              >  1:25572615‑M1/1‑68 (MQ=255)
                                      GCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGCATAGCGt                              >  1:33324475‑M1/1‑68 (MQ=255)
                                        cgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcg                            <  1:32739496/68‑1 (MQ=40)
                                         gAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                           <  1:1022897/68‑1 (MQ=40)
                                             cgggtttttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGgaa                       >  1:4352511‑M2/10‑67 (MQ=255)
                                              ggttttttACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                      <  1:19462539‑M2/60‑1 (MQ=255)
                                                 tttttACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                   <  1:13838711‑M2/63‑1 (MQ=255)
                                                 tttttACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGt                    <  1:2279178‑M2/62‑1 (MQ=255)
                                                 tttttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGt                    <  1:10014924‑M2/62‑1 (MQ=255)
                                                 tttttACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGt                    <  1:18154878‑M2/62‑1 (MQ=255)
                                                 ttgttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                   <  1:9313805‑M2/63‑1 (MQ=255)
                                                 ttattACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGt                    <  1:35780545‑M2/62‑1 (MQ=255)
                                                 ttattACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGt                    <  1:24628768‑M2/62‑1 (MQ=255)
                                                 ttattACGTAAAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                   <  1:32241331‑M2/63‑1 (MQ=255)
                                                 ttattACGTAAAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                   <  1:15174294‑M2/63‑1 (MQ=255)
                                                  ttttACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                   <  1:28148917‑M2/63‑1 (MQ=255)
                                                  ttttACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                   <  1:7938627‑M2/63‑1 (MQ=255)
                                                  ttttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                   <  1:2092073‑M2/63‑1 (MQ=255)
                                                  ttttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                   <  1:16174871‑M2/63‑1 (MQ=255)
                                                  ttttACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                   <  1:30917533‑M2/63‑1 (MQ=255)
                                                      aCGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGtcattca              >  1:12208724/1‑68 (MQ=255)
                                                      aCGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattca              >  1:33136669/1‑68 (MQ=40)
                                                      aCGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAGCGTCGTGGAAAGtcattca              >  1:5110085/1‑68 (MQ=255)
                                                        gTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATc            <  1:4420745/68‑1 (MQ=255)
                                                        gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGtcattca              >  1:27099094/1‑66 (MQ=255)
                                                        gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCAAGtcaatca              >  1:6699400/1‑66 (MQ=14)
                                                            aCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCCTTCATCttta        >  1:10674622/1‑67 (MQ=255)
                                                             ccAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCATTAATATTTCt       >  1:3368036/1‑68 (MQ=11)
                                                             ccAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCAATCATATTTcg       >  1:8481845/1‑67 (MQ=11)
                                                                  tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGTCATTCATCGTTCAGcacc  >  1:2433999/1‑68 (MQ=11)
                                                                  tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTGTAcccc  >  1:17018838/1‑68 (MQ=255)
                                                                  tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGcc    >  1:37047526/1‑66 (MQ=35)
                                                                   |                                                                  
AGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCC  >  CP009273/3263399‑3263532

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: