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

CCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAA‑AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGCGCGAGGCATAA  >  NC_000913/3511618‑3511882
                                                                                                                                  |                                                                                                                                        
agcatggtcagttcgttttccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATg                                                                                                                               >  1:742813‑M1/122‑139 (MQ=255)
           ttcgttttccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATc                                                                                                                    <  2:1254949‑M1/29‑1 (MQ=255)
           ttcgttttccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATc                                                                                                                    <  1:996653‑M1/29‑1 (MQ=255)
                tttccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATccgccg                                                                                                               >  2:223502‑M1/106‑139 (MQ=255)
                  tccgccaccacgttccagatcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGaa                                                                                                             >  2:93786‑M1/104‑139 (MQ=255)
                                    atcaccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                           >  2:996653‑M1/86‑139 (MQ=255)
                                        ccgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCt                                                                                       <  1:763063‑M1/58‑1 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑A‑ACCGGTCATGATTTTTTGATCCCCCGAACGGTTATCCAcccccaccccctt                                                                                      <  2:761968‑M1/59‑7 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTAt                                                                                                      <  1:1228023‑M1/43‑1 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTAt                                                                                                      >  2:1228023‑M1/81‑123 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCtt                                                                                      <  2:119338‑M1/59‑1 (MQ=255)
                                         cgcatcgacatcgcctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCtt                                                                                      <  2:157946‑M1/59‑1 (MQ=255)
                                           cGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAA‑AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          >  2:468477/1‑118 (MQ=255)
                                           cGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAA‑AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          <  1:468477/118‑1 (MQ=255)
                                                      cctttaacaatgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGAc                                                                         <  2:243554‑M1/72‑1 (MQ=255)
                                                       tccTTTAAGAAAGGCTTCTAAAGTCTCGCGATAATATAGATCGACGCGTTCCACATCCCTCCCGCCAAAAA‑AACACATCCGTCATGATTTTCTTATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          <  1:971759/138‑1 (MQ=255)
                                                                tgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGt                                                                                                                                   <  2:696425‑M1/14‑1 (MQ=255)
                                                                tgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGt                                                                                                                                   >  1:696425‑M1/58‑71 (MQ=255)
                                                                tgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                               <  2:710566‑M1/82‑1 (MQ=255)
                                                                tgcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                               <  1:375069‑M1/82‑1 (MQ=255)
                                                                 gcgttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACg                                                              <  2:876208‑M1/83‑1 (MQ=255)
                                                                   gttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGc                                                            >  1:415909‑M1/55‑139 (MQ=255)
                                                                    ttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATc                                                                                                     <  1:296487‑M1/44‑1 (MQ=255)
                                                                    ttgtaaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATc                                                                                                     >  2:296487‑M1/54‑97 (MQ=255)
                                                                       taaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAAt                                                        <  2:169318‑M1/89‑1 (MQ=255)
                                                                        aaactctcgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATc                                                       >  2:375069‑M1/50‑139 (MQ=255)
                                                                           ctctcgtgataagagagatctactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAAc                                                    >  1:581122‑M1/47‑139 (MQ=255)
                                                                               cgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGCCTCGCCTTTACGGCAAATCAACTGGt                                                <  2:1408308‑M1/97‑1 (MQ=255)
                                                                               cgtgataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGAGCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGt                                                <  1:129425‑M1/97‑1 (MQ=255)
                                                                                   ataagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            >  2:129425‑M1/39‑139 (MQ=255)
                                                                                     aagagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGc                                          >  1:169318‑M1/37‑139 (MQ=255)
                                                                                      tAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAA‑AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGc                                          <  1:1164658/139‑1 (MQ=255)
                                                                                      agagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACt                                                                        >  2:3310‑M1/36‑108 (MQ=255)
                                                                                      agagagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACt                                                                        <  1:3310‑M1/73‑1 (MQ=255)
                                                                                        agagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            <  1:81678‑M1/101‑1 (MQ=255)
                                                                                        agagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            >  2:81678‑M1/34‑134 (MQ=255)
                                                                                        agagatcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCaa                                       <  2:786886‑M1/106‑1 (MQ=255)
                                                                                            atcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            <  2:592415‑M1/101‑1 (MQ=255)
                                                                                            atcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            >  1:592415‑M1/30‑130 (MQ=255)
                                                                                            atcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTa                                   <  2:757096‑M1/110‑1 (MQ=255)
                                                                                            atcgactcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTa                                   <  1:1093479‑M1/110‑1 (MQ=255)
                                                                                                  tcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGtgtg                             <  2:1426347‑M1/116‑1 (MQ=255)
                                                                                                  tcgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGtgtg                             >  1:1120666‑M1/24‑139 (MQ=255)
                                                                                                   cgttccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGc                            >  1:1254949‑M1/23‑139 (MQ=255)
                                                                                                   cTCGTTCCACATCACTACCGCCAAAAATAA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGgtgt                              >  2:738695/1‑139 (MQ=255)
                                                                                                    tCGTTCCACATCACTACCGCCAAAAATAA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGtgtg                             <  1:738695/139‑1 (MQ=255)
                                                                                                      tccacatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGc                         <  1:733835‑M1/120‑1 (MQ=255)
                                                                                                         acatcactaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGcccc                      >  2:763063‑M1/17‑139 (MQ=255)
                                                                                                              actaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGcc                        <  1:460939‑M1/121‑1 (MQ=255)
                                                                                                              actaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGcc                        >  2:460939‑M1/12‑132 (MQ=255)
                                                                                                               ctgccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTcc                <  1:1047162‑M1/129‑1 (MQ=255)
                                                                                                               ctaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          <  1:132783‑M1/71‑1 (MQ=255)
                                                                                                               ctaccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          >  2:132783‑M1/11‑81 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑ACATCGGTCATGATTTTCTGATccgccg                                                                                                               >  1:124748‑M1/10‑44 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑ACATCGGTCATGATTTTCTGATccgccg                                                                                                               <  2:124748‑M1/35‑1 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            >  2:1112728‑M1/10‑111 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                                            <  1:1112728‑M1/102‑1 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          <  1:894194‑M1/39‑1 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          >  2:894194‑M1/10‑48 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                           >  1:386838‑M1/10‑63 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                           <  2:386838‑M1/54‑1 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          >  2:1100343‑M1/10‑80 (MQ=255)
                                                                                                                taccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGa                                                                          <  1:1100343‑M1/71‑1 (MQ=255)
                                                                                                                  ccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          <  1:731785‑M1/39‑1 (MQ=255)
                                                                                                                  ccgccaaAAAAA‑AA‑‑CATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                                          >  2:731785‑M1/8‑46 (MQ=255)
                                                                                                                     gccaAAAAA‑AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAgcgc           <  2:309067‑M1/135‑1 (MQ=255)
                                                                                                                      ccaAAAAA‑AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGcgcg          >  2:1060765‑M1/4‑139 (MQ=255)
                                                                                                                               AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGCGCGAGGCATaa  >  1:757096‑M1/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                        
CCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAA‑AA‑ACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGCGCGAGGCATAA  >  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: ‑