Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,186,039	G→A	7.5%	V66I (GTC→ATC)	pepT →	peptidase T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,186,039	0	G	A	7.5%	160.6 / 2.5	53	V66I (GTC→ATC)	pepT	peptidase T
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base A (2/2);  total (26/27)
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.77e-01

CCAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATC  >  NC_000913/1185891‑1186181                                                                                                                                                     |                                                                                                                                               ccAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGCAGGGCACTTTGATGGCGACGTTACCGGCTAACg                                                                                                                                                >  6:371533/1‑149 (MQ=255) ccAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACg                                                                                                                                                >  1:110518/1‑149 (MQ=255) ccAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACg                                                                                                                                                >  8:301769/1‑149 (MQ=255) ccAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACg                                                                                                                                                >  8:299531/1‑149 (MQ=255) ccAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACg                                                                                                                                                >  3:107229/1‑149 (MQ=255)  cAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGt                                                                                                                                               >  4:64142/1‑149 (MQ=255)  cAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGt                                                                                                                                               >  4:189581/1‑149 (MQ=255)    aTCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTcc                                                                                                                                             <  2:110518/149‑1 (MQ=255)    aTCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGCCACTTTGATGGCGACGTTACCGGCTAACGTcc                                                                                                                                             >  1:37903/1‑149 (MQ=255)     tCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTccc                                                                                                                                            <  8:204415/149‑1 (MQ=255)             aGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGCCGaaa                                                                                                                                    >  5:139154/1‑147 (MQ=255)                 gTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATccc                                                                                                                                >  5:289656/1‑149 (MQ=255)                   gagaCAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCgg                                                                                                                              <  7:301769/149‑1 (MQ=255)                               cAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTa                                                                                                                  <  5:371533/149‑1 (MQ=255)                                         ggCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                        <  3:189581/149‑1 (MQ=255)                                         ggCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                        <  7:299531/149‑1 (MQ=255)                                           ccAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGCCCTTAAGTGAGAAGGTCACTTTGATTTCGACGTTACCGGATAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTgg                                                                                                      <  3:221674/149‑1 (MQ=255)                                             aaTGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCTATGTGACCTTAAGTGAGAAGGGCACTGTGATGGCGCCGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGa                                                                                                     >  3:208906/1‑148 (MQ=255)                                                    gTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAc                                                                                             <  4:107229/149‑1 (MQ=255)                                                           cATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTg                                                                                      <  3:64142/149‑1 (MQ=255)                                                             tCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAATGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAt                                                                                                          <  1:36868/127‑1 (MQ=255)                                                             tCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAATGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAt                                                                                                          >  2:36868/1‑127 (MQ=255)                                                             tCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAt                                                                                                          >  4:312851/1‑127 (MQ=255)                                                                gctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGAt                                                                                                    <  5:279095/130‑1 (MQ=255)                                                                gctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGAt                                                                                                    >  6:279095/1‑130 (MQ=255)                                                                gctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCg                                                                                 >  7:32364/1‑149 (MQ=255)                                                                 ctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGc                                                                                  <  6:289656/147‑1 (MQ=255)                                                                  tgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACATCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCt                                                                                                <  8:65722/132‑1 (MQ=255)                                                                  tgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACATCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCt                                                                                                <  8:59692/132‑1 (MQ=255)                                                                  tgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACATCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCt                                                                                                >  7:65722/1‑132 (MQ=255)                                                                  tgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACATCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCt                                                                                                >  7:59692/1‑132 (MQ=255)                                                                   gAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGc                                                                               >  7:123706/1‑148 (MQ=255)                                                                        agcagcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGCCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAt                                                                         <  2:37903/149‑1 (MQ=255)                                                                            gcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGa                                                                     <  8:367053/149‑1 (MQ=255)                                                                             cTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGaa                                                                    >  7:122484/1‑149 (MQ=255)                                                                                gAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATcc                                                                 >  3:111166/1‑149 (MQ=255)                                                                                  agagATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTAt                                                                                                                 <  7:381978/99‑1 (MQ=255)                                                                                  agagATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTAt                                                                                                                 >  8:381978/1‑99 (MQ=255)                                                                                    agaTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCa                                                              >  7:223821/1‑148 (MQ=255)                                                                                          ggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttgtt                                                       <  8:123706/149‑1 (MQ=255)                                                                                          ggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttgtt                                                       <  4:111166/149‑1 (MQ=255)                                                                                          ggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttgtt                                                       <  4:352438/149‑1 (MQ=255)                                                                                              tATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaa                                                    >  2:384398/1‑148 (MQ=255)                                                                                                  aaagagACCTTAAGTGAGAAGGGCACTTTGATTGCGACGATACAGCCTAAAGACCAACGCGATAACCCGGAGATTGGATTAATATCTCATGTGGATACATCACCTGATTTCAGCGGCAAAAATGTGAATCGGCTAATTGTTGAAAACt                                                >  6:263102/6‑148 (MQ=255)                                                                                                   aTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTAt                                              >  5:192509/1‑149 (MQ=255)                                                                                                     gtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcg                                            >  6:368455/1‑149 (MQ=255)                                                                                                          cTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                        >  4:44400/1‑84 (MQ=255)                                                                                                          cTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                        <  3:44400/84‑1 (MQ=255)                                                                                                          cTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGa                                                                     >  2:111018/1‑119 (MQ=255)                                                                                                          cTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGa                                                                     <  1:111018/119‑1 (MQ=255)                                                                                                             aaGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTgg                                      <  6:192509/147‑1 (MQ=255)                                                                                                                    aaGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgc                             <  8:32364/149‑1 (MQ=255)                                                                                                                      gggCACTTTGATGGAGACGTTACAGGCTAACGTCCATGGCGATATCCCGTCGATTGGCTTAATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGatatagagc                           >  6:338128/1‑144 (MQ=255)                                                                                                                         cACTTTGATGGCGACGTACCCGGCCAACGTCCCTGGCGATATCCCGGCGATTGGCTCTACTTCTCATGAGGATACCTCACCGGATTGCAGCGGCAAAAATCCCAATCCGCAAATTGTGGAAAACTATCGCGCGGGCGATATTGCGCTgg                        <  4:208906/149‑1 (MQ=255)                                                                                                                          aCTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTggg                       >  6:161588/1‑149 (MQ=255)                                                                                                                                       cGTTACCGGCTAACGTCCCTGTCGTTATCCCTGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGGTGGGTATCGGCGATGaa          <  5:263102/149‑1 (MQ=255)                                                                                                                                           aCCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGttt      >  1:382063/1‑149 (MQ=255)                                                                                                                                               gCTAAGGTCCCTGGCGATATCCCGGCTATTGTCTTTATTTCTCATTTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTTCGCTGGGTATCGGCGATGAAGTTTTATc  <  5:338128/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               CCAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATC  >  NC_000913/1185891‑1186181

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