Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,900,610	C→T	W537* (TGG→TGA)	bglH ←	carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,900,610	0	C	T	100.0%	92.2 / NA	28	W537* (TGG→TGA)	bglH	carbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/10);  total (18/10)

GTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCC  >  NC_000913/3900365‑3900841                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        gTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGAt                                                                                                                                                                                                                                    <  1:56501/251‑1 (MQ=255)   ccggaaccggaaCACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGTCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGAtt                                                                                                                                                                                                                                   <  2:229281/250‑1 (MQ=255)          cggaacACCCACCACTACAGAGACATTTTTTGCACCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTAACCGGGACAGTGGGTGCCGGGCTGGCTGGCATATCAGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCCATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTg                                                                                                                                                                                                                           >  1:292251/1‑251 (MQ=255)           ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTgg                                                                                                                                                                                                                          >  1:106451/1‑251 (MQ=255)                 accaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACAAACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCAGATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCCGCCTGGGCACCg                                                                                                                                                                                                                    >  2:110514/2‑251 (MQ=255)                                                         aGCGAAAGGTCCCGCTGTTATCGGCATTGCCCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGACATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGc                                                                                                                                                                            >  2:125433/1‑251 (MQ=255)                                                                 gTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGt                                                                                                                                                                    <  2:49435/251‑1 (MQ=255)                                                                                 cATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCtt                                                                                                                                                    >  1:84923/1‑251 (MQ=255)                                                                                       acTGAGTCACATACTGTTTCACAGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGatat                                                                                                                                              >  2:195306/2‑251 (MQ=255)                                                                                                   aCTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGCAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGt                                                                                                                                     >  2:132469/1‑248 (MQ=255)                                                                                                              gggATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTTTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTcc                                                                                                                       <  2:149636/251‑1 (MQ=255)                                                                                                              gggATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTcc                                                                                                                       >  1:292822/1‑251 (MQ=255)                                                                                                                   aGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGAcg                                                                                                                  >  1:132045/1‑251 (MQ=255)                                                                                                                          gacGGGCTGGCTAGCATATACGCTGCTATGGCACACTGAGGAGAAATACCGCTGGCAATAACCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGCGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGAGTAGGAACGAATTTCCGGActcgctgtc                                                                                                           >  2:77066/3‑243 (MQ=255)                                                                                                                                            tCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACt                                                                                         >  1:171800/1‑251 (MQ=255)                                                                                                                                               gctgctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACttt                                                                                       >  2:43557/1‑250 (MQ=255)                                                                                                                                                  gctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAg                                                                                   <  1:153577/251‑1 (MQ=255)                                                                                                                                                     aTGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGt                                                                                 <  1:132469/250‑1 (MQ=255)                                                                                                                                                                    gAAATACCGCTGACAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCAACAGCAAACTGGCCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGAATGGAAGAGCGTGGtt                                                                 >  2:285777/1‑251 (MQ=255)                                                                                                                                                                      aatgccgttGGCAATAGCCAGTGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTGCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTGATTTTCTGCAAAGGTGAAGCTATCCTGTTCGTTGTCCAGGGCTTTGATATACGTGGCGTAGAAGCGAATTTCCTGACTCGAGGACAACGGAGGGGATTTGAGTTTAAAGGTATGGAAGAGCGTGGtttt                                                               <  1:22699/243‑1 (MQ=255)                                                                                                                                                                               ggCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAAcc                                                         >  2:105937/1‑248 (MQ=255)                                                                                                                                                                                           cGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAAttt                                          >  2:262252/1‑251 (MQ=255)                                                                                                                                                                                                  gcagcaATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACt                                    >  2:105994/1‑250 (MQ=255)                                                                                                                                                                                                  gcagcaATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTa                                   <  2:304848/251‑1 (MQ=255)                                                                                                                                                                                                    agcaATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTAtt                                 <  1:103666/251‑1 (MQ=255)                                                                                                                                                                                                     gcaATTTTTATATTAATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTc                                >  1:206498/1‑251 (MQ=255)                                                                                                                                                                                                        aTTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCa                             <  2:84923/251‑1 (MQ=255)                                                                                                                                                                                                                                    cATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGccc  >  1:282152/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        GTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCC  >  NC_000913/3900365‑3900841

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