Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,761,889	C→A	D388Y (GAT→TAT)	folC ←	tetrahydrofolate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,761,889	0	C	A	94.7%	50.2 / NA	19	D388Y (GAT→TAT)	folC	tetrahydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (8/10);  minor base G (1/0);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATC  >  CP000730/1761756‑1762019                                                                                                                                      |                                                                                                                                   ccATCGCTATCCTGGTTTTATTTTTATAATGCTTCAAAGTCTAACTTTGATCTAACTTCACTTGTTAAATACTGACGACCGGTAATTACTAATGTATCACCTTGTTAATTTTTTATAAATTCAACGTAGTCATATACTAAt                                                                                                                             <  2:359239/140‑1 (MQ=255)        tATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTAttt                                                                                                                      >  2:320508/1‑141 (MQ=255)                   tagtttcaTAATGGCTGAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCa                                                                                                           <  2:401729/141‑1 (MQ=255)                         tataATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCA‑TTTCAAtactac                                                                                                    >  2:483270/1‑141 (MQ=255)                               gCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACa                                                                                               >  2:436405/1‑141 (MQ=255)                                    aaaGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATttc                                                                                          <  2:402794/141‑1 (MQ=255)                                      aGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAGTGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATttctt                                                                                        <  2:40298/141‑1 (MQ=255)                                        tCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCAGTTTCAATACTACCTACAATTtcttct                                                                                      <  2:373410/141‑1 (MQ=255)                                                                         aGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATGTACTAATTGTATTTCATCATTTTCAATACTACCCACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGtt                                                     >  2:367903/1‑141 (MQ=255)                                                                         aGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGtt                                                     >  2:318543/1‑141 (MQ=255)                                                                                          tAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATgcgc                                    <  2:403483/141‑1 (MQ=255)                                                                                                 tCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTaaa                             <  2:170181/141‑1 (MQ=255)                                                                                                     cTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCTTAACGCTTTCGGAAAATTAAATTCAGTTGCATAATACGTATGCGCAATTAAACTTa                         >  2:60462/1‑141 (MQ=255)                                                                                                     cTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTa                         >  2:346753/1‑141 (MQ=255)                                                                                                       tGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTaaa                       <  2:54966/141‑1 (MQ=255)                                                                                                       tGATAATTTTTTAGAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTaaa                       <  2:235480/141‑1 (MQ=255)                                                                                                                    taAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTTCTATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCAtct          <  2:272933/141‑1 (MQ=255)                                                                                                                          cAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGGATGCGCAATTAAACTTAAATGTTTGACCATATCGTTat    >  2:394325/1‑140 (MQ=255)                                                                                                                            aCGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATc  >  2:220239/1‑141 (MQ=255)                                                                                                                                      |                                                                                                                                   TCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATC  >  CP000730/1761756‑1762019

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