Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	973,216	C→T	L62L (CTG→CTA)	elyC ←	envelope biogenesis factor; DUF218 superfamily protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	973,216	0	C	T	100.0%	132.5 / NA	36	L62L (CTG→CTA)	elyC	envelope biogenesis factor; DUF218 superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/15);  total (21/15)

TCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGG  >  NC_000913/973070‑973359                                                                                                                                                   |                                                                                                                                                tcattcaGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGc                                                                                                                                             >  1:156717/1‑151 (MQ=255)   attcaGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCgg                                                                                                                                           <  1:662371/151‑1 (MQ=255)     tcaGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGt                                                                                                                                          <  1:248667/150‑1 (MQ=255)      caGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCg                                                                                                                                        <  2:546534/151‑1 (MQ=255)           atGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGGTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTa                                                                                                                                    >  1:600954/2‑150 (MQ=255)              gCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCgg                                                                                                                                <  2:338615/151‑1 (MQ=255)                  aCTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGc                                                                                                                            >  1:491035/1‑151 (MQ=255)                     gTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCaa                                                                                                                          >  1:268026/1‑150 (MQ=255)                              tAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATaa                                                                                                                 <  1:634116/150‑1 (MQ=255)                                        cTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCgcc                                                                                                       >  1:473932/1‑150 (MQ=255)                                        cTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCgcc                                                                                                       >  2:500348/1‑150 (MQ=255)                                          cGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGTTTGCAAACTCAATAACAAGAGCgcca                                                                                                      <  1:487545/149‑1 (MQ=255)                                           ggTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCgccagcc                                                                                                   <  2:763351/151‑1 (MQ=255)                                                  aCTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAAt                                                                                            >  2:827691/1‑151 (MQ=255)                                                        ggTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCtaat                                                                                      <  2:748445/151‑1 (MQ=255)                                                         gTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCtaata                                                                                     <  2:794106/151‑1 (MQ=255)                                                                  ataACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCtt                                                                             <  1:143127/150‑1 (MQ=255)                                                                  ataACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTc                                                                            >  1:724338/1‑151 (MQ=255)                                                                    aaCCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTccc                                                                          >  2:138941/1‑151 (MQ=255)                                                                        accTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGtt                                                                      >  2:338052/1‑151 (MQ=255)                                                                        accTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGtt                                                                      >  2:30347/1‑151 (MQ=255)                                                                            ccGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGttttt                                                                   <  2:206855/150‑1 (MQ=255)                                                                                 caccacGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGaaa                                                              <  1:748885/150‑1 (MQ=255)                                                                                   ccacGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACg                                                           <  2:42035/151‑1 (MQ=255)                                                                                                 aCTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCAcagc                                              >  2:823975/1‑150 (MQ=255)                                                                                                      ttGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCgcc                                        >  1:747970/1‑151 (MQ=255)                                                                                                          gAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCgccaggc                                    >  2:235842/1‑151 (MQ=255)                                                                                                          gAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCgccaggc                                    >  1:444708/1‑151 (MQ=255)                                                                                                            attattCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCgccaggcc                                   >  2:267448/1‑150 (MQ=255)                                                                                                               attCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAgcg                               >  2:768060/1‑151 (MQ=255)                                                                                                                  ccATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCa                            <  1:335697/151‑1 (MQ=255)                                                                                                                            tAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATcaac                   >  2:39636/1‑150 (MQ=255)                                                                                                                             aGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATcaacaa                 <  2:851728/151‑1 (MQ=255)                                                                                                                             aGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATcaaca                  >  2:796276/1‑150 (MQ=255)                                                                                                                                tGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACAt              >  1:570262/1‑151 (MQ=255)                                                                                                                                            gACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTgg  >  1:582367/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                TCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGG  >  NC_000913/973070‑973359

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