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%	171.8 / NA	54	*578* (TAG→TAA)	recJ	ssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/28);  total (26/28)

CGTGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTA  >  NC_000913/3036238‑3036502                                                                                                                                        |                                                                                                                                  cGTGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaa                                                                                                                                <  2:583386/139‑1 (MQ=255)       gTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGc                                                                                                                          <  1:645669/139‑1 (MQ=255)         ccTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCa                                                                                                                        >  1:1058345/1‑139 (MQ=255)             aaTGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGc                                                                                                                          <  1:1226081/133‑1 (MQ=255)             aaTGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGc                                                                                                                          >  2:1226081/1‑133 (MQ=255)                                  aaTTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAgg                                                                                               <  1:923908/139‑1 (MQ=255)                                              ggTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcg                                                                                   >  1:761346/1‑139 (MQ=255)                                              ggTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcg                                                                                   >  2:1118615/1‑139 (MQ=255)                                                 cTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGa                                                                                <  2:1058345/139‑1 (MQ=255)                                                    aTTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCgatgat                                                                                                         >  2:642972/1‑111 (MQ=255)                                                    aTTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCgatgat                                                                                                         <  1:642972/111‑1 (MQ=255)                                                                   cTAGTCAAAATGCGGTGATGAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGc                                                              >  1:897021/1‑139 (MQ=255)                                                                   cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGtt                                                                        <  2:1045818/129‑1 (MQ=255)                                                                   cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGtt                                                                        >  1:1045818/1‑129 (MQ=255)                                                                           aaTGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGt                                                                                        >  1:642275/1‑105 (MQ=255)                                                                           aaTGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGt                                                                                        <  2:642275/105‑1 (MQ=255)                                                                                ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGc                                                                                                                          >  1:678761/1‑66 (MQ=255)                                                                                ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGc                                                                                                                          <  2:678761/66‑1 (MQ=255)                                                                                ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGtt                                                 >  1:622928/1‑139 (MQ=255)                                                                                 gTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                               <  2:1282593/60‑1 (MQ=255)                                                                                 gTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                               >  1:1282593/1‑60 (MQ=255)                                                                                 gTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                               <  2:79048/60‑1 (MQ=255)                                                                                 gTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                               >  1:79048/1‑60 (MQ=255)                                                                                            gggATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     <  2:622928/139‑1 (MQ=255)                                                                                            gggATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     <  2:819882/139‑1 (MQ=255)                                                                                              gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGcc                                                                                                                         >  2:1053824/1‑53 (MQ=255)                                                                                              gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGcc                                                                                                                         <  1:1053824/53‑1 (MQ=255)                                                                                              gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGata                                                                                                                    <  1:347349/58‑1 (MQ=255)                                                                                              gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGata                                                                                                                    >  2:347349/1‑58 (MQ=255)                                                                                                ttGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     <  1:352931/135‑1 (MQ=255)                                                                                                ttGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     >  2:352931/1‑135 (MQ=255)                                                                                                   tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATaa                                                        <  2:1134539/113‑1 (MQ=255)                                                                                                   tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATaa                                                        >  1:1134539/1‑113 (MQ=255)                                                                                                      ccAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTGGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATc                           <  2:897021/139‑1 (MQ=255)                                                                                                             aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                              >  1:1216811/1‑81 (MQ=255)                                                                                                             aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                              <  2:1216811/81‑1 (MQ=255)                                                                                                             aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGc                                                      <  1:56850/105‑1 (MQ=255)                                                                                                             aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGc                                                      >  2:56850/1‑105 (MQ=255)                                                                                                             aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCgcgc                                      >  2:1240666/1‑121 (MQ=255)                                                                                                             aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCgcgc                                      <  1:1240666/121‑1 (MQ=255)                                                                                                             aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCag                    >  2:590837/1‑139 (MQ=255)                                                                                                              cGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGGGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCaga                   <  1:590837/139‑1 (MQ=255)                                                                                                                       ttATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTAt          <  2:761346/139‑1 (MQ=255)                                                                                                                         ataGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGtt                                                                                       <  1:974499/60‑1 (MQ=255)                                                                                                                         ataGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGtt                                                                                       >  2:974499/1‑60 (MQ=255)                                                                                                                             agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     >  2:46647/1‑106 (MQ=255)                                                                                                                             agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     <  1:46647/106‑1 (MQ=255)                                                                                                                             agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACAt    >  2:163741/1‑139 (MQ=255)                                                                                                                              gaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACAtt   <  1:163741/139‑1 (MQ=255)                                                                                                                               aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGa                                                                                >  2:774052/1‑61 (MQ=255)                                                                                                                               aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGa                                                                                <  1:774052/61‑1 (MQ=255)                                                                                                                               aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa  >  1:982119/1‑139 (MQ=255)                                                                                                                                  aTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCg                                                                          >  1:325241/1‑64 (MQ=255)                                                                                                                                  aTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCg                                                                          <  2:325241/64‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  CGTGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTA  >  NC_000913/3036238‑3036502

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: ‑