Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,082,863	C→T	T36T (ACC→ACT)	sbcB →	exodeoxyribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,082,863	0	C	T	100.0%	88.3 / NA	28	T36T (ACC→ACT)	sbcB	exodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/9);  total (19/9)

GCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCA  >  NC_000913/2082666‑2083095                                                                                                                                                                                                      |                                                                                                                                                                                                                                          gCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCg                                                                                                                                                                                       <  1:89186/250‑1 (MQ=255)       cGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCtgatga                                                                                                                                                                                <  1:490256/251‑1 (MQ=255)               aaTTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGAACAGTTCGCAGCCATTCGCACTGCTAGCGAATTCAATGTAATCGGCGAACCAGCGGACTTTTACTGCAAGCCCGCTGATGACAATTTAc                                                                                                                                                                        >  2:318002/1‑251 (MQ=255)                                                 tttGATACGCTGGCAGACAGCCGAAATAAAGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGCAGACTTTTATTGCAAGCCCGCTGATTACTATTTACCC‑CAGCCAGAAGCCGTAATAATTACCggatttac                                                                                                                                     >  2:123551/1‑245 (MQ=255)                                                                 aCAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACgg                                                                                                                     >  1:15995/1‑251 (MQ=255)                                                                           tAACGTATTTAAAATAATGCTAAATGACGGTAAGCCACAATCTACCTATTTGTTTCACGATTACGAAACCTTTGGCACGCACCACGCGTTAGATCTCCCGGCCCAGTTCGCAACCATTAGAACTGATAGCGAATTAAATGCAATAGACGAACCCGAACTCTTTTACTGAAATACCGCTGATGGCACTTTACCCTCA‑CCAGGACACTTATAAATTACCGGTATTGCCCGGCAGGaacaacagcccatacga                                                                                                            >  2:559060/1‑235 (MQ=255)                                                                            aaCGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAgg                                                                                                             >  1:339470/1‑248 (MQ=255)                                                                              cGGATTTAACCTAATGATGAATGACGTTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGaaa                                                                                                        <  1:100259/251‑1 (MQ=255)                                                                                                taatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTgccg                                                                                      >  1:537144/2‑251 (MQ=255)                                                                                                   tgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGccc                                                                                   >  1:228922/1‑251 (MQ=255)                                                                                                          aaGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGAACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCCAGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTCTTc                                                                             >  2:477820/1‑250 (MQ=255)                                                                                                                      aCCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtac                                                                >  1:332403/1‑251 (MQ=255)                                                                                                                              gTTTCACGATTACGAAACCGTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CCGCCAGGCGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGCAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtaccgtaccg                                                         >  2:71168/1‑250 (MQ=255)                                                                                                                                    cGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCt                                                  <  1:318002/251‑1 (MQ=255)                                                                                                                                           gAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCt                                           >  1:16414/1‑251 (MQ=255)                                                                                                                                                    ggccacccccccTCGTTAGATCGCCCTGCACAGTTCTCAGCCATTCGCACTGATAGTGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTacaa                                  <  2:156935/243‑1 (MQ=255)                                                                                                                                                    ggCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAGCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTacaa                                  <  1:372500/251‑1 (MQ=255)                                                                                                                                                         gcacCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CGGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAAtg                             >  1:396210/1‑251 (MQ=255)                                                                                                                                                         gcacCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAAGTCAATGGCAGCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAAtg                             <  1:349444/251‑1 (MQ=255)                                                                                                                                                          cacCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAAtgt                            >  1:433420/1‑251 (MQ=255)                                                                                                                                                               cgcgTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGtt                       >  1:426056/1‑251 (MQ=255)                                                                                                                                                                    tAGATCGCACTGCACAGTTAGCAGCCATTAGCACTGATAGCGAATTCAATGTCATAGGCGAACCCGAAGTCATTTAATGCAAGACCGATGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCCCAGGAAGCACGGTCGAAAGGCAAAAACGACGCAGCGTTTGCCGCTCGTATTCAATGGCTTTTTACCGTACCGAAGACCTGAATTCAGGGCTACAACAATATGCGTTTctcc                  >  2:556278/1‑248 (MQ=255)                                                                                                                                                                     aGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACg                                                                                                                      >  2:129564/1‑150 (MQ=255)                                                                                                                                                                     aGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACg                                                                                                                      <  1:129564/150‑1 (MQ=255)                                                                                                                                                                             cTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTACTTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTAACGTACCGAAGCCCTGTATTCGGGGCTACAACAATGTGCGTTTCGCCGACGAAGTc         >  2:143565/1‑251 (MQ=255)                                                                                                                                                                              tgctcaGTTCGCAGCCATTCGTTCTGATAGCGAATGCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCTAGGAGCCGTATTAATTACCGGTATTACCCCGAAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTCCCGAAGACCTGTATTCTGGGCTACAAAAATGTGCGTTTCGACGACGAAGTc         <  2:486787/246‑1 (MQ=255)                                                                                                                                                                                cacaGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACACTGTGCGTTTCGACGACGAAGTcac       >  2:450936/1‑250 (MQ=255)                                                                                                                                                                                    gTTCGCAGCCATTCGCACTGAAAGCGAATTAAATGTCATCGGCGAACCAGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACAA‑CAGCCAGGAGCCGTATTAATTACCGGTATTAACCCGCAAGAAGCAAGGGCGAAAGGAGAAAAAGAGACAGTGATTGCCGCCCGTATGCAGTCACTTTATACCGCTCCGAAGACTGGTTTACTAGGCTAAAACCATGTGCGtttaagataagaaatcccacgcc  >  2:368752/1‑231 (MQ=255)                                                                                                                                                                                                      |                                                                                                                                                                                                                                          GCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCC‑CAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCA  >  NC_000913/2082666‑2083095

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