Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,036,373	C→T	*578* (TAG→TAA)	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,036,373	0	C	T	100.0%	219.8 / NA	68	*578* (TAG→TAA)	recJ	ssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/31);  total (37/31)

GTGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTAAAAG  >  NC_000913/3036239‑3036506                                                                                                                                       |                                                                                                                                      gTGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                                   <  2:1489205/139‑1 (MQ=255) gTGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                                   <  2:394673/139‑1 (MQ=255)     ggTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                               <  2:811448/139‑1 (MQ=255)       tCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGcc                                                                                                                             >  2:1311168/1‑139 (MQ=255)                   ttattTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCga                                                                                                                 >  2:1239141/1‑139 (MQ=255)                         aCCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATaa                                                                                                           >  1:1438066/1‑139 (MQ=255)                                               tCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcgcg                                                                                     >  1:1035994/1‑139 (MQ=255)                                                   aTTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCa                                                                                                     <  1:246665/119‑1 (MQ=255)                                                   aTTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCa                                                                                                     >  2:246665/1‑119 (MQ=255)                                                         tttATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                               <  2:640989/87‑1 (MQ=255)                                                         tttATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                               >  1:640989/1‑87 (MQ=255)                                                         tttATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                               <  1:249992/87‑1 (MQ=255)                                                         tttATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                               >  2:249992/1‑87 (MQ=255)                                                             ttGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATTCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGatat                                                                       >  1:1417413/1‑139 (MQ=255)                                                              tGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATc                                                                      <  1:1239141/139‑1 (MQ=255)                                                                aGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCga                                                                                                                 <  2:983973/94‑1 (MQ=255)                                                                aGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCga                                                                                                                 >  1:983973/1‑94 (MQ=255)                                                                aGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                    >  1:533192/1‑139 (MQ=255)                                                                          aaTGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGc                                                          <  1:1311168/139‑1 (MQ=255)                                                                          aaTGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGc                                                          <  2:1417413/139‑1 (MQ=255)                                                                               ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                                   <  2:1054711/61‑1 (MQ=255)                                                                               ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                               <  2:248627/65‑1 (MQ=255)                                                                               ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                               >  1:248627/1‑65 (MQ=255)                                                                               ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                                  <  2:804196/110‑1 (MQ=255)                                                                               ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                                  >  1:804196/1‑110 (MQ=255)                                                                               ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTGTTTTATAGAGGAGATGACGTTaaa                                                                                                                                   >  1:1054711/1‑61 (MQ=255)                                                                                    tAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACt                                                >  1:272926/1‑139 (MQ=255)                                                                                           gggATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                         <  2:533192/139‑1 (MQ=255)                                                                                             gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCgcgc                                          >  1:137349/1‑136 (MQ=255)                                                                                             gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCgcgc                                          <  2:137349/136‑1 (MQ=255)                                                                                                  tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                                  <  1:1358743/91‑1 (MQ=255)                                                                                                  tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                                  >  2:1358743/1‑91 (MQ=255)                                                                                                  tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGtt                                  >  1:611497/1‑139 (MQ=255)                                                                                                       aaTCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCgatg                                                                                                               <  2:567019/57‑1 (MQ=255)                                                                                                       aaTCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCgatg                                                                                                               >  1:567019/1‑57 (MQ=255)                                                                                                         tCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCga                                                                                                                 >  2:935808/1‑53 (MQ=255)                                                                                                         tCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCga                                                                                                                 <  1:935808/53‑1 (MQ=255)                                                                                                         tCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGc                           <  2:1551063/139‑1 (MQ=255)                                                                                                            aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGt                                                                                                                    <  2:540718/47‑1 (MQ=255)                                                                                                            aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGt                                                                                                                    >  1:540718/1‑47 (MQ=255)                                                                                                            aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCa                                                                                                     <  2:960226/62‑1 (MQ=255)                                                                                                            aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCa                                                                                                     >  1:960226/1‑62 (MQ=255)                                                                                                                     tttATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTa               >  2:134978/1‑139 (MQ=255)                                                                                                                          agagAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAggg                    <  1:845831/129‑1 (MQ=255)                                                                                                                          agagAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAggg                    >  2:845831/1‑129 (MQ=255)                                                                                                                            agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGc                                                                                                  >  1:285429/1‑49 (MQ=255)                                                                                                                            agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGc                                                                                                  <  2:285429/49‑1 (MQ=255)                                                                                                                            agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGata                                                                        >  2:573068/1‑75 (MQ=255)                                                                                                                            agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGata                                                                        <  1:573068/75‑1 (MQ=255)                                                                                                                            agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATaa                                                            <  2:260733/87‑1 (MQ=255)                                                                                                                            agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATaa                                                            >  1:260733/1‑87 (MQ=255)                                                                                                                            agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACAt        >  2:667271/1‑139 (MQ=255)                                                                                                                            agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACAt        >  1:40350/1‑139 (MQ=255)                                                                                                                            agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACAt        >  2:879740/1‑139 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                                  <  1:1173346/63‑1 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                    >  1:1418031/1‑77 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                    <  2:1418031/77‑1 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAg                                                                   >  1:212750/1‑78 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAg                                                                   <  2:212750/78‑1 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATc                               >  2:417500/1‑114 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATc                               <  1:417500/114‑1 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa      >  1:1500548/1‑139 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa      >  1:480434/1‑139 (MQ=255)                                                                                                                              aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCACGCTGCGGTTGCCGCGaaa                                                                                  >  2:1173346/1‑63 (MQ=255)                                                                                                                                gaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTaaa    >  2:51481/1‑139 (MQ=255)                                                                                                                                 aTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCg            <  2:1100679/130‑1 (MQ=255)                                                                                                                                 aTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCg            >  1:1100679/1‑130 (MQ=255)                                                                                                                                  tGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTAAAAg  <  2:480434/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      GTGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTAAAAG  >  NC_000913/3036239‑3036506

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