Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,551,012	G→T	L95L (CTC→CTA)	gntP ←	fructuronate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,551,012	0	G	T	98.2%	193.7 / ‑3.9	56	L95L (CTC→CTA)	gntP	fructuronate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (29/26);  minor base C (0/1);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCGAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCATAGTGTG  >  NC_000913/4550876‑4551146                                                                                                                                         |                                                                                                                                       tCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCtaga                                                                                                                                     <  1:61276/140‑1 (MQ=255)        cGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCca                                                                                                                              >  2:13278/1‑140 (MQ=255)                    ttCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACAGCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCCAGACGCGCCAGCAGCGTATgcg                                                                                                                  <  2:622229/140‑1 (MQ=255)                                aacaataacCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAg                                                                                                      <  2:575861/140‑1 (MQ=255)                                 acaataacCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGc                                                                                                     <  1:1578/140‑1 (MQ=255)                                       acCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGcc                                                                                               >  1:292786/1‑140 (MQ=255)                                          aGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGa                                                                                            >  2:457093/1‑140 (MQ=255)                                            cggcgCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGt                                                                                          >  2:703867/1‑140 (MQ=255)                                                   aaCATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCAt                                                                                   <  1:75204/140‑1 (MQ=255)                                                            aaGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTAcc                                                                          <  2:339160/140‑1 (MQ=255)                                                              ggCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGa                                                                        >  2:75310/1‑140 (MQ=255)                                                               gCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGTGTCGACCATCAATTTACCGAt                                                                       >  1:113460/1‑140 (MQ=255)                                                                ccACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCAt                                                                                   >  1:258429/1‑127 (MQ=255)                                                                ccACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCAt                                                                                   <  2:258429/127‑1 (MQ=255)                                                                ccACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATg                                                                      >  2:376381/1‑140 (MQ=255)                                                                   cTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGAcc                                                                   >  1:287546/1‑140 (MQ=255)                                                                    ttCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCg                                                                  >  1:314655/1‑140 (MQ=255)                                                                        tAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCAcc                                                              <  2:292786/140‑1 (MQ=255)                                                                         aGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCg                                                             >  1:640044/1‑140 (MQ=255)                                                                           aaCATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCGTCAATTTACCGATGACCGCACCGaa                                                           >  1:277486/1‑140 (MQ=255)                                                                           aaCATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGaa                                                           <  2:214074/140‑1 (MQ=255)                                                                                cGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGa                                                                     >  2:191248/1‑125 (MQ=255)                                                                                cGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGa                                                                     <  1:191248/125‑1 (MQ=255)                                                                                cGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAAcacca                                                      >  1:56756/1‑140 (MQ=255)                                                                                cGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAAcacca                                                      >  1:239094/1‑140 (MQ=255)                                                                                cGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAAcacca                                                      >  2:108212/1‑140 (MQ=255)                                                                                cGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAAcacca                                                      >  1:522326/1‑140 (MQ=255)                                                                                cGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAAcacca                                                      >  2:687207/1‑140 (MQ=255)                                                                                cGCCAGACCAAAAATCAGGCCGATGATAAACACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAAcacca                                                      <  2:719650/140‑1 (MQ=255)                                                                                          aaaaTCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGaa                                                           <  1:350358/125‑1 (MQ=255)                                                                                          aaaaTCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGaa                                                           >  2:350358/1‑125 (MQ=255)                                                                                          aaaaTCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGaa                                                           >  2:504604/1‑125 (MQ=255)                                                                                          aaaaTCAGGCCGAGGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTCGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGaa                                                           <  1:504604/125‑1 (MQ=255)                                                                                                  gCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTccc                                    >  2:52698/1‑140 (MQ=255)                                                                                                   ccGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTcccc                                   <  1:687207/140‑1 (MQ=255)                                                                                                   ccGATGATAATCACCGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTcccc                                   <  2:640044/140‑1 (MQ=255)                                                                                                                ccGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCg                      >  2:351892/1‑140 (MQ=255)                                                                                                                 cGACAGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGa                     >  2:20320/1‑140 (MQ=255)                                                                                                                     aGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATCGCCAGTTCCCCCAGCGTGTTGCCGAAGcc                 <  2:113460/140‑1 (MQ=255)                                                                                                                     aGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGcc                 <  2:98583/140‑1 (MQ=255)                                                                                                                     aGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGcc                 <  2:291845/140‑1 (MQ=255)                                                                                                                     aGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGcc                 <  1:13278/140‑1 (MQ=255)                                                                                                                     aGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGcc                 <  1:20320/140‑1 (MQ=255)                                                                                                                     aGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGcc                 <  1:351892/140‑1 (MQ=255)                                                                                                                     aGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGcc                 <  1:376381/140‑1 (MQ=255)                                                                                                                     aGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGcc                 <  2:239094/140‑1 (MQ=255)                                                                                                                     aGCTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGcc                 <  2:314655/140‑1 (MQ=255)                                                                                                                       cTGTACATAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCg               <  1:457093/140‑1 (MQ=255)                                                                                                                             aTAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCa         <  2:277486/140‑1 (MQ=255)                                                                                                                              tAGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCAt        >  2:705337/1‑140 (MQ=255)                                                                                                                               aGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTtata       >  1:333931/1‑140 (MQ=255)                                                                                                                               aGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCATa       >  1:624662/1‑140 (MQ=255)                                                                                                                               aGCGCAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCATa       <  1:339304/140‑1 (MQ=255)                                                                                                                                 cgcAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCATAgt     >  1:498812/1‑140 (MQ=255)                                                                                                                                 cgcAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCATAgt     >  1:406493/1‑140 (MQ=255)                                                                                                                                  gcAGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCATAgtg    <  2:624662/140‑1 (MQ=255)                                                                                                                                    aGACCTAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCATAGtgtg  >  2:291585/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       TCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGGCCGATGATAATCACCGACAGCTGTACATAGCGCAGACCGAGACGCGCCAGCAGCGTATGCGCTATCTGGTGAGCCGCGCCGGAGTCGACCATCAATTTACCGATGACCGCACCGAACACCACGATGATAGCCAGTTCCCCCAGCGTGTTGCCGAAGCCCGCTTTCATAGTGTG  >  NC_000913/4550876‑4551146

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