Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,093,773	A→C	K92Q (AAG→CAG)	yqgF →	putative Holliday junction resolvase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,093,773	0	A	C	100.0%	185.2 / NA	48	K92Q (AAG→CAG)	yqgF	putative Holliday junction resolvase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/28);  total (20/28)

GAACATTATCGAGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCG  >  NC_000913/3093625‑3093920                                                                                                                                                     |                                                                                                                                                    gAACATTATCGAGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAg                                                                                                                                                   <  1:187532/151‑1 (MQ=255)        aTCGAGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATg                                                                                                                                            <  2:906570/151‑1 (MQ=255)                      gAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTa                                                                                                                                <  2:353951/149‑1 (MQ=255)                      gAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGc                                                                                                                              >  1:107876/1‑151 (MQ=255)                          gagTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACtg                                                                                                                          <  2:282766/151‑1 (MQ=255)                             tGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTgg                                                                                                                       >  1:119676/1‑151 (MQ=255)                                                    cgGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGa                                                                                                 <  2:310271/150‑1 (MQ=255)                                                      gTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAAc                                                                                               >  1:207679/1‑150 (MQ=255)                                                            cGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGggcgg                                                                                        <  1:294283/151‑1 (MQ=255)                                                              cTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGggc                                                                                          <  1:832554/147‑1 (MQ=255)                                                              cTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCt                                                                                      >  2:510861/1‑151 (MQ=255)                                                               tGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCt                                                                                      >  1:437458/1‑150 (MQ=255)                                                                    atGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCggg                                                                                >  2:473231/1‑151 (MQ=255)                                                                    atGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCggg                                                                                >  1:254657/1‑151 (MQ=255)                                                                      ggACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGgcg                                                                              <  1:454487/151‑1 (MQ=255)                                                                                gAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAg                                                                    <  1:861089/151‑1 (MQ=255)                                                                                     aCCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCaaa                                                               <  2:821155/151‑1 (MQ=255)                                                                                        aTTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGtt                                                            <  2:545875/151‑1 (MQ=255)                                                                                        aTTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGtt                                                            >  2:545080/1‑151 (MQ=255)                                                                                         ttGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGtt                                                            <  2:228598/150‑1 (MQ=255)                                                                                               gCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActctg                                                     <  1:556766/151‑1 (MQ=255)                                                                                                       gcgTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCg                                             >  2:866339/1‑151 (MQ=255)                                                                                                        cgTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCgg                                            >  2:290500/1‑151 (MQ=255)                                                                                                           aaaTTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGtt                                          <  2:463466/150‑1 (MQ=255)                                                                                                           aaaTTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGt                                           >  2:882455/1‑149 (MQ=255)                                                                                                                tGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTtattat                                     <  2:676625/150‑1 (MQ=255)                                                                                                                tGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTtattat                                     <  1:686158/150‑1 (MQ=255)                                                                                                                tGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTattatt                                    >  2:869294/1‑151 (MQ=255)                                                                                                                   cAACCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATtctc                                 >  2:158012/1‑151 (MQ=255)                                                                                                                    aaCCGTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCg                                <  2:295378/151‑1 (MQ=255)                                                                                                                        gTATTCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGaa                              >  2:246854/1‑149 (MQ=255)                                                                                                                            tCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTa                         <  2:875258/150‑1 (MQ=255)                                                                                                                            tCATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTAt                        >  1:723339/1‑151 (MQ=255)                                                                                                                             cATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTAtt                       <  1:281920/151‑1 (MQ=255)                                                                                                                             cATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTAtt                       >  1:32696/1‑151 (MQ=255)                                                                                                                             cATGGCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTAt                        <  2:598311/150‑1 (MQ=255)                                                                                                                                ggCCGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTc                     <  2:639071/150‑1 (MQ=255)                                                                                                                                  ccGTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCg                    >  1:785843/1‑149 (MQ=255)                                                                                                                                    gTTTCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCa                <  2:460275/151‑1 (MQ=255)                                                                                                                                       tCGGTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAgg              <  2:100468/150‑1 (MQ=255)                                                                                                                                         ggTGTTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGat           <  2:732935/151‑1 (MQ=255)                                                                                                                                          gtgtTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGata          <  2:514575/151‑1 (MQ=255)                                                                                                                                          gtgtTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGata          <  2:458244/151‑1 (MQ=255)                                                                                                                                           tgtTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGatat         >  2:33731/1‑151 (MQ=255)                                                                                                                                            gtTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATAtt        <  2:570256/151‑1 (MQ=255)                                                                                                                                            gtTGAAGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATAtt        >  2:175720/1‑151 (MQ=255)                                                                                                                                                 aGTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGc   >  1:875544/1‑151 (MQ=255)                                                                                                                                                  gTACAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCg  <  1:292129/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GAACATTATCGAGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCG  >  NC_000913/3093625‑3093920

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑