Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,389,688	T→G	G253G (GGT→GGG)	ycjG →	L‑Ala‑D/L‑Glu epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,389,688	0	T	G	95.2%	149.5 / ‑0.1	42	coding (759/966 nt)	ycjG	L‑Ala‑D/L‑Glu epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (23/17);  minor base . (1/1);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGCTGGAGAATTTTATTCATCCGTTGCCGATTTGTGCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGTCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGTTTCG  >  NC_000913/1389564‑1389815                                                                                                                             |                                                                                                                                cgcTGGAGAATTTTATTCATCCGTTGCCGATTTGTGCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAgcgctg                                                                                                                  >  2:2656/1‑140 (MQ=255)      gagaATTTTATTCATCCGTTGCCGATTTGTGCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAgcgctggcgct                                                                                                             <  1:595330/140‑1 (MQ=255)                     ccGTTGCCGATTTGTGCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGcgcg                                                                                              >  2:449627/1‑140 (MQ=255)                         tGCCGATTTGTGCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGca                                                                                          <  1:96034/140‑1 (MQ=255)                           ccGATTTGTGCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGcaca                                                                                        <  1:247663/140‑1 (MQ=255)                              aTTTGTGCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAgg                                                                                     >  1:846222/1‑140 (MQ=255)                              aTTTGTGCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAgg                                                                                     >  1:649143/1‑140 (MQ=255)                                 tgtgCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGttt                                                                                  <  2:355685/140‑1 (MQ=255)                                   tgCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCa                                                                                >  2:661594/1‑140 (MQ=255)                                             aGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCt                                                                      <  2:649143/140‑1 (MQ=255)                                               ttGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTgg                                                                    <  1:76567/140‑1 (MQ=255)                                                  tCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCt                                                                 >  2:643063/1‑140 (MQ=255)                                                  tCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCt                                                                 >  1:158973/1‑140 (MQ=255)                                                    aTACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGc                                                               <  2:636591/140‑1 (MQ=255)                                                           gaGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGATTCAGTCTGATGCTGGGCTGCAtgttgt                                                        >  2:320383/2‑140 (MQ=255)                                                            aGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCAtg                                                            >  2:268086/1‑135 (MQ=255)                                                            aGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCAtg                                                            <  1:268086/135‑1 (MQ=255)                                                            aGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGG‑CTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTgtgt                                                      >  2:892902/1‑140 (MQ=255)                                                              cAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTa                                                     <  1:832073/140‑1 (MQ=255)                                                              cAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTa                                                     <  1:448476/140‑1 (MQ=255)                                                                aTTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTAcc                                                   <  1:449627/140‑1 (MQ=255)                                                                   tGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACctct                                                <  1:2656/140‑1 (MQ=255)                                                                      aGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCa                                                              <  1:297059/123‑1 (MQ=255)                                                                      aGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCa                                                              >  2:297059/1‑123 (MQ=255)                                                                               aaGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGG‑CTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAgcg                                   <  1:892902/140‑1 (MQ=255)                                                                                 gggCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGGGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAgcgc                                  <  1:233441/140‑1 (MQ=255)                                                                                             aTGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCt                      >  2:787623/1‑140 (MQ=255)                                                                                             aTGGTTAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCt                      >  1:13525/1‑140 (MQ=255)                                                                                                 ttAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTgg                    <  2:279351/138‑1 (MQ=255)                                                                                                 ttAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTgg                    >  1:279351/1‑138 (MQ=255)                                                                                                 ttAATATTAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTg                  >  1:61961/1‑140 (MQ=255)                                                                                                       ttAAGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGcgg            >  2:335070/1‑138 (MQ=255)                                                                                                          aGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTc         >  1:712386/1‑140 (MQ=255)                                                                                                          aGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTc         >  1:654148/1‑140 (MQ=255)                                                                                                          aGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTc         >  1:126258/1‑140 (MQ=255)                                                                                                          aGCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGCGCACAAGGTTTCAGACTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTc         >  2:521149/1‑140 (MQ=255)                                                                                                           gCTCGATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCa        >  2:33861/1‑140 (MQ=255)                                                                                                               gATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACTGCTGGTGCCGCAGGTCAGttt    >  2:67177/1‑140 (MQ=255)                                                                                                               gATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGttt    >  1:65613/1‑140 (MQ=255)                                                                                                               gATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGttt    >  2:735727/1‑140 (MQ=255)                                                                                                               gATAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGttt    >  1:479563/1‑140 (MQ=255)                                                                                                                 tAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGTCGCAGGTCAGTTTCg  >  2:847346/1‑140 (MQ=255)                                                                                                                 tAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGTTTCg  <  2:13525/140‑1 (MQ=255)                                                                                                                 tAAAACCGGGGGGCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGTTTCg  <  2:846222/140‑1 (MQ=255)                                                                                                                             |                                                                                                                                CGCTGGAGAATTTTATTCATCCGTTGCCGATTTGTGCTGATGAAAGTTGTCATACTCGTAGCAATTTGAAGGCGCTGAAAGGGCGCTATGAGATGGTTAATATTAAGCTCGATAAAACCGGGGGTCTGACGGAAGCGCTGGCGCTGGCGACTGAAGCGCGTGCACAAGGTTTCAGTCTGATGCTGGGCTGCATGTTGTGTACCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGTTTCG  >  NC_000913/1389564‑1389815

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