Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,511,746 (A)8→7 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,7460A.81.8% 157.4 / 50.2 66coding (599/906 nt)yhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (6/6);  new base . (24/30);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGCGCGAGGCATAA  >  NC_000913/3511618‑3511882
                                                                                                                                  |                                                                                                                                        
agcatggtcagttcgttttccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATg                                                                                                                               >  1:742813‑M1/122‑139 (MQ=255)
           ttcgttttccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATc                                                                                                                    <  2:1254949‑M1/29‑1 (MQ=255)
           ttcgttttccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATc                                                                                                                    <  1:996653‑M1/29‑1 (MQ=255)
                tttccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                               >  2:223502‑M1/106‑139 (MQ=255)
                  tccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGaa                                                                                                             >  2:93786‑M1/104‑139 (MQ=255)
                                    atcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                           >  2:996653‑M1/86‑139 (MQ=255)
                                        ccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCt                                                                                       <  1:763063‑M1/58‑1 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAAAACCGGTCATGATTTTTTGATCCCCCGAACGGTTATCCAcccccaccccctt                                                                                      <  2:761968‑M1/59‑7 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTAt                                                                                                      <  1:1228023‑M1/43‑1 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTAt                                                                                                      >  2:1228023‑M1/81‑123 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCtt                                                                                      <  2:119338‑M1/59‑1 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCtt                                                                                      <  2:157946‑M1/59‑1 (MQ=255)
                                           cGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          >  2:468477/1‑118 (MQ=255)
                                           cGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          <  1:468477/118‑1 (MQ=255)
                                                      cctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGAc                                                                         <  2:243554‑M1/72‑1 (MQ=255)
                                                       tccTTTAAGAAAGGCTTCTAAAGTCTCGCGATAATATAGATCGACGCGTTCCACATCCCTCCCGCCAAAAAAACACATCCGTCATGATTTTCTTATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          <  1:971759/138‑1 (MQ=255)
                                                                tgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGt                                                                                                                                   <  2:696425‑M1/14‑1 (MQ=255)
                                                                tgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGt                                                                                                                                   >  1:696425‑M1/58‑71 (MQ=255)
                                                                tgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                               <  2:710566‑M1/82‑1 (MQ=255)
                                                                tgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                               <  1:375069‑M1/82‑1 (MQ=255)
                                                                 gcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACg                                                              <  2:876208‑M1/83‑1 (MQ=255)
                                                                   gttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGc                                                            >  1:415909‑M1/55‑139 (MQ=255)
                                                                    ttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATc                                                                                                     <  1:296487‑M1/44‑1 (MQ=255)
                                                                    ttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATc                                                                                                     >  2:296487‑M1/54‑97 (MQ=255)
                                                                       taaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAAt                                                        <  2:169318‑M1/89‑1 (MQ=255)
                                                                        aaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATc                                                       >  2:375069‑M1/50‑139 (MQ=255)
                                                                           ctctcgtgataagagagatctactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAAc                                                    >  1:581122‑M1/47‑139 (MQ=255)
                                                                               cgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGCCTCGCCTTTACGGCAAATCAACTGGt                                                <  2:1408308‑M1/97‑1 (MQ=255)
                                                                               cgtgataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGAGCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGt                                                <  1:129425‑M1/97‑1 (MQ=255)
                                                                                   ataagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            >  2:129425‑M1/39‑139 (MQ=255)
                                                                                     aagagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGc                                          >  1:169318‑M1/37‑139 (MQ=255)
                                                                                      tAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGc                                          <  1:1164658/139‑1 (MQ=255)
                                                                                      agagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACt                                                                        >  2:3310‑M1/36‑108 (MQ=255)
                                                                                      agagagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACt                                                                        <  1:3310‑M1/73‑1 (MQ=255)
                                                                                        agagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            <  1:81678‑M1/101‑1 (MQ=255)
                                                                                        agagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            >  2:81678‑M1/34‑134 (MQ=255)
                                                                                        agagatcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCaa                                       <  2:786886‑M1/106‑1 (MQ=255)
                                                                                            atcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            <  2:592415‑M1/101‑1 (MQ=255)
                                                                                            atcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            >  1:592415‑M1/30‑130 (MQ=255)
                                                                                            atcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTa                                   <  2:757096‑M1/110‑1 (MQ=255)
                                                                                            atcgactcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTa                                   <  1:1093479‑M1/110‑1 (MQ=255)
                                                                                                  tcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGtgtg                             <  2:1426347‑M1/116‑1 (MQ=255)
                                                                                                  tcgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGtgtg                             >  1:1120666‑M1/24‑139 (MQ=255)
                                                                                                   cgttccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGc                            >  1:1254949‑M1/23‑139 (MQ=255)
                                                                                                   cTCGTTCCACATCACTACCGCCAAAAATAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGgtgt                              >  2:738695/1‑139 (MQ=255)
                                                                                                    tCGTTCCACATCACTACCGCCAAAAATAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGtgtg                             <  1:738695/139‑1 (MQ=255)
                                                                                                      tccacatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGc                         <  1:733835‑M1/120‑1 (MQ=255)
                                                                                                         acatcactaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGcccc                      >  2:763063‑M1/17‑139 (MQ=255)
                                                                                                              actaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGcc                        <  1:460939‑M1/121‑1 (MQ=255)
                                                                                                              actaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGcc                        >  2:460939‑M1/12‑132 (MQ=255)
                                                                                                               ctgccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTcc                <  1:1047162‑M1/129‑1 (MQ=255)
                                                                                                               ctaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          <  1:132783‑M1/71‑1 (MQ=255)
                                                                                                               ctaccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          >  2:132783‑M1/11‑81 (MQ=255)
                                                                                                                taccgccaaAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                               >  1:124748‑M1/10‑44 (MQ=255)
                                                                                                                taccgccaaAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                               <  2:124748‑M1/35‑1 (MQ=255)
                                                                                                                taccgccaaAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            >  2:1112728‑M1/10‑111 (MQ=255)
                                                                                                                taccgccaaAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            <  1:1112728‑M1/102‑1 (MQ=255)
                                                                                                                taccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          <  1:894194‑M1/39‑1 (MQ=255)
                                                                                                                taccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          >  2:894194‑M1/10‑48 (MQ=255)
                                                                                                                taccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                           >  1:386838‑M1/10‑63 (MQ=255)
                                                                                                                taccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                           <  2:386838‑M1/54‑1 (MQ=255)
                                                                                                                taccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          >  2:1100343‑M1/10‑80 (MQ=255)
                                                                                                                taccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          <  1:1100343‑M1/71‑1 (MQ=255)
                                                                                                                  ccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          <  1:731785‑M1/39‑1 (MQ=255)
                                                                                                                  ccgccaaAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          >  2:731785‑M1/8‑46 (MQ=255)
                                                                                                                     gccaAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAgcgc           <  2:309067‑M1/135‑1 (MQ=255)
                                                                                                                      ccaAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGcgcg          >  2:1060765‑M1/4‑139 (MQ=255)
                                                                                                                               AAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGCGCGAGGCATaa  >  1:757096‑M1/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                        
CCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGCGCGAGGCATAA  >  NC_000913/3511618‑3511882

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