Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,655,443	T→C	T306A (ACC→GCC)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,655,443	0	T	C	96.8%	98.1 / ‑4.0	31	T306A (ACC→GCC)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (14/16);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GGAACTCCGCCAGCGGCAGACGCCAGAACGGTTCGTTCTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTACCGTTGCGATAGGTGATGATAT  >  NC_000913/2655308‑2655577                                                                                                                                        |                                                                                                                                       ggAACTCCGCCAGCGGCAGACGCCAGAACGGTTCGTTCTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAgg                                                                                                                                    >  1:584107/1‑140 (MQ=255)     cTCCGCCAGCGGCAGACGCCAGAACGGTTCGTTCTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTcg                                                                                                                                >  1:306775/1‑140 (MQ=255)     cTCCGCCAGCGGCAGACGCCAGAACGGTTCGTTCTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTcg                                                                                                                                <  1:357271/140‑1 (MQ=255)      tCCGCCAGCGGCAGACGCCAGAACGGTTCGTTCTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTcgc                                                                                                                               <  2:302014/140‑1 (MQ=255)      tCCGCCAGCGGCAGACGCCAGAACGGTTCGTTCTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTcgc                                                                                                                               <  2:253480/140‑1 (MQ=255)       ccGCCAGCGGCAGACGCCAGAACGGTTCGTTCTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTcgcc                                                                                                                              <  2:306775/140‑1 (MQ=255)                gCAGACGCCAGAACGGTTCGTTCTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGACCCGGCGAGGGTCGCCGCATCAATg                                                                                                                     >  2:673277/1‑140 (MQ=255)                                ttcgttcTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCTCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCt                                                                                                     >  2:670041/1‑140 (MQ=255)                                          tGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGc                                                                                           <  1:392035/140‑1 (MQ=255)                                            cgcgGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCat                                                                                         >  2:624930/1‑140 (MQ=255)                                                 cACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCatcaatc                                                                                    >  1:678603/1‑140 (MQ=255)                                                 cACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCatcaatc                                                                                    >  2:99343/1‑140 (MQ=255)                                                 cACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCatcaatc                                                                                    >  2:545732/1‑140 (MQ=255)                                                           cAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCgg                                                                          >  1:584978/1‑140 (MQ=255)                                                               cGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAg                                                                      <  1:209796/140‑1 (MQ=255)                                                                             aTCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCtt                                                        <  1:670041/140‑1 (MQ=255)                                                                                   aaaaCTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCAt                                                  <  1:99343/140‑1 (MQ=255)                                                                                   aaaaCTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCAt                                                  <  1:545732/140‑1 (MQ=255)                                                                                    aaaCTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGACCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATAAGACCATCGGCAAGCACCAGACGCCCTTCGGCATc                                                 <  1:673277/140‑1 (MQ=255)                                                                                               cgcgTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTtcatca                                      >  1:604245/1‑140 (MQ=255)                                                                                                      tAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAAc                               <  2:604245/140‑1 (MQ=255)                                                                                                            ttACCCAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACtttttt                         >  1:622822/1‑140 (MQ=255)                                                                                                               cccAGCGCAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTAcc                      <  1:624930/140‑1 (MQ=255)                                                                                                                    cgcAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTAcc                      >  1:242911/1‑135 (MQ=255)                                                                                                                    cgcAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTAcc                      <  2:242911/135‑1 (MQ=255)                                                                                                                    cgcAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTACCGTTGc                 >  2:49236/1‑140 (MQ=255)                                                                                                                    cgcAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTACCGTTGc                 >  2:638532/1‑140 (MQ=255)                                                                                                                      cAGTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTACCGTTGCGa               <  1:49236/140‑1 (MQ=255)                                                                                                                        gTTTTCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTACCGTTGCGATa             <  1:675350/140‑1 (MQ=255)                                                                                                                         ttttCGCCGCCCCGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTACCGTTGCGATAg            <  2:584978/140‑1 (MQ=255)                                                                                                                                   cccGGCGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTACCGTTGCGATAGGTGATGatat  <  2:678603/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       GGAACTCCGCCAGCGGCAGACGCCAGAACGGTTCGTTCTCCTGCGCGGCACTCGCCAGCAAGCGACCGGCCAGCGCATCGTCAAAACTGAACAGCGCGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTGTTCATCACTTCAACTTTTTTACCGTTGCGATAGGTGATGATAT  >  NC_000913/2655308‑2655577

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