Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	822,421	C→T	L26L (TTG→TTA)	ybhN ←	UPF0104 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	822,421	0	C	T	100.0%	143.9 / NA	41	L26L (TTG→TTA)	ybhN	UPF0104 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/21);  total (20/21)

GTAAAAGCGGGCGAGCAGGTCATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGA  >  NC_000913/822274‑822567                                                                                                                                                    |                                                                                                                                                   gTAAAAGCGGGCGAGCAGGTCATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTaa                                                                                                                                                  >  2:680770/1‑150 (MQ=255)                  ggTCATAGCAGCCGTAACTCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATaaa                                                                                                                                >  1:456658/1‑151 (MQ=255)                  ggTCATAGCAGCCGTAAATCAGATAGCTGGCGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATaaa                                                                                                                                >  1:626837/1‑151 (MQ=255)                  ggTCATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATaaa                                                                                                                                >  2:191817/1‑151 (MQ=255)                  ggTCATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATaaa                                                                                                                                >  2:431136/1‑151 (MQ=255)                  ggTCATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATaaa                                                                                                                                >  1:410185/1‑151 (MQ=255)                  ggTCATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATaaa                                                                                                                                >  1:529601/1‑151 (MQ=255)                    tCATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATaaaaa                                                                                                                              <  2:720005/151‑1 (MQ=255)                    tCATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATaaaaa                                                                                                                              <  2:440375/151‑1 (MQ=255)                      aTAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAAc                                                                                                                            <  2:266877/151‑1 (MQ=255)                                    tCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGAt                                                                                                              <  2:722536/151‑1 (MQ=255)                                    tCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGAt                                                                                                              <  1:17232/151‑1 (MQ=255)                                    tCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGAt                                                                                                              <  1:83513/151‑1 (MQ=255)                                       gATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGAtctt                                                                                                           >  2:102543/1‑151 (MQ=255)                                           gCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCttt                                                                                                       >  2:5706/1‑151 (MQ=255)                                           gCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCttt                                                                                                       >  2:212161/1‑151 (MQ=255)                                               acgacCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTa                                                                                                   <  1:124276/151‑1 (MQ=255)                                                  accaccAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAgc                                                                                                <  1:206103/151‑1 (MQ=255)                                                            cGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCgg                                                                                      >  1:273174/1‑151 (MQ=255)                                                                  cACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGa                                                                                >  2:582979/1‑151 (MQ=255)                                                                    cTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGAtt                                                                              <  2:108839/151‑1 (MQ=255)                                                                      aaGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTa                                                                            >  2:421837/1‑151 (MQ=255)                                                                         cagcGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTAc                                                                           <  2:404298/149‑1 (MQ=255)                                                                         cagcGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTAc                                                                           <  1:313735/149‑1 (MQ=255)                                                                         cagcGCAACGCGATTGCAGTCGAGGATGACCTTCCAGTCCTCTTCCCAGTCCACTATTTTGGAGAAGAGCACCAATAACACAGTCACCGCTATACAGAACAGCCAGGTGAGGATCTTCTTTGCCAAGCGCCAGCGCGGGTGTGATTTAc                                                                           <  1:36513/149‑1 (MQ=255)                                                                            cgcAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACtc                                                                         <  1:138234/148‑1 (MQ=255)                                                                            cgcAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACtc                                                                         <  2:534499/148‑1 (MQ=255)                                                                                          aGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTAcccc                                                         >  2:120919/1‑150 (MQ=255)                                                                                                     aCCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAg                                              <  2:376910/150‑1 (MQ=255)                                                                                                              aCCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTAcc                                      <  1:398675/149‑1 (MQ=255)                                                                                                              aCCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTAcc                                      <  1:365814/149‑1 (MQ=255)                                                                                                               ccTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGa                                   >  2:332927/1‑151 (MQ=255)                                                                                                                ctctTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGAt                                  <  2:266352/151‑1 (MQ=255)                                                                                                                 tctTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATc                                 <  1:354772/151‑1 (MQ=255)                                                                                                                  ctTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATcc                                <  1:680752/151‑1 (MQ=255)                                                                                                                  ctTCCCAGTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATc                                 >  2:385422/1‑150 (MQ=255)                                                                                                                         gTCCACTTTTTTGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATCCTGCGttt                         >  1:286558/1‑151 (MQ=255)                                                                                                                                    tGGCGTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATCCTGCGTTTCCATTGTCGGt              <  1:634361/151‑1 (MQ=255)                                                                                                                                        gTAGACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATCCTGCGTTTCCATTGTCGGTtgtg          >  1:724851/1‑151 (MQ=255)                                                                                                                                           gACCACCAGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATCCTGCGTTTCCATTGTCGGTTGTgcgg       >  2:94572/1‑151 (MQ=255)                                                                                                                                                ccaGTAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGa  >  2:626348/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GTAAAAGCGGGCGAGCAGGTCATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGA  >  NC_000913/822274‑822567

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