Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,174,083	C→T	R315H (CGC→CAC)	coaA ←	pantothenate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,174,083	0	C	T	97.5%	117.8 / NA	40	R315H (CGC→CAC)	coaA	pantothenate kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/21);  minor base A (1/0);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGC  >  NC_000913/4173944‑4174219                                                                                                                                            |                                                                                                                                         ttgcACGCGGAATAGACAAACAGGGGGCTTTATTACTGGAGCAGGATGGAAGAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTTCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTa                                                                                                                                       <  1:427619/140‑1 (MQ=255)     aCGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGAGTAGCCt                                                                                                                                   >  2:543370/1‑143 (MQ=255)      cgcgGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTg                                                                                                                                  <  1:461079/143‑1 (MQ=255)        cgGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGAc                                                                                                                                <  1:838152/143‑1 (MQ=255)         ggAATAGACAAACAGGGGGCTTTATTGCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGAGGAGTATTCGCTACCCTGCAAATTATTTGTGTAGTCTGAcc                                                                                                                               >  1:758672/1‑143 (MQ=255)                           gCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTa                                                                                                             <  1:17544/143‑1 (MQ=255)                              ttattaCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTAGAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCa                                                                                                          <  2:772377/143‑1 (MQ=255)                                  taCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAAGAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCGACCGCATGGCTAGCCCttt                                                                                                      >  1:841033/1‑143 (MQ=255)                                      tGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGt                                                                                                  >  2:690934/1‑143 (MQ=255)                                        aGCAGGATGGAATAATAAAGCCCTGGATGGGCGGTGAAAGACCCCGGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCGCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCa                                                                                                <  1:173329/143‑1 (MQ=255)                                        aGCAGGATGGAATAATAAAACCCCGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTACTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCa                                                                                                <  2:778728/143‑1 (MQ=255)                                               tGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACATCTTCTACCGCAGAATTAGCACTTTTCGTCAGGATTaa                                                                                         >  1:811861/1‑143 (MQ=255)                                                  aataagAGAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGCGCAGTCCTCGCTCCCCTTCAAATTATTTGTGTACTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACt                                                                                      <  1:726866/143‑1 (MQ=255)                                                     aataaAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGgc                                                                                   <  1:594161/143‑1 (MQ=255)                                                       tagatctcTGGATGTGCGGTGAAATATCCCTGCGTAGTGCAGAAGAATGAGAAAGGGGAGTGTTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGgcgc                                                                                 <  2:577169/136‑1 (MQ=255)                                                             ccTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGACGCGCTCAc                                                                           >  1:223169/1‑143 (MQ=255)                                                               tGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGa                                                                         <  2:260761/143‑1 (MQ=255)                                                                 gATGGTCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgt                                                                       <  1:78982/143‑1 (MQ=255)                                                                      gCGGTGAAGTCTCCCTGCGTAGTGCAGAAAAATGAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtaggta                                                                  <  2:824765/143‑1 (MQ=255)                                                                          tGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAat                                                              >  1:50872/1‑143 (MQ=255)                                                                                   ccTGAGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCACCTGCAAATTATGTGTGTAGTCTGACCTCTCCTACCGCATGATTAGCACTTTTGGGCAGGAGTAAACTGGCGAGCTCGCGAGTAGGAAGCATACGTTGCtt                                                     >  2:91066/1‑143 (MQ=255)                                                                                      gCGTAGTGCAGAAAAATAAGAAAGGGGCGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTaa                                                  <  1:370053/143‑1 (MQ=255)                                                                                               aGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTCAACTGGCGCGCTCACGAGTAGGGAGAATAGTTTGCTTTAAGTTCAGCCa                                         >  2:397247/1‑143 (MQ=255)                                                                                               aGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGGGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCa                                         >  1:487413/1‑143 (MQ=255)                                                                                                gAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCCCTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTCGGTAGAATATTTTGCTTTAAGTTCAGcccg                                        >  2:25390/1‑141 (MQ=255)                                                                                                 aaaaaTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTGACGCATAGGTAGACTATTTTGCTTTAACTTCAACCCGt                                       >  2:429307/1‑143 (MQ=255)                                                                                                 aaaaaTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCCGCCAGt                                       >  1:681164/1‑143 (MQ=255)                                                                                                         gAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGCCCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAACCAGTTGATctct                               >  2:691046/1‑143 (MQ=255)                                                                                                          aaaGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTCGGTAGAATGTTGTGCTTTAAGTTCAGCCAGTTGATCTCtt                              >  1:259335/1‑143 (MQ=255)                                                                                                           aaGGGGAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCGTGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGctctctat                             >  2:32802/1‑141 (MQ=255)                                                                                                             ggggAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCAAGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGAAGGTAGAATATTTTACTTTAAGTTCAGCCCGTTATTCTCTTTcc                           >  1:589791/1‑143 (MQ=255)                                                                                                              gggAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTCCCTCATGCTTAGCACTACTCGCCAGGATTAAACCGGCGCGCTCACGAGTCGGTAGAACATTTTGCTTTACGTTAGGCAACTTGATCTCATTCca                          >  2:613948/1‑143 (MQ=255)                                                                                                              gggAGTATTCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCca                          <  2:536191/143‑1 (MQ=255)                                                                                                                      tCGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCa                  <  1:690934/143‑1 (MQ=255)                                                                                                                       cGCTCCCCTGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCAt                 <  1:298076/143‑1 (MQ=255)                                                                                                                               tGCAAATTATTTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCt         <  1:522162/143‑1 (MQ=255)                                                                                                                                   gattattTGTGTAGTTTGTCCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAAt     <  2:499209/142‑1 (MQ=255)                                                                                                                                    attattTGTGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTTATCTCTTTCCACAATGTCATGGCAGTCTTAATc    >  1:168987/1‑143 (MQ=255)                                                                                                                                      tattTGTGTCGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGc  <  1:314826/143‑1 (MQ=255)                                                                                                                                      tattTGTGTAGTCTGACGTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGc  <  2:390622/143‑1 (MQ=255)                                                                                                                                            |                                                                                                                                         TTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGC  >  NC_000913/4173944‑4174219

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