Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,511,746 (A)8→9 coding (599/906 nt) yhfZ ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,511,7381.A98.2% 181.1 / ‑3.7 57coding (607/906 nt)yhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (29/27);  total (30/27)
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.94e-01

CGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑AAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTC  >  NC_000913/3511604‑3511867
                                                                                                                                       |                                                                                                                                    
cgTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑aaaa                                                                                                                                 >  2:850245/1‑139 (MQ=255)
                aGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATg                                                                                                                 >  2:1513690/1‑139 (MQ=255)
                aGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATg                                                                                                                 <  1:1513690/139‑1 (MQ=255)
                  cATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑AAAAAAACACATCGGTCATGAt                                                                                                              >  1:431363/1‑139 (MQ=255)
                       tCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCt                                                                                                         >  2:12873/1‑139 (MQ=255)
                                tttCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                <  2:431363/139‑1 (MQ=255)
                                tttCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                <  1:850245/139‑1 (MQ=255)
                                     gCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                           <  2:78289/139‑1 (MQ=255)
                                                    aTCACCGCATCGACATCGCTTTTAACAATGCGTTGTAAACTCTTGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑aaaaaaa                                                                                                                              <  1:584769/90‑1 (MQ=255)
                                                    aTCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑aaaaaaa                                                                                                                              >  2:584769/1‑90 (MQ=255)
                                                       aCCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACgcgc                                                                         >  2:989534/1‑139 (MQ=255)
                                                         cGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCtt                                                                       >  1:1069516/1‑139 (MQ=255)
                                                           cATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCCCTCCCCCCAAAAAAAAACACCGGCCTGATTTTTTTGACCCCCCGACCGGTTTTCAACCCCCACCCCCTTCa                                                                     >  2:664277/1‑139 (MQ=255)
                                                                        tttAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                           <  2:615191/104‑1 (MQ=255)
                                                                        tttAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                           >  1:615191/1‑104 (MQ=255)
                                                                        tttAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTc                                                        <  2:1069516/139‑1 (MQ=255)
                                                                                 gCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTc                                                                                                                    >  2:454771/1‑71 (MQ=255)
                                                                                 gCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTc                                                                                                                    <  1:454771/71‑1 (MQ=255)
                                                                                                 tGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                            >  1:1158698/1‑94 (MQ=255)
                                                                                                 tGATAAGAGAGATCGACTCGTTCCACAGCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                            <  2:1158698/94‑1 (MQ=255)
                                                                                                   aTAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                             <  1:12873/139‑1 (MQ=255)
                                                                                                      agagagATCGACTCGTTCCACATCACTACCGCCAAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                            >  2:1355123/1‑90 (MQ=255)
                                                                                                      agagagATCGACTCGTTCCACATCACTACCGCCAAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                            <  1:1355123/90‑1 (MQ=255)
                                                                                                      agagagATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGcc                          >  2:31024/1‑139 (MQ=255)
                                                                                                           gaTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGt                     >  1:1226747/1‑139 (MQ=255)
                                                                                                            aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGt                     >  2:1295924/1‑139 (MQ=255)
                                                                                                            aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccg                                                                                                   >  1:608725/1‑60 (MQ=255)
                                                                                                            aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccg                                                                                                   <  2:608725/60‑1 (MQ=255)
                                                                                                            aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                           <  2:335455/68‑1 (MQ=255)
                                                                                                            aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                           >  1:335455/1‑68 (MQ=255)
                                                                                                            aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                            >  1:716828/1‑83 (MQ=255)
                                                                                                            aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                            <  2:716828/83‑1 (MQ=255)
                                                                                                            aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACt                                                         >  1:1480412/1‑102 (MQ=255)
                                                                                                            aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACt                                                         <  2:1480412/102‑1 (MQ=255)
                                                                                                             tCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGc                                             <  1:1547043/113‑1 (MQ=255)
                                                                                                             tCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGc                                             >  2:1547043/1‑113 (MQ=255)
                                                                                                               gACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCaa                                           <  2:1507826/113‑1 (MQ=255)
                                                                                                               gACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCaa                                           >  1:1507826/1‑113 (MQ=255)
                                                                                                                aCTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGAt                                                                                                      <  1:804553/53‑1 (MQ=255)
                                                                                                                aCTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGAt                                                                                                      >  2:804553/1‑53 (MQ=255)
                                                                                                                  tCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAAc                                     >  1:727974/1‑116 (MQ=255)
                                                                                                                  tCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAAc                                     <  2:727974/116‑1 (MQ=255)
                                                                                                                  tCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCa                         >  2:440628/1‑128 (MQ=255)
                                                                                                                  tCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCa                         <  1:440628/128‑1 (MQ=255)
                                                                                                                  tCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTCTACGGCAAATGAACTGGTGCTCGCCAACGTAGGtgtg              >  1:1063448/1‑139 (MQ=255)
                                                                                                                    gTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTc                                                        <  1:1173282/95‑1 (MQ=255)
                                                                                                                    gTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTc                                                        >  2:1173282/1‑95 (MQ=255)
                                                                                                                           aTCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                            <  2:390411/68‑1 (MQ=255)
                                                                                                                           aTCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                            >  1:390411/1‑68 (MQ=255)
                                                                                                                               cTACCGCCAAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTc  <  1:1295924/139‑1 (MQ=255)
                                                                                                                                tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                           >  2:499625/1‑48 (MQ=255)
                                                                                                                                tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                           <  1:499625/48‑1 (MQ=255)
                                                                                                                                tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTTCGGACTCGCCTTTACGGCAAATCAACTGGTg                                >  1:1397884/1‑107 (MQ=255)
                                                                                                                                tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGc                                             <  2:231433/94‑1 (MQ=255)
                                                                                                                                tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGc                                             >  1:231433/1‑94 (MQ=255)
                                                                                                                                tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTg                                <  2:1397884/107‑1 (MQ=255)
                                                                                                                                tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGtgtg              <  2:1113607/125‑1 (MQ=255)
                                                                                                                                tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGtgtg              >  1:1113607/1‑125 (MQ=255)
                                                                                                                                  ccGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                <  2:348006/89‑1 (MQ=255)
                                                                                                                                  ccGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                >  1:348006/1‑89 (MQ=255)
                                                                                                                                       |                                                                                                                                    
CGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑AAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTC  >  NC_000913/3511604‑3511867

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: