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%	91.3 / NA	28	T36T (ACC→ACT)	sbcB	exodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/10);  total (18/10)

GTGCAAGAAAATGTGACTACCCTCTCATTTTTATCTGACATGATCTGTTGCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTAT  >  NC_000913/2082617‑2083106                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    gTGCAAGAAAATGTGACTACCCTCTCATTTTTATCTGACATGATCTGTTGCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGAt                                                                                                                                                                                                                                                  >  1:125467/1‑250 (MQ=255)               gACTACCCTCTCATTTTTATCTGACATGATCTGTTGCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGt                                                                                                                                                                                                                                    <  2:132379/250‑1 (MQ=255)                 cTACCCTCTCATTTTTATCTGACATGATCTGTTGCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCAt                                                                                                                                                                                                                                 >  1:220154/1‑251 (MQ=255)                   aCCCTCTCATTTTTATCTGACATGATCTGTTGCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATc                                                                                                                                                                                                                                <  2:137293/250‑1 (MQ=255)                           ctttttATCTGCCATGATCTGTTGCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGGTCGCAGCCATTAGCACTGATCGCGAATTCAATGTCATCGGCTAAccc                                                                                                                                                                                                                       >  2:216401/2‑251 (MQ=255)                             ttttATCTGACATGATCTGTTGCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCg                                                                                                                                                                                                                      >  1:282169/1‑250 (MQ=255)                                            tCTGTTGCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGCGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCAAAGTCTTTTACTGCa                                                                                                                                                                                                       >  2:255550/1‑250 (MQ=255)                                                                                                 atTTGATCCTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAACCCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGCACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGAAtt                                                                                                                                                  >  2:192562/1‑250 (MQ=255)                                                                                                   ttGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGCACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTAcc                                                                                                                                               >  1:123086/1‑251 (MQ=255)                                                                                                     gATACTCTGGGAGACAGCAGAAATAACGGATTTAACCTAATGAGGAATGACGGTAATCAACAATCTACCTGTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTGAGAGCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGAGGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGTTATTAccc                                                                                                                                              <  1:197177/250‑1 (MQ=255)                                                                                                                                     tAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGa                                                                                                               >  1:63573/1‑249 (MQ=255)                                                                                                                                        ccTAATGATGAATGACGGTAAGCAACAATCTGCCTTTTTGTTGCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGGCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTGCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCcg                                                                                                          <  1:192562/251‑1 (MQ=255)                                                                                                                                                tgaatgaCGGTAAGCAACCATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAAccc                                                                                                                                                                                                                       >  2:142045/1‑134 (MQ=255)                                                                                                                                                tgaatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAAccc                                                                                                                                                                                                                       <  1:142045/134‑1 (MQ=255)                                                                                                                                                  aatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTAtt                                                                                                                                                  >  2:87013/1‑201 (MQ=255)                                                                                                                                                  aatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTAtt                                                                                                                                                  <  1:87013/201‑1 (MQ=255)                                                                                                                                                           aaGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCTCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGAGTTTGCCGCCCGTATTc                                                                                        <  2:8097/250‑1 (MQ=255)                                                                                                                                                              caacaaTCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTc                                                                                    >  2:109330/1‑251 (MQ=255)                                                                                                                                                                            tttGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTGTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtaccgta                                                                       <  2:84803/250‑1 (MQ=255)                                                                                                                                                                             ttGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtaccgtacc                                                                     >  1:88270/1‑251 (MQ=255)                                                                                                                                                                             tagTTTCAAGCTTACGAAACATTTGGGACGCACAACGAGTTAGATCGCCCTGTATAGGTCGCAGCGATTAGCCATGACAGCGAATTACATGTCCTCGGCGGACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTATTTAACGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGACGCCGCGTTTGACGCACGTGTGCACTCGCTTTTtaacttacc                                                                     >  2:140789/3‑244 (MQ=255)                                                                                                                                                                                     cGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCt                                                             >  2:176676/1‑251 (MQ=255)                                                                                                                                                                                      gATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTg                                                            >  1:273688/1‑251 (MQ=255)                                                                                                                                                                                                             ccccGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGc                                     <  2:220154/251‑1 (MQ=255)                                                                                                                                                                                                                                   caGTTCGCAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACTACGAAGTCACACg               <  2:46658/251‑1 (MQ=255)                                                                                                                                                                                                                                          cAGCCATTCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAtt        >  1:248199/1‑251 (MQ=255)                                                                                                                                                                                                                                           aGCCATTCGCACTGATAGCGAATTCAATGTCATCAGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTAAAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAttt       >  1:185111/1‑251 (MQ=255)                                                                                                                                                                                                                                                ttCGCACTGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTAt  >  1:247177/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    GTGCAAGAAAATGTGACTACCCTCTCATTTTTATCTGACATGATCTGTTGCCACTCGCTGCCAAATTGTGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTAT  >  NC_000913/2082617‑2083106

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