Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,416,676	A→C	I232L (ATC→CTC)	unknown →	Monofunctional biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,416,676	0	A	C	100.0%	96.2 / NA	35	I232L (ATC→CTC)	unknown	Monofunctional biosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (16/19);  total (16/19)

CAAAGCTGGTTACGTAAAGCCATTGAAGTGCCGACGACCTTGATGCTGGAAACCTTGTGGTCAAAAAAACGCATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGATCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAAAATAAATGACTTAAAATTGGAATCACCATTAATTTAAAACCAATACAACTTAGTGTTGCTATAGGCAATTCTACTTTCCAACTCTTCCAACAACATCTGAAAAGTCAGAAGTGGTTTCAA  >  1/1416433‑1416899                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                cAAAGCTGGTTACGTAAAGCCATTGAAGTGCCGACGACCTTGATGCTGGAAACCTTGTGGTCAAAAAAACGCATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTacac                                                                                                                                                                                                                          <  1:729444/251‑1 (MQ=255) cAAAGCTGGTTACGTAAAGCCATTGAAGTGCCGACGACCTTGATGCTGGAAACCTTGTGGTCAAAAAAACGCATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTacac                                                                                                                                                                                                                          <  1:577236/251‑1 (MQ=255) cAAAGCTGGTTACGTAAAGCCATTGAAGTGCCGACGACCTTGATGCTGGAAACCTTGTGGTCAAAAAAACGCATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTacac                                                                                                                                                                                                                          <  1:270970/251‑1 (MQ=255)                      aTTGAAGTGCCGACGACCTTGATGCTGGAAACCTTGTGGTCAAAAAAACGCATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGAt                                                                                                                                                                                                     <  1:891595/251‑1 (MQ=255)                      aTTGAAGTGCCGACGACCTTGATGCTGGAAACCTTGTGGTCAAAAAAACGCATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGAt                                                                                                                                                                                                     <  1:292351/251‑1 (MQ=255)                             tGCCGACGACCTTGATGCTGGAAACCTTGTGGTCAAAAAAACGCATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAAc                                                                                                                                                                                                <  1:410980/249‑1 (MQ=255)                              gCCGACGACCTTGATGCTGGAAACCTTGTGGTCAAAAAAACGCATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGc                                                                                                                                                                                              >  1:293075/1‑250 (MQ=255)                                                                        cATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGaaa                                                                                                                                                            >  1:657077/1‑242 (MQ=255)                                                                        cATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCaatt                                                                                                                                                   <  1:254501/251‑1 (MQ=255)                                                                          ttttAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGa                                                                                                                                                                            <  1:579319/224‑1 (MQ=255)                                                                          ttttAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCaatta                                                                                                                                                  <  1:94709/250‑1 (MQ=255)                                                                            ttAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCaattaatt                                                                                                                                               >  1:366913/1‑251 (MQ=255)                                                                                    gtATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAGAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCaatta                                                                                                                                                  <  1:57581/240‑1 (MQ=255)                                                                                    gtATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTacac                                                                                                                                                                                                                          >  1:886807/1‑168 (MQ=255)                                                                                    gtATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCGAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCaatta                                                                                                                                                  >  1:441817/1‑240 (MQ=255)                                                                                    gtATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGcagac                                                                                                                                                                                                               >  1:360955/1‑179 (MQ=255)                                                                                    gtATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATcc                                                                                                                                                                      >  1:93529/1‑220 (MQ=255)                                                                                    gtATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGaaa                                                                                                                                                            >  1:464234/1‑230 (MQ=255)                                                                                    gtATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCattta                                                                                                                                       >  1:499840/1‑251 (MQ=255)                                                                                     tATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATc                                                                                                                                                                       <  1:278299/218‑1 (MQ=255)                                                                                      aTTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTAttt                                                                                                                                                                 >  1:47420/1‑223 (MQ=255)                                                                                      aTTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCatttatt                                                                                                                                     >  1:629195/1‑251 (MQ=255)                                                                                                    gAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGaaa                                                                                                                                                            >  1:129105/1‑214 (MQ=255)                                                                                                                ggCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCatt                                                                                                                                         <  1:278804/221‑1 (MQ=255)                                                                                                                  cATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAaaat                                                                                                                       <  1:702447/237‑1 (MQ=255)                                                                                                                  cATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAAAATAAATGACTTAAAAt                                                                                                         >  1:792229/1‑251 (MQ=255)                                                                                                                                  gcagcaAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAAt                                                                                                                                                                                         <  1:798503/155‑1 (MQ=255)                                                                                                                                      caAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAAAATAAATGACTTAAAATTGGAATCACCATTAAtt                                                                                        <  1:15064/248‑1 (MQ=255)                                                                                                                                               tttATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATTCCCCCAGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAAAATAAATGACTTAAAATTGGAATCACCATTAATTTAAAACCAATAc                                                                            <  1:562530/251‑1 (MQ=255)                                                                                                                                                  aTTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAAAATAAATGACTTAAAATTGGAATCACCATTAATTTAAAACCAATACAAc                                                                         <  1:751997/251‑1 (MQ=255)                                                                                                                                                                                              cTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTAttt                                                                                                                                                                 <  1:657283/119‑1 (MQ=255)                                                                                                                                                                                                     ccgccgCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAAAATAAATGACTTAAAATTGGAATCACCATTAATTTAAAACCAATACAACTTAGTGTAGCTATAGGCAATTCTACTTTCCAACTCTTCCAACAACATCTg                       <  1:890124/250‑1 (MQ=255)                                                                                                                                                                                                          gcgTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTAttt                                                                                                                                                                 >  1:159621/1‑107 (MQ=255)                                                                                                                                                                                                               aCCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAAAATAAATGACTTAAAATTGGAATCACCATTAATTTAAAACCAATACAACTTAGTGTTGCTATAGGCAATTCTACTTTCCAActc                                      >  1:354990/1‑225 (MQ=255)                                                                                                                                                                                                                         attattTACAAAGTCAATGCCCCCGGGCTCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAAAATAAATGACTTAAAATTGGAATCACCATTAATTTAAAACCAATACAACTTAGTGTTGCTATAGGCAATTCTACTTTCCAACTCTTCCAACAACATCTGAAAAGTCAGAAGTGGTTTCaa  >  1:793983/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                CAAAGCTGGTTACGTAAAGCCATTGAAGTGCCGACGACCTTGATGCTGGAAACCTTGTGGTCAAAAAAACGCATTTTAGAAGTGTATTTAAATATTGCCGAATTCGGGAATGGCATTTATGGTGTGGAAGCAGCAAGCCGTTTTTATTTCAAAAAGTCGGCAAAAAATTTATCGCAAAATGAAGCCGCACTTTTAGCCGCCGCGTTACCCAATCCGATTATTTACAAAGTCAATGCCCCCGGGATCTACACCAAACGCAGACAGATTTGGATTCAACGCCAAATGAATCAACTGGGAAAATCCTATTTAGAAAATATCAATTAATTCTCATTTATTTTAAAATCTAAAATAAATGACTTAAAATTGGAATCACCATTAATTTAAAACCAATACAACTTAGTGTTGCTATAGGCAATTCTACTTTCCAACTCTTCCAACAACATCTGAAAAGTCAGAAGTGGTTTCAA  >  1/1416433‑1416899

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