Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,654,144 T→A 7.5% noncoding (290/304 nt) ryfA → small regulatory RNA RyfA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,654,1440TA 7.5% 140.6 / 4.4 54noncoding (290/304 nt)ryfAsmall regulatory RNA RyfA
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base A (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

TGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGA  >  NC_000913/2654002‑2654284
                                                                                                                                              |                                                                                                                                             
tGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCttt                                                                                                                                               <  1:525568/143‑1 (MQ=255)
    gATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggc                                                                                                                                           <  1:155348/143‑1 (MQ=255)
       aCGGTGCTTTTGTACCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGggcggc                                                                                                                                        <  1:131162/143‑1 (MQ=255)
          gTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGCGCCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCatt                                                                                                                                     <  2:59715/143‑2 (MQ=255)
          gTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGACTCGCCGGCGACAAGTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCata                                                                                                                                     <  1:396839/143‑1 (MQ=255)
                 tGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGCCAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGACTGCCGCCCTTTGCGCGTCATATCTTCGa                                                                                                                             >  1:7617/1‑143 (MQ=255)
                 tGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGa                                                                                                                             <  1:581027/143‑1 (MQ=255)
                           ttATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGAtt                                                                                                                   >  1:465909/1‑143 (MQ=255)
                           ttATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGAtt                                                                                                                   >  2:617646/1‑143 (MQ=255)
                            tATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTc                                                                                                                  >  1:528694/1‑143 (MQ=255)
                                gTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGtt                                                                                                              >  2:416733/1‑143 (MQ=255)
                                 tttCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGCATGCCGCCCTTTGGGCGGCATATTTTCGATTATCCGCTTCTGttt                                                                                                             >  1:100696/1‑143 (MQ=255)
                                  ttCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTa                                                                                                            <  2:559270/143‑1 (MQ=255)
                                    cTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaaa                                                                                                          <  1:282006/143‑1 (MQ=255)
                                     tGGTGCGCTGTTAAGCGAGGTAGATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAg                                                                                                         <  2:471925/143‑1 (MQ=255)
                                        tGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCa                                                                                                      <  1:315756/143‑1 (MQ=255)
                                        tGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCAGGTTTTAGATTAGCCGGTTCTGTTTACAGTCa                                                                                                      >  2:753269/1‑143 (MQ=255)
                                             tGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaa                                                                                                 >  2:124864/1‑143 (MQ=255)
                                               ttAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaaa                                                                                               <  2:458973/143‑1 (MQ=255)
                                                tggccGAAGTAAATAATGACAGGAGTCGCTCCGGCGACAATGTACTGGTGGTTAACAACGTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaaaa                                                                                              <  2:772426/140‑1 (MQ=255)
                                                 aaCCGAGGTAAATAATAACCGGAGTCTTTCCGGCGCCACTGCACTTGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCTCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaaaaa                                                                                             <  2:187711/143‑1 (MQ=255)
                                                  aCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCGAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTGGAGATTATCCGATTCTGCTTAAAGTCACGAAAAAAac                                                                                            >  2:915362/1‑143 (MQ=255)
                                                  aCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAac                                                                                            <  2:416483/143‑1 (MQ=255)
                                                        gTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAGAGTCACGCAAAACACCAcccc                                                                                      >  1:52675/1‑143 (MQ=255)
                                                                taaCCGGAGTCTCTGCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGGCCTGTGGGCGGCATATTTTAGATGAGCCGATTCTGTATAAAGTCACGCAAAAAACCACCCCAGCGACGt                                                                              <  2:362458/143‑1 (MQ=255)
                                                                   ccGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGACTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCa                                                                           >  1:814643/1‑143 (MQ=255)
                                                                     ggAGTTTCTCCGGCGACATTTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTGTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATa                                                                         <  2:855346/143‑1 (MQ=255)
                                                                           tctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCAAAGAATgg                                                                   >  2:75473/1‑143 (MQ=255)
                                                                            ctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATGCGATTCTGTTTAAAGTCACGCAAAAAACCACCGCAGCGACGTTCATAGAATGgc                                                                  >  1:19824/1‑143 (MQ=255)
                                                                            ctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGgc                                                                  >  1:710497/1‑143 (MQ=255)
                                                                                       aTTTACTGGTGTTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAACAAACCACCCCAGCGACGTTCATAGCATGGCGCCATCAGTTc                                                       >  1:254703/1‑143 (MQ=255)
                                                                                         ttACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATCTTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCACTAAGTTCgg                                                     >  2:245260/1‑143 (MQ=255)
                                                                                         tgaCTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCACGAATGCAGCACTTTGGGCGGCATATTTTACACTAACCGTTTCTGTTCAAACTCACGCAAAAAACCTCCCCAGCGACGTTCATAGAATGGCGCAATATGTTCgg                                                     <  2:403868/141‑1 (MQ=255)
                                                                                           aCTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTTGATTATCCTATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGta                                                   <  1:305227/143‑1 (MQ=255)
                                                                                             tggtggTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTAGGAGAAATATGGTGAGATATCCCATCCTGTTTAGAGTCGCGCAAAGAACCACCCCAGCGACGTGCATAGAATGGCGCAATATGTTCgcgact                                                 >  2:647528/1‑138 (MQ=255)
                                                                                             tggtggTTAACAACCATCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaat                                                 >  2:772008/1‑143 (MQ=255)
                                                                                                tggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATaaa                                              <  1:124864/143‑1 (MQ=255)
                                                                                                    taccaacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTgg                                          <  1:213889/140‑1 (MQ=255)
                                                                                                    tAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTgg                                          >  2:520927/1‑143 (MQ=255)
                                                                                                    tAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCAATGAATTGCGCCATATGTTCGGTAATAACGtcg                                          >  1:707305/1‑141 (MQ=255)
                                                                                                            ttCAGAGCAGCAAGTAAGCACGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTACATAGAATGGCGCAATATATTCGGTAATAAAGTGACGAATccc                                  >  2:268485/1‑143 (MQ=255)
                                                                                                               agagCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATCATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCttt                               <  1:546272/143‑1 (MQ=255)
                                                                                                                     gcaAGTAAGCCCGAATGCCGCCCTTTGGGCCGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTcccc                         <  2:7617/143‑1 (MQ=255)
                                                                                                                      caAGTAAGCCAGAATGCCGCCCTTTGGGCGGCATATTTTAGGTTATCAGATTCTGTTTAAAGTCACGCAAAACCCCACACCAGCGCCGTTCATAGAATGGCGCAAGACGTTCGGTAATACAGCGGCTAAGTCCtcttgcgtct                        >  2:106081/1‑137 (MQ=255)
                                                                                                                        aGTACGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCttt                      <  1:450200/143‑1 (MQ=255)
                                                                                                                        aGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTGAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTTATTCCTTTTTCGCCttt                      >  1:546197/1‑143 (MQ=255)
                                                                                                                         gTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGACCCTGCCTAAAGTCACGCAAAAACCCACCCCAGCGACGTCCATACAACGGCGCAATATGTTCGGTAATCAAGTGGCTCATTCCTCTCCCCCCttct                     >  1:408649/1‑141 (MQ=255)
                                                                                                                         gTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGACTATACGATGCTGTCTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttt                     <  1:662931/143‑1 (MQ=255)
                                                                                                                            aGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATGATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTc                  <  1:617646/143‑1 (MQ=255)
                                                                                                                             taaaataTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCa                 >  2:592259/7‑143 (MQ=255)
                                                                                                                               ccGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCAcc               <  2:814643/143‑1 (MQ=255)
                                                                                                                                  taTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGa            >  1:69505/2‑143 (MQ=255)
                                                                                                                                  aaTGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGa            <  1:520927/143‑1 (MQ=255)
                                                                                                                                     gccgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACaa         >  2:282472/1‑143 (MQ=255)
                                                                                                                                        gccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTGCCCTGACAAata      >  2:667591/1‑143 (MQ=255)
                                                                                                                                            tttGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGa  >  1:628234/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGA  >  NC_000913/2654002‑2654284

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: