Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	964,086	C→T	R87C (CGC→TGC)	ihfB →	integration host factor (IHF), DNA‑binding protein, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	964,086	0	C	T	100.0%	152.9 / NA	44	R87C (CGC→TGC)	ihfB	integration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/23);  total (21/23)

TCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAA  >  NC_000913/963961‑964230                                                                                                                              |                                                                                                                                                 tCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTgtgc                                                                                                                                                >  4:21404/1‑128 (MQ=255) tCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTgtgc                                                                                                                                                <  3:21404/128‑1 (MQ=255) tCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGt                                                                                                                           >  5:97820/1‑149 (MQ=255)             gCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTAc                                                                                                               <  6:97820/149‑1 (MQ=255)              cAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATcgc                                                                                                                                          <  5:132489/121‑1 (MQ=255)              cAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATcgc                                                                                                                                          >  6:132489/1‑121 (MQ=255)              cAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACt                                                                                                              >  5:27510/1‑149 (MQ=255)                gTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTa                                                                                                                               >  8:107904/1‑130 (MQ=255)                gTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTa                                                                                                                               <  7:107904/130‑1 (MQ=255)                   tctcTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCaaa                                                                                                          <  6:27510/148‑1 (MQ=255)                    ctctTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAAc                                                                                                         <  3:19749/148‑1 (MQ=255)                              tACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGtttt                                                                                                                   <  6:143311/128‑1 (MQ=255)                              tACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGtttt                                                                                                                   >  5:143311/1‑128 (MQ=255)                              tACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACga                                                                                                 >  2:163776/1‑146 (MQ=255)                                ccGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgagaga                                                                                             >  3:2227/1‑148 (MQ=255)                                ccGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGaa                                                                                            <  4:2227/149‑1 (MQ=255)                                           gTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAg                                                                                                                       >  1:72067/1‑111 (MQ=255)                                           gTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAg                                                                                                                       <  2:72067/111‑1 (MQ=255)                                           gTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTc                                                                                 >  6:148378/1‑149 (MQ=255)                                              ccGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                               <  5:148378/148‑1 (MQ=255)                                                 gACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGc                                                                           >  4:73292/1‑149 (MQ=255)                                                      aaTCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTc                                                                                                             <  8:142008/110‑1 (MQ=255)                                                      aaTCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTc                                                                                                             >  7:142008/1‑110 (MQ=255)                                                            aaGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGATAAGCACCTGTCGGGTGCTTTTTTCAttt                                                                <  1:163776/149‑1 (MQ=255)                                                              gACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                               >  1:161976/1‑132 (MQ=255)                                                              gACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                               <  2:161976/132‑1 (MQ=255)                                                                          aGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                               <  4:24501/120‑1 (MQ=255)                                                                          aGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                               >  3:24501/1‑120 (MQ=255)                                                                           gTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaact                                                 <  1:44557/149‑1 (MQ=255)                                                                               aaCTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactggaa                                             <  3:73292/149‑1 (MQ=255)                                                                                       ggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTc                                     <  4:4132/149‑1 (MQ=255)                                                                                        gaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctg                                                      >  2:70880/1‑131 (MQ=255)                                                                                        gaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctg                                                      <  1:70880/131‑1 (MQ=255)                                                                                              aCGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGt                              >  6:161424/1‑149 (MQ=255)                                                                                                    cTCACTTTAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGaa                        <  5:161424/149‑1 (MQ=255)                                                                                                    cTCACTTTAAACCTGGTAAAGAACTGTGCGATAGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgaga                                                                                               >  1:161494/1‑78 (MQ=255)                                                                                                         tttAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctg                                                      <  5:5341/114‑1 (MQ=255)                                                                                                         tttAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctg                                                      >  6:5341/1‑114 (MQ=255)                                                                                                         tttAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctg                                                      >  3:57511/1‑114 (MQ=255)                                                                                                         tttAAACCTGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctg                                                      <  4:57511/114‑1 (MQ=255)                                                                                                                 tGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACga                                                                                                 >  1:5883/1‑63 (MQ=255)                                                                                                                 tGGTAAAGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACga                                                                                                 <  2:5883/63‑1 (MQ=255)                                                                                                                     aaaGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCg       >  2:33441/1‑149 (MQ=255)                                                                                                                      aaGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaa                                                   <  2:34172/104‑1 (MQ=255)                                                                                                                      aaGAACTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaa                                                   >  1:34172/1‑104 (MQ=255)                                                                                                                          aCTGTGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTaaa  <  1:33441/149‑1 (MQ=255)                                                                                                                              |                                                                                                                                                 TCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAA  >  NC_000913/963961‑964230

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