Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,357:1 +T 17.9% intergenic (+151/+21) ydbL → / ← feaR DUF1318 family protein/transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,3571.T17.9% 200.4 / 49.0 67intergenic (+151/+21)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base . (28/27);  new base T (8/4);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCA  >  NC_000913/1446258‑1446456
                                                                                                    |                                                                                                   
ttttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCa                                                                                                      <  2:462558/100‑1 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                    >  2:75185/1‑100 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                    <  2:171291/100‑1 (MQ=255)
  ttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGt                                                                                                   <  1:981515/100‑1 (MQ=255)
    ccTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTAc                                                                                                 >  1:2110550/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                                >  1:1131065/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                                >  1:1131321/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                                >  2:2298229/1‑100 (MQ=255)
      tGAAACAAGGTTAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACAt                                                                                               <  1:1956654/100‑1 (MQ=255)
        aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGc                                                                                             <  2:2199655/100‑1 (MQ=255)
              ggTGAATATTCAAAAACTCCTGTCGAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                       <  1:1770857/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                       <  1:1420608/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                       <  1:764464/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGcc                                                                                       <  2:1648916/100‑1 (MQ=255)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGcc                                                                                       <  1:1029318/100‑1 (MQ=255)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGcc                                                                                       >  2:490321/1‑100 (MQ=255)
               gTGAATATTCAAAAACTCCTGTCAAATTGTCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTTATTTTCGTTGCATGTACTGCAACg                                                                                         <  2:1360148/97‑1 (MQ=255)
                tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCAGATTCAAATGCCTGTTTTAACTAATTGGCGTTGAAGTAGATGCAAAGCCa                                                                                      <  2:1556902/99‑1 (MQ=255)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAt                                                                                    >  1:1222556/1‑100 (MQ=255)
                  aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAtt                                                                                   <  1:1047426/100‑1 (MQ=255)
                     atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagt                                                                                >  2:2324804/1‑100 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                               >  1:182841/1‑100 (MQ=255)
                         aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGc                                                                            >  2:38250/1‑100 (MQ=255)
                          aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCg                                                                           >  1:1941508/1‑100 (MQ=255)
                               tCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAtt                                                                      <  1:933169/100‑1 (MQ=255)
                                     cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                >  1:1092690/1‑100 (MQ=255)
                                         ttGCATTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGAt                                                             <  2:2086057/99‑1 (MQ=255)
                                          tGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATAc                                                           >  2:860062/1‑100 (MQ=255)
                                            ccgtttGCGCTGAAAGATGCATAGTCAAATATCTGTTTTAAGTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATAc                                                           <  2:2298718/95‑1 (MQ=39)
                                             cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                        <  1:1177822/100‑1 (MQ=255)
                                               tttGCCCTGAAAAATGCATAGTCAAATATCTGGTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGcc                                                      <  2:1676643/100‑1 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                 >  1:2013069/1‑100 (MQ=255)
                                                     cTGATAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGTAATGATAATTTGTTACCGGAATTTACGTCGATACTCGCCTgg                                                   >  2:2005219/1‑97 (MQ=255)
                                                     cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                                >  2:1713256/1‑100 (MQ=255)
                                                      tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                               >  1:2357601/1‑100 (MQ=255)
                                                          aaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATcc                                           <  1:1706154/100‑1 (MQ=255)
                                                           aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATccc                                          >  2:296619/1‑100 (MQ=255)
                                                                   gTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                  >  1:719648/1‑100 (MQ=255)
                                                                   gTCAAATATCTGGTTTAACTAATTGGCGTTGCAGTACATTCAACGCAAATTAGTTAGCGGAATTTACGTAGATACTCGCCTGGCGTCATCCCAAAGc                                     >  2:1618782/1‑97 (MQ=255)
                                                                    tCAAATATCTGTTTTAATTAATTGGCGTTGCAGTATATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                 <  1:549768/100‑1 (MQ=255)
                                                                    tCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                 >  1:1532270/1‑100 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGTGTCATCCCAAAGCGTTGc                                <  2:979857/100‑1 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGACAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                <  1:947591/100‑1 (MQ=255)
                                                                      aaaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCt                               <  1:1326919/100‑1 (MQ=255)
                                                                       aaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                              >  2:1016717/1‑100 (MQ=255)
                                                                        atatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTa                             <  1:1713269/100‑1 (MQ=255)
                                                                          atCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                           <  2:1176477/100‑1 (MQ=255)
                                                                           tCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                          >  1:155466/1‑100 (MQ=255)
                                                                            cTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATa                         <  2:2275074/100‑1 (MQ=255)
                                                                               ttttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCg                      <  1:1087741/100‑1 (MQ=255)
                                                                               tcttccgatcttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCg                      <  1:1426869/90‑1 (MQ=255)
                                                                                  tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                   >  1:1962191/1‑100 (MQ=255)
                                                                                   aaCTAATTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                  >  2:579385/1‑100 (MQ=37)
                                                                                   aaCTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                  <  2:1707881/100‑1 (MQ=255)
                                                                                     gtAATTGGCGTTGCATGTTCTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACACGCCTGGCTTCATCCCAAAGCATAGCGTAAACACCGTTGaa                 >  2:1650212/2‑99 (MQ=21)
                                                                                     cTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTCAATACCGTTGata                <  1:1614564/100‑3 (MQ=255)
                                                                                       aaTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa              >  1:5907/1‑100 (MQ=38)
                                                                                       aaTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa              <  1:1536325/100‑1 (MQ=255)
                                                                                       aaTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa              >  2:2328569/1‑100 (MQ=255)
                                                                                       aaTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa              >  2:1771632/1‑100 (MQ=255)
                                                                                         ttGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg            >  1:2314046/1‑100 (MQ=38)
                                                                                         ttGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg            >  1:1715407/1‑100 (MQ=255)
                                                                                          tGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa           >  2:88613/1‑100 (MQ=38)
                                                                                           ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc          >  1:1212155/1‑100 (MQ=38)
                                                                                           ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc          >  1:2000323/1‑100 (MQ=38)
                                                                                           ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc          >  2:2206036/1‑100 (MQ=38)
                                                                                           ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc          >  1:177207/1‑100 (MQ=255)
                                                                                            gCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAct         >  2:378959/1‑100 (MQ=255)
                                                                                             cGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc        <  1:2159664/100‑1 (MQ=255)
                                                                                                  cAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc   <  2:27801/100‑1 (MQ=255)
                                                                                                   aGTACATGCAAGGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc   <  2:2202392/99‑1 (MQ=255)
                                                                                                   aGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc   >  2:674521/1‑99 (MQ=255)
                                                                                                   aGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa  <  1:1402752/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCA  >  NC_000913/1446258‑1446456

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