Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,037,738	A→T	D123E (GAT→GAA)	recJ ←	ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,037,738	0	A	T	100.0%	167.3 / NA	42	D123E (GAT→GAA)	recJ	ssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (21/21);  total (21/21)

CCGCTGCGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCG  >  NC_000913/3037593‑3037876                                                                                                                                                  |                                                                                                                                           ccGCTGCGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGAc                                                                                                                                       <  2:355092/151‑1 (MQ=255)   gCTGCGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGAcc                                                                                                                                      >  2:299708/1‑150 (MQ=255)   gCTGCGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCa                                                                                                                                     <  1:124347/151‑1 (MQ=255)           tAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCggg                                                                                                                             >  1:679135/1‑151 (MQ=255)                          cAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCg                                                                                                              >  2:557196/1‑151 (MQ=255)                            aaTGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTc                                                                                                            >  2:36374/1‑151 (MQ=255)                               ggTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCtt                                                                                                          <  1:272034/150‑1 (MQ=255)                                  gATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGa                                                                                                       <  2:79586/150‑1 (MQ=255)                                   aTCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGaa                                                                                                      >  2:157940/1‑150 (MQ=255)                                      ggTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACgg                                                                                                  >  1:512208/1‑151 (MQ=255)                                          aCAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTg                                                                                              <  2:424742/151‑1 (MQ=255)                                           cAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTgg                                                                                             >  2:330224/1‑151 (MQ=255)                                             aTAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTa                                                                                           <  1:523022/151‑1 (MQ=255)                                             aTAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTa                                                                                           <  1:144760/151‑1 (MQ=255)                                               aaCCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTAcc                                                                                         >  1:465949/1‑151 (MQ=255)                                                   gggATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGACTTCGAAACGGTTTGGTACCAGGt                                                                                     >  2:449500/1‑151 (MQ=255)                                                      aTGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGt                                                                                  <  1:68220/151‑1 (MQ=255)                                                      aTGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGt                                                                                  <  1:542129/151‑1 (MQ=255)                                                        gCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCg                                                                                >  2:666320/1‑151 (MQ=255)                                                                 gcgagcgTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATtgctg                                                                       >  1:480331/1‑151 (MQ=255)                                                                   gagcgTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATtgctgc                                                                      >  1:105192/1‑150 (MQ=255)                                                                                   cGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATc                                                     <  2:117095/151‑1 (MQ=255)                                                                                         ggaggaAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGc                                               <  1:343569/151‑1 (MQ=255)                                                                                              aaaTACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCt                                          >  2:417115/1‑151 (MQ=255)                                                                                                    cGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGc                                    <  2:347718/151‑1 (MQ=255)                                                                                                        aTCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTg                                <  1:185976/151‑1 (MQ=255)                                                                                                              cGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCGCGCTTAGAGCCGTGCTGGTc                          <  2:404636/151‑1 (MQ=255)                                                                                                               gTGACAATTAACTGCGCGCCACGGGCATGGGCCTGTCCGACCACTTCCGCGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCCGGTAGGCGATATTGCTGCAACCAACCGACCGCATCGCCACCACCCTCAGAGCCGTGCTGGTc                          >  2:518087/1‑150 (MQ=255)                                                                                                                 gACAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTcgc                        <  1:266369/150‑1 (MQ=255)                                                                                                                  aCAATTAACTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCgcgc                      >  2:272391/1‑151 (MQ=255)                                                                                                                          cTGCGCGCCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCg              >  2:340988/1‑151 (MQ=255)                                                                                                                            gcgcgcCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTc            >  2:459552/1‑151 (MQ=255)                                                                                                                            gcgcgcCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTc            <  1:124249/151‑1 (MQ=255)                                                                                                                             cgcgcCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCg           <  1:332858/151‑1 (MQ=255)                                                                                                                              gcgcCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGa          <  2:200844/151‑1 (MQ=255)                                                                                                                                gcCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGaaa        <  2:244333/151‑1 (MQ=255)                                                                                                                                gcCACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGaa         >  1:624952/1‑150 (MQ=255)                                                                                                                                 ccACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAAt       >  1:144297/1‑151 (MQ=255)                                                                                                                                  cACGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGCCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATc      >  1:395137/1‑151 (MQ=255)                                                                                                                                   aCGGGCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCa     <  2:406967/151‑1 (MQ=255)                                                                                                                                     gggCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCAcc   <  2:168411/151‑1 (MQ=255)                                                                                                                                       gCATGGGCCTGTTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCg  >  2:287683/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                           CCGCTGCGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCG  >  NC_000913/3037593‑3037876

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