Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	131,443	Δ1 bp	intergenic (‑183/‑172)	yacH ← / → acnB	DUF3300 domain‑containing protein YacH/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	131,442	0	A	.	100.0%	129.0 / NA	31	intergenic (‑182/‑173)	yacH/acnB	DUF3300 domain‑containing protein YacH/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (14/17);  total (14/17)

GGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCC  >  NC_000913/131363‑131516                                                                                |                                                                           ggATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTac                                                                           <  1:31874/81‑2 (MQ=255)  gATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTacg                                                                          <  2:31875/81‑3 (MQ=255)       gTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTc                                                                     <  1:31876/80‑1 (MQ=255)       gTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTc                                                                     <  2:31877/80‑1 (MQ=255)         tGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGc                                                                  >  1:31886/1‑81 (MQ=255)            gaAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCAc                                                               <  1:31878/81‑1 (MQ=255)             aaaGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACt                                                              >  1:31909/1‑81 (MQ=255)             aaaGGTGATTATCACCATGCGAATTAACGAAATTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACt                                                              <  2:31879/81‑1 (MQ=255)                  tGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACtttttt                                                         <  1:31880/81‑1 (MQ=255)                  tGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACttttt                                                          >  2:31915/1‑80 (MQ=255)                    aTTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTAt                                                       <  1:31881/81‑1 (MQ=255)                        tCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTg                                                   <  2:31882/81‑1 (MQ=255)                           ccATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCtt                                                <  2:31883/81‑1 (MQ=255)                                 gAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGt                                          <  2:31884/81‑1 (MQ=255)                                  aaTTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGt                                          >  1:31920/1‑80 (MQ=255)                                     tAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTaaa                                       >  2:31916/1‑80 (MQ=255)                                          aaGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGAtt                                 <  1:31885/81‑1 (MQ=255)                                                   cGGAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTc                         <  2:31886/80‑1 (MQ=255)                                                    ggAGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGt                       <  1:31887/81‑1 (MQ=255)                                                      aGGGAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTc                      >  1:31921/1‑80 (MQ=255)                                                       gggAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCaa                    <  1:31888/81‑1 (MQ=255)                                                       gggAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCa                     >  1:31942/1‑80 (MQ=255)                                                        ggAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTcga                    >  1:31911/1‑78 (MQ=255)                                                         gAAACAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCaaaa                  >  2:31914/1‑81 (MQ=255)                                                             cAATCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCt              >  1:31918/1‑81 (MQ=255)                                                              aaTCTCTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTg             >  1:31925/1‑81 (MQ=255)                                                                tctcTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGc            >  1:31938/1‑80 (MQ=255)                                                                tctcTAGACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGc            >  1:31926/1‑80 (MQ=255)                                                                      gACCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTAtt        <  1:31889/78‑1 (MQ=255)                                                                       aCCATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGc    <  1:31890/81‑1 (MQ=255)                                                                        ccATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGcc   <  2:31891/81‑1 (MQ=255)                                                                        ccATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGcc   >  1:31919/1‑81 (MQ=255)                                                                         cATCCTT‑ACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGccc  <  2:31892/81‑1 (MQ=255)                                                                                |                                                                           GGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCC  >  NC_000913/131363‑131516

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