Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,956 T→C 8.4% intergenic (+354/+260) 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
*NC_0009133,269,9560TC 8.4% 148.6 / 4.4 62intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base C (2/3);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTT  >  NC_000913/3269811‑3270093
                                                                                                                                                 |                                                                                                                                         
gcATCCGTCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATACTGTCGTGGACAGTGATTAATAATACTGCCCCTCCAAAAGTaa                                                                                                                                         >  7:35197/1‑148 (MQ=21)
    ccGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTGAATCATCTTTCTGCCGCTGCAAAATTTAAAAcc                                                                                                                                    >  7:323413/1‑149 (MQ=21)
        cATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgcc                                                                                                                                <  3:401513/149‑1 (MQ=255)
           gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgccg                                                                                                                               >  7:417427/1‑147 (MQ=255)
           gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgccg                                                                                                                               >  2:394047/1‑147 (MQ=255)
              tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTAGATGACAATAAAACAACTATGAACACCACAAAAATTAAAAAACCTCCGAAGCg                                                                                                                          >  5:33331/1‑149 (MQ=18)
                aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                         <  8:417427/148‑1 (MQ=255)
                    gcACTCAACAATCTGAAACCCGCCGGAACGGGTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGTGATAGTGTCGTGGACCTTCAGGCATCTTTCTGCCACTCCAAAAGCAAAACCCCGCCGACGCGGGtttg                                                                                                                    >  6:301715/1‑148 (MQ=17)
                    gcACTCAACAATCTGAAACCCGCCGGAACGGGTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtttt                                                                                                                    <  1:185048/149‑1 (MQ=255)
                     cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                   <  8:260027/149‑1 (MQ=255)
                      aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                  >  1:56361/1‑149 (MQ=255)
                        tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTGTAGTGTAGAGGACAGTAATTAAGAGTTCTGACCGTCAAAAAGTAAAAAACCGCAGAAGAGGTTTTTTaa                                                                                                                 >  5:370969/1‑147 (MQ=18)
                         caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                               <  5:342483/149‑1 (MQ=255)
                          aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAATACGAAGCCCAAGCGGGttataacca                                                                                                               >  5:234026/1‑141 (MQ=17)
                          aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGACGTAGAGAGTCGTTAATGTTTATAGACATACAAAAGTAAAAACACGCCGTAGAGGCGTTTTAAGt                                                                                                               >  2:150199/1‑148 (MQ=17)
                           acaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAGAGGCTTTGAGCTATCTGCGCCTCCAAAATTAAAAACACGCAGAAGCGGTTTGTTAAGtt                                                                                                              >  5:387623/1‑147 (MQ=17)
                            caaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                  <  8:69401/143‑1 (MQ=255)
                            caaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAACTAAAAACCCGCCGAAGCGGGTTTTTAcaga                                                                                                              >  7:69401/1‑144 (MQ=255)
                                  gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                  <  7:371186/137‑1 (MQ=255)
                                  gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                  >  8:371186/1‑137 (MQ=255)
                                   aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                      <  2:56361/148‑1 (MQ=255)
                                   aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTCACCGATAAGCCGCATTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCCGGTTTTTACGTAAAACAg                                                                                                        <  3:186811/146‑1 (MQ=255)
                                       ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                 <  1:334346/149‑1 (MQ=255)
                                       ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                 <  4:66732/149‑1 (MQ=255)
                                         gccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgcc                                                                                                                                >  1:79741/1‑116 (MQ=255)
                                         gccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgcc                                                                                                                                <  2:79741/116‑1 (MQ=255)
                                         gccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                      >  4:132057/1‑142 (MQ=255)
                                         gccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                      <  3:132057/142‑1 (MQ=255)
                                          ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                                                                              <  4:72483/149‑1 (MQ=255)
                                             gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                                                                                           <  1:394047/149‑1 (MQ=255)
                                             gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                         <  7:453641/119‑1 (MQ=255)
                                             gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                         >  8:453641/1‑119 (MQ=255)
                                              aaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                  <  8:259007/125‑1 (MQ=255)
                                                   tttATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                      <  5:96939/132‑1 (MQ=255)
                                                   tttATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                      >  6:96939/1‑132 (MQ=255)
                                                       ttaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                  >  3:350070/1‑116 (MQ=255)
                                                       ttaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                  <  4:350070/116‑1 (MQ=255)
                                                              cttCAGGTTAATGTGACCGTTATGCCGCTATCTTTTTTTTATAGTGTCGTCGACTGTCTTTCATCTTTCTCCCCCTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                          <  7:380635/147‑1 (MQ=11)
                                                              cATCACGCAAAACTGACCGATAACCCGCTTCCTTATCCGTCTATTTTCGGCTACTGTCATTCATCTATCTGCCACTCAAAAAGTTCAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCtttt                                                                          <  8:323413/149‑1 (MQ=2)
                                                                tCACGTATAACTGACCGATAAGCCGCTTTCTTTTGGCTATTGTGTCGTGGACTGTCTTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                        <  8:35197/149‑1 (MQ=11)
                                                                  tggttacaCTGACCGCTCAGCTGCGTTCTTGTGTGTCTCGTGTCGGGGACTGTGCTTGCTCTTTCTGGCCCTCCAAAAGTATAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGGGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                        <  4:152776/140‑1 (MQ=2)
                                                                  attaaaaaCTTACCGATAAGCCGCATGCATTTGGTTATAGTGTCGACGACACGCCATCATCTATCTGCCACTCCACAAGAACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                                                      <  5:414889/145‑1 (MQ=11)
                                                                    ttataaCTGACCTATACGCATCTTTCTTTTGTGTATAGTCGCGTTGACATGAATTCATCTTTCTGGCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGta                                                                     <  1:150199/144‑1 (MQ=11)
                                                                    gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                      >  5:36299/1‑115 (MQ=39)
                                                                    gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                      <  6:36299/115‑1 (MQ=39)
                                                                      aaaaCTCTCAGATACGACTCTATCTTTTTGCTATTGTGTCGTGGACAGTCTTTCATCTTCCTTCCCCTCAAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                   <  1:389684/148‑1 (MQ=11)
                                                                                                            tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                       <  5:193677/122‑1 (MQ=38)
                                                                                                            tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                       >  6:193677/1‑122 (MQ=38)
                                                                                                                    aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt                                                                            >  7:182934/1‑93 (MQ=37)
                                                                                                                    aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt                                                                            <  8:182934/93‑1 (MQ=37)
                                                                                                                        attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGATTCTGGTTTTTACGTAAAACAGGTGAAACTGACCGAAAAGCCCCTGTCTTTTTGGGTTAGTGTCGTCGTCAGTCATTCATCTTTCTGCCCCTCCaaaa                            <  1:15709/137‑1 (MQ=11)
                                                                                                                        attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCg                                                                                   >  7:248504/1‑82 (MQ=255)
                                                                                                                        attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                       <  5:397610/110‑1 (MQ=17)
                                                                                                                        attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                       >  6:397610/1‑110 (MQ=17)
                                                                                                                            atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                                         >  5:417388/1‑104 (MQ=17)
                                                                                                                                      cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                       >  8:31374/1‑64 (MQ=40)
                                                                                                                                      cccGCCACAAGCACAAACCCGCCGGACCGGGTTTTTACGTACAACAGGGCAAACTTACCGATAAGCCGCTTTCTTTTGGGTCTAGTGTCGTGGACAGGCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt  <  5:301715/149‑1 (MQ=11)
                                                                                                                                       ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                        <  8:17322/78‑1 (MQ=39)
                                                                                                                                       ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                        >  7:17322/1‑78 (MQ=39)
                                                                                                                                       ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                   >  1:302883/1‑83 (MQ=38)
                                                                                                                                       ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                            <  7:386995/122‑1 (MQ=38)
                                                                                                                                                gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtca                                                   >  6:37435/1‑90 (MQ=17)
                                                                                                                                                 |                                                                                                                                         
GCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTT  >  NC_000913/3269811‑3270093

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: