Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,980,583	T→C	Y18C (TAC→TGC)	ygeA ←	amino acid racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,980,583	0	T	C	96.2%	169.1 / ‑1.9	52	Y18C (TAC→TGC)	ygeA	amino acid racemase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (18/32);  total (19/33)
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.62e-01

TCCCCGGTTTTATCCCATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAAT  >  NC_000913/2980451‑2980713                                                                                                                                     |                                                                                                                                   tCCCCGGTTTTATCCCATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAAt                                                                                                                             <  1:348698/140‑1 (MQ=255)      ggTTTTATCCCATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAAtggtgg                                                                                                                        <  1:253529/140‑1 (MQ=255)           tATCCCATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTcc                                                                                                                   >  2:242718/1‑140 (MQ=255)             tCCCATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCa                                                                                                                 >  1:522305/1‑140 (MQ=255)                  ttCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCa                                                                                                            >  1:45050/1‑140 (MQ=255)                  ttCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCa                                                                                                            >  1:372171/1‑140 (MQ=255)                  ttCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCa                                                                                                            >  1:381890/1‑140 (MQ=255)                    ccccGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCAtt                                                                                                          >  1:615755/1‑140 (MQ=255)                          gACGCTGGCACTCTTCTGTTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTccc                                                                                                    >  2:455688/1‑140 (MQ=255)                                       ttctATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTg                                                                                       <  1:82369/140‑1 (MQ=255)                                               cATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCAtt                                                                               <  1:524243/140‑1 (MQ=255)                                                 tGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTcc                                                                             >  1:695247/1‑140 (MQ=255)                                                 tGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTcc                                                                             >  2:713898/1‑140 (MQ=255)                                                  gAAAATCGACGCTATGTAGCAGCACTTGTGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCct                                                                            <  2:344975/140‑1 (MQ=255)                                                  gAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCct                                                                            <  1:314849/140‑1 (MQ=255)                                                  gAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCct                                                                            <  2:45050/140‑1 (MQ=255)                                                  gAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCct                                                                            <  2:372171/140‑1 (MQ=255)                                                  gAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCct                                                                            <  2:615755/140‑1 (MQ=255)                                                  gAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCct                                                                            <  1:203219/140‑1 (MQ=255)                                                  gAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCct                                                                            <  1:729429/140‑1 (MQ=255)                                                   aaaaTCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCctc                                                                           <  1:295636/140‑1 (MQ=255)                                                   aaaaTCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCctc                                                                           <  1:469197/140‑1 (MQ=255)                                                    aaaTCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCctct                                                                          <  1:514573/140‑1 (MQ=255)                                                    aaaTCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCctct                                                                          <  2:300811/140‑1 (MQ=255)                                                            gCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGAt                                                                  <  2:585209/140‑1 (MQ=255)                                                                 gTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGAtgtgtg                                                             <  2:522305/140‑1 (MQ=255)                                                                  tAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATgtgtgt                                                            <  2:577795/140‑1 (MQ=255)                                                                          cTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGcaa                                                    >  1:163638/1‑140 (MQ=255)                                                                          cTTGCGCAGAGTGAAACCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGcaa                                                    >  2:284136/1‑140 (MQ=255)                                                                            tGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGcaaca                                                  <  1:402313/140‑1 (MQ=255)                                                                               gcAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAgg                                               <  2:163638/140‑1 (MQ=255)                                                                                  gagTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAgg                                               <  1:698233/137‑1 (MQ=255)                                                                                  gagTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAgg                                               >  2:698233/1‑137 (MQ=255)                                                                                   agTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGa                                           >  1:379135/1‑140 (MQ=255)                                                                                    gTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGaa                                          <  2:112078/140‑1 (MQ=255)                                                                                    gTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGaa                                          >  1:112078/1‑140 (MQ=255)                                                                                    gTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGaa                                          <  2:695247/140‑1 (MQ=255)                                                                                             accaAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACg                                 >  1:687270/1‑140 (MQ=255)                                                                                                             aTGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAAc                 >  1:665687/1‑140 (MQ=255)                                                                                                             aTGCCTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAAc                 <  1:245067/140‑1 (MQ=255)                                                                                                              tGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACg                >  1:459051/1‑140 (MQ=255)                                                                                                                 cTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGcc             <  2:379135/140‑1 (MQ=255)                                                                                                                 cTTCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGcc             <  1:455688/140‑1 (MQ=255)                                                                                                                  ttCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCa            <  2:665687/140‑1 (MQ=255)                                                                                                                  ttCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCa            <  1:713898/140‑1 (MQ=255)                                                                                                                  ttCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCa            <  1:284136/140‑1 (MQ=255)                                                                                                                   tCATTTATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCaa           >  1:561731/1‑140 (MQ=255)                                                                                                                        tATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGcc      <  2:459051/140‑1 (MQ=255)                                                                                                                        tATCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGcc      <  1:463522/140‑1 (MQ=255)                                                                                                                          tCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCa    <  2:381890/140‑1 (MQ=255)                                                                                                                          tCAAACGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCa    >  1:53705/1‑140 (MQ=255)                                                                                                                            aaaCGATAGCAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAAt  <  1:242718/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                   TCCCCGGTTTTATCCCATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAAT  >  NC_000913/2980451‑2980713

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