Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,385,869	(A)7→6	coding (10/972 nt)	yfbL →	putative M28A family peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,385,863	0	A	.	100.0%	130.6 / NA	32	coding (4/972 nt)	yfbL	putative M28A family peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (18/14);  total (18/14)

TAAATAAACTCAGGCACCAGACATAAAATATACTCTTTGTATAATCCTGCTAATCTTAATAAAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATGAAAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGCAAGTCCTGAGCAATTAGAAAAAACAGTTCG  >  NC_000913/2385735‑2385996                                                                                                                                 |                                                                                                                                      taaataaaCTCAGGCACCAGACATAAAATATACTCTTTGTATAATCCTGCTAATCTTAATAAAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATc                                                                                                                            <  1:88701/139‑1 (MQ=255)                                 ctctTTGTATAATCCTGCTAATCTTAATAAAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACt                                                                                            >  1:991401/1‑139 (MQ=255)                                      tGTATAATCCTGCTAATCTTAATAAAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTttatta                                                                                                               >  1:1411746/1‑115 (MQ=255)                                      tGTATAATCCTGCTAATCTTAATAAAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTttatta                                                                                                               <  2:1411746/115‑1 (MQ=255)                                            aTCCTGCTAATCTTAATAAAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGatta                                                                                 >  2:325221/1‑139 (MQ=255)                                                tGCTAATCTTAATAAAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTg                                                                                                                        <  1:1290845/96‑1 (MQ=255)                                                tGCTAATCTTAATAAAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTg                                                                                                                        >  2:1290845/1‑96 (MQ=255)                                                           ataaAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCt                                                                  >  1:525534/1‑139 (MQ=255)                                                                      tatCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTa                                                                                              <  1:63010/100‑1 (MQ=255)                                                                      tatCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTa                                                                                              >  2:63010/1‑100 (MQ=255)                                                                          gAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCtttt                                                                                                                   >  1:447975/1‑75 (MQ=255)                                                                          gAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCtttt                                                                                                                   <  2:447975/75‑1 (MQ=255)                                                                           aaaCCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTAcc                                                  >  1:820233/1‑139 (MQ=255)                                                                             aCCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTAc                                                                         >  1:697103/1‑114 (MQ=255)                                                                             aCCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTAc                                                                         <  2:697103/114‑1 (MQ=255)                                                                             aCCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGt                                                              <  2:1020155/125‑1 (MQ=255)                                                                             aCCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGt                                                              >  1:1020155/1‑125 (MQ=255)                                                                                  tAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTttat                                                                                                                 <  1:1413068/69‑1 (MQ=255)                                                                                  tAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTttat                                                                                                                 >  2:1413068/1‑69 (MQ=255)                                                                                             gATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGc                                <  2:525534/139‑1 (MQ=255)                                                                                             gATTCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGc                                >  2:1282979/1‑139 (MQ=255)                                                                                                tCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGa                                     >  1:787752/1‑131 (MQ=255)                                                                                                tCAACAACGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGa                                     <  2:787752/131‑1 (MQ=255)                                                                                                  aacaacGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCa                                                            >  2:20504/1‑106 (MQ=255)                                                                                                  aacaacGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCa                                                            <  1:20504/106‑1 (MQ=255)                                                                                                       cGCAACCCATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGCAAGTCCTGAg                      >  2:430843/1‑139 (MQ=255)                                                                                                            cccATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGc                                >  2:331547/1‑124 (MQ=255)                                                                                                            cccATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGc                                <  1:331547/124‑1 (MQ=255)                                                                                                            cccATTAATAATTTGCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGCAAGTCCTGAGCAAtt                 >  2:527454/1‑139 (MQ=255)                                                                                                                          gCATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGCAAGTCCTGAGCAATTAGAAAAAACAGTTc   >  2:806650/1‑139 (MQ=255)                                                                                                                           cATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGCAAGTCCTGAGCAATTAGAAAAAACAGTTCg  <  1:325221/139‑1 (MQ=255)                                                                                                                           cATATG‑AAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGCAAGTCCTGAGCAATTAGAAAAAACAGTTCg  <  2:991401/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                      TAAATAAACTCAGGCACCAGACATAAAATATACTCTTTGTATAATCCTGCTAATCTTAATAAAATACAATATCGAAACCCATAGCAAATAATGATTCAACAACGCAACCCATTAATAATTTGCATATGAAAAAAATAATCTTTGCTTTTATTATATTATTTGTGTTTTTACTCCCTATGATTATTTTTTACCAACCCTGGGTCAATGCATTACCGTCAACGCCCCGACACGCAAGTCCTGAGCAATTAGAAAAAACAGTTCG  >  NC_000913/2385735‑2385996

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