Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,577,059	A→G	T187T (ACT→ACC)	yddB ←	putative TonB‑dependent receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,577,059	0	A	G	100.0%	120.4 / NA	34	T187T (ACT→ACC)	yddB	putative TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/8);  total (26/8)

TTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCT  >  NC_000913/1576927‑1577193                                                                                                                                     |                                                                                                                                       ttCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATGCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGAt                                                                                                                                 <  2:668258/140‑1 (MQ=255)   ccTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATaa                                                                                                                               <  1:112166/140‑1 (MQ=255)           aaCGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAAttt                                                                                                                       >  2:61101/1‑140 (MQ=255)           aaCGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAAttt                                                                                                                       >  1:52880/1‑140 (MQ=255)                 aaGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCtt                                                                                                                 >  1:77096/1‑140 (MQ=255)                   ggTATAAAAGTTCTTTTTAAAATATGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGc                                                                                                               >  1:537192/1‑140 (MQ=255)                               cTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGt                                                                                                   >  1:75857/1‑140 (MQ=255)                               cTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGCCCGAACGCCTGGCGCGTTAACCCAATTTCACCTTGCTATCATCAGCGt                                                                                                   >  2:275546/1‑140 (MQ=255)                                                  taataaGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGc                                                                                <  1:550101/140‑1 (MQ=255)                                                  taataaGTACTTCCTGAAGAGCCTTCATTGAATGCGCACTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGc                                                                                <  1:275546/140‑1 (MQ=255)                                                   aataaGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGAt                                                                                      <  1:33038/133‑1 (MQ=255)                                                   aataaGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGAt                                                                                      >  2:33038/1‑133 (MQ=255)                                                            ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcc                                                                      >  1:676841/1‑140 (MQ=255)                                                              ccTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATTATCAGCGATGAAGCGTATGATATTTGCATCAATTAcccc                                                                    >  2:345822/1‑140 (MQ=255)                                                                   agagCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCAt                                                               >  1:289542/1‑140 (MQ=255)                                                                     agCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTg                                                             >  1:561604/1‑140 (MQ=255)                                                                      gCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGa                                                            >  1:443932/1‑140 (MQ=255)                                                                       ccTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa                                                           >  2:23138/1‑140 (MQ=255)                                                                                 aTGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTTATCGTTGCATCAATTAGCCCGGCATTGAAGCGACCAAAt                                                 >  2:201333/1‑140 (MQ=255)                                                                                      ctctTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc                                            >  1:543801/1‑140 (MQ=255)                                                                                           gTTATTCTCATCGATCTGCGATTTTGACCAGTCCGAAGGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGca                                       <  1:302736/140‑1 (MQ=255)                                                                                           gTTATTCTCATCGATCTGCGATCTTGACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGca                                       <  1:345822/140‑1 (MQ=255)                                                                                               ttCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaa                                   >  2:466498/1‑140 (MQ=255)                                                                                                tctcATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa                                  <  2:75857/140‑1 (MQ=255)                                                                                                   cATCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTg                               >  2:532879/1‑140 (MQ=255)                                                                                                    aTCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAGGCGCCCAAATTCAACCGGCACAAAACTGt                              >  1:334404/1‑140 (MQ=255)                                                                                                    aTCGATATGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGt                              >  1:709392/1‑140 (MQ=255)                                                                                                         tatGCGATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATaa                         >  2:62362/1‑140 (MQ=255)                                                                                                             cGATGTTACCCAGTCCGAACGGTGGGTGCGATAACGCAAGTTCCCCTTCCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCCACCAAATTCAACCGGCACAAAACTGTCATAAAGCg                     <  2:289542/140‑1 (MQ=255)                                                                                                              gATGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGt                    >  1:25881/1‑140 (MQ=255)                                                                                                               aTGTTAACCAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCAACCAAATTCAAGCGGCACAAACCTGTCATAAAGCGTc                   >  1:253824/1‑140 (MQ=255)                                                                                                                  ttAACCAGTCCGAGCGGGTGTGGCAATACCCCAATTTCACCTTGCTATCATGAGCGTTGAAGCTTTAGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca                <  2:443932/140‑1 (MQ=255)                                                                                                                       cAGTCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTc           >  1:418717/1‑140 (MQ=255)                                                                                                                          tCCGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg        >  2:123203/1‑140 (MQ=255)                                                                                                                            cGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTa      >  2:558776/1‑140 (MQ=255)                                                                                                                            cGAACGCGTGGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTa      >  2:735287/1‑140 (MQ=255)                                                                                                                                cgcgTGGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt  <  2:77096/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       TTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCT  >  NC_000913/1576927‑1577193

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