Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,075,528	C→T	A1119T (GCA→ACA)	putA ←	fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,075,528	0	C	T	94.1%	88.4 / ‑1.8	34	A1119T (GCA→ACA)	putA	fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (15/17);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e-01

GCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGACCATAACC  >  NC_000913/1075293‑1075750                                                                                                                                                                                                                                            |                                                                                                                                                                                                                               gCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCgg                                                                                                                                                                                                                 <  1:32852/251‑1 (MQ=255)                   tCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGCGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACg                                                                                                                                                                                               >  1:49053/1‑251 (MQ=255)                               agcGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATag                                                                                                                                                                                     >  1:134484/1‑249 (MQ=255)                                    ccAGGTGTGGCGTTCACCCGTCGGCCCCGGCAGCAATCGTGTTGTTCCTGCCTGCGCCAGCGCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGAc                                                                                                                                                                              <  2:364655/251‑1 (MQ=255)                                         tgtTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCt                                                                                                                                                                         >  1:199293/1‑251 (MQ=255)                                          gtTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAGCTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTg                                                                                                                                                                        <  1:340758/251‑1 (MQ=255)                                                             cccGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGTCGATTCGCCAGCAGACGGTAGAGGACGAGCGGCACGCCTGCTTTCGGCCCGGTACCggg                                                                                                                                                     >  2:135110/1‑250 (MQ=255)                                                                            gttgtccctgccgtcgcCAGCTCGCCTTATTTCGTACATACCTGCTGCAATTCTGGACGATTTTCTGCCCATTCCCGCAGTGCATTTAGCTGCTGAGTCAAGTCGGCTTTCAACTGCGCATCGACCGGATCCTTTGCATCCGGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTATAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTcgccgc                                                                                                                                       <  2:269575/236‑1 (MQ=255)                                                                                               gCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACaccaacc                                                                                                                   <  1:216057/251‑1 (MQ=255)                                                                                                        atTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGCGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACaccaacagctaggaa                                                                                                           >  2:109344/1‑241 (MQ=255)                                                                                                        atTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACaccaac                                                                                                                    <  1:109344/241‑1 (MQ=255)                                                                                                                    actcctgcatgtcgggccgagttGCTGCCCCTGCCCGCAGTGCATTTAGCGGCTGAGTCAATTCGGCTTTCAACTGCGCATCGCCCGGATACTTGGCATCCTGACGCGCGAGCGTCACTGTCAGCGCCCTTTCCGGGCGATTCGCCATCAGACGGAAGGGATAGAGCTGACCGCCTGCTTGCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACgg                                                                                              <  2:171827/230‑1 (MQ=255)                                                                                                                       ccTGCAATTCTGGACGATTTGCTGCCCATTCTCGCAGTGCATTTAGCGGCAGAGTCAATGCGGCTTTCAACTGCGCATCAACCGAATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCAGGCGATTCGCCAGCAGACGGTAGAGATAGTGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGACCGGCTGCACACCAACCAATGCGCCCCCCATATTACTGtt                                                                                            >  1:356522/1‑250 (MQ=255)                                                                                                                          gCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACa                                                                                        <  1:135110/251‑1 (MQ=255)                                                                                                                             aTTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCTCTGTCAGCGCACTTTCCGGGCGATTCGACAGCAGACGGTAGAGATAGAGCGGCCCGCCTGCTTTCGGCGCGGTCCCGGACAACACTTCGCCGCCGTACGGCTACACGACAACCACTGCGCCCCCCCTACTACGGTTaaattac                                                                                     >  2:108047/1‑246 (MQ=255)                                                                                                                                   gACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCAGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCCCCCTATTACGGTTAACATACAggat                                                                                >  1:289620/1‑248 (MQ=255)                                                                                                                                   gACGATTTGCTGCCCATTCCCACAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGCCCGGATACTTTTCATCCTGACGCGCGAGCGTCACTGTCAGCGCGCTTTCCGGGCGATTCGCGAGCGGACGGTAGAGATAGAGCGGACCGCCTGCTGTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCTCTGCGCCCACCAAATTACGGTTAACATACAGGTTa                                                                               <  1:136894/251‑1 (MQ=255)                                                                                                                                     cGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTAcc                                                                             <  1:315948/251‑1 (MQ=255)                                                                                                                                        tttGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAAc                                                                          >  1:129285/1‑251 (MQ=255)                                                                                                                                             tgcCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCATCATggg                                                                     >  1:145026/1‑251 (MQ=255)                                                                                                                                              gcCCATCCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGc                                                                    <  1:197520/251‑1 (MQ=255)                                                                                                                                                 cATTCCCGCAGAGAATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGAGCATCGCCCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCACACACTGCGCCCACCATATTACGGTTAACATACCGGTTACAAACATgggaca                                                                  >  2:163187/1‑247 (MQ=255)                                                                                                                                                        gCAGTGCATTTAGCGGCTGAGTCAAGGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTg                                                          <  1:9555/251‑1 (MQ=255)                                                                                                                                                                tttAGCGGCTGAGTCAATGCGGCTTTCAAATGCGAATCGACCGGCTACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGTCCCGGTACCAGACAACCCTTCGGCGACGAACGGATGCCGACCAAACACTGAACCCACCATATAACGTTTAACATAAAGGTTACCAACACGGGCCAAGCCAGTGACCGGGGc                                                  >  2:360759/1‑251 (MQ=255)                                                                                                                                                                          gAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGttt                                          <  1:307310/249‑1 (MQ=255)                                                                                                                                                                           aGTCAATGCGGCTTTCAACTGCGCATCGACCGTATACTTTGCATCCGGACGCGCGAGCGTCACTGTCAGCGCACTTCCCGGGCGATTCGCCAGCAGTCGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTtcat                                       <  1:318143/251‑1 (MQ=255)                                                                                                                                                                           aGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTtcat                                       >  1:240115/1‑251 (MQ=255)                                                                                                                                                                            gTCAATTCGGCTTTCAACTGCACCTCGACCGGATACTTTGCATCCTGACGCGCGAGCGGCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGGAGAGATAGAGCGGACCGCCTGCTTTCGGACCGGTACCGGACAACCCTTCGCCGCCGAACCTCTGCACACCAACCACTGCTCCCACCATATAACGGTTAACATAAAGGTTACCAAAATGGACCGCGCCCGATACCTGGGCGATGGTTtcct                                       >  2:40921/1‑248 (MQ=255)                                                                                                                                                                               aaTGCGGATTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTGCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAAt                                   <  1:175733/251‑1 (MQ=255)                                                                                                                                                                                   cGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGTCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATgcgc                               >  1:56739/1‑251 (MQ=255)                                                                                                                                                                                                         cGGATACTTTGCATCCTGCCGCGCGAGCGTCCCTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCTCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGAc         <  1:328142/251‑1 (MQ=255)                                                                                                                                                                                                         cGGATACTTTGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGAc         >  1:199357/1‑251 (MQ=255)                                                                                                                                                                                                           gATACTTTGCATCCTGACGCGCGAGCGACACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGAcc        >  1:128777/1‑250 (MQ=255)                                                                                                                                                                                                                 ttGCATCCTGACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGACCATAAcc  <  1:161/250‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                               GCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGACCATAACC  >  NC_000913/1075293‑1075750

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