Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,900	C→G	V146V (GTG→GTC)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,900	0	C	G	96.7%	90.5 / NA	31	V146V (GTG→GTC)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base G (15/14);  minor base T (0/1);  total (16/15)
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.91e-01

GCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACC‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTC  >  NC_000913/67653‑68142                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    gcgcGTTAACATCGCCGTCGCTGATGGAGTTCCCCCCCTGCCCCAGACCGCCAACCCCCCAGTTATTGAGGGAGAAACCGAACTTGATCTGTGCGGCAACTTTACCCCCATCGGGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAacg‑acg                                                                                                                                                                                                                                                 <  1:266446/251‑1 (MQ=255)        aaCATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGtt                                                                                                                                                                                                                                          >  1:164562/1‑251 (MQ=255)        aaCATCGCCGTCGCTGATAACGTCCACCACCTGCAGCAGATCGCCAACCGAACAGGTATTGACGGAGACACCGAACTTGATCCGTGCGGCAACTTTATCACCATCGGTGCCCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGCCGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGCCTTGTTTATCCTGCCAGTGACTGCTAACG‑ACGGCCTGtt                                                                                                                                                                                                                                          >  2:165524/1‑251 (MQ=255)            tCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTg                                                                                                                                                                                                                                      >  1:97540/1‑251 (MQ=255)               ccGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACg                                                                                                                                                                                                                                   >  1:176291/1‑251 (MQ=255)                     gggagggAGTCCACCACCTCCCCCAGACCGCCAACCGCCCAGGTATTGAGGGAGAAACCGAACTTGATCTGTGCGGCAACTTTTTCACCATCGTTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATAcg                                                                                                                                                                                                                             <  2:265479/246‑1 (MQ=255)                                                              ggTATTGACGGAGAAACCGAACTTGCTCTGTGCGGCAACGTTATCACAATCGGAGACCGCCACTTAACGCATGTTATCGCCCAAACGGCAGACTTTCAGATGAAGGGTATCCTGTTTAGCCACCGCCTGACGCATACAGGAGCAGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCt                                                                                                                                                                                    >  1:263205/1‑251 (MQ=255)                                                                      cGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcagg                                                                                                                                                                            >  2:192663/1‑251 (MQ=255)                                                                          gaAACCGCACTTGATCCGTGAGGCAACTTTATCACCCTCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGCATCATGTTTAGCAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCAGAACTCGCGACCGCCATGTGCAGTCTggttcaggttc                                                                                                                                                                         >  2:157240/1‑250 (MQ=255)                                                                                                     tttaccccctccgtgaaCCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCAc                                                                                                                                              <  1:197152/234‑1 (MQ=255)                                                                                                                            cTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGgtg                                                                                                                      >  2:216802/1‑251 (MQ=255)                                                                                                                                    aTGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTg                                                                                                              <  1:237129/251‑1 (MQ=255)                                                                                                                                         aTCGCCAAAACGGCAGACTTTCAGATGACTGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTTTCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTgcagc                                                                                                          <  2:238944/250‑1 (MQ=255)                                                                                                                                            gCCAAAGCGGCAGACTTTCAGATGACGGGTATCCTGTTAAGAAACCGCCTGATGCATCCAGGAGCCGATACTCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGGCCGCCATGTGCAGTCTGGTTCAGTTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACg                                                                                                      <  2:190014/251‑1 (MQ=255)                                                                                                                                            ccaaaaacagCAGACGTGCAGAGGACGGGTATCCTGTTGAGAAACCGCCTGATGCATCAAGGAGCCCAGACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGCAACT‑ACGGCATGTTGCTGACGCATACGCTCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAACGCCGCGTTGAACTTGGTGTGGAATTGCAGCAACg                                                                                                      <  2:190023/242‑1 (MQ=255)                                                                                                                                                 aacGGCAGACTTTCAGATGACGGGTATCCTGTTTCGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTttg                                                                                                 >  2:96615/3‑251 (MQ=255)                                                                                                                                                     gCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGCCGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGa                                                                                             <  2:18233/251‑1 (MQ=255)                                                                                                                                                              cAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTc                                                                                     >  1:83259/1‑250 (MQ=255)                                                                                                                                                               aGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAg                                                                                   <  2:247833/251‑1 (MQ=255)                                                                                                                                                               aGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAg                                                                                   <  2:176291/251‑1 (MQ=255)                                                                                                                                                                         tATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTTTGGTTCAGGTTCATAAAGTCCATATCGATACTGTccc                                                                                                                                                <  2:115826/180‑1 (MQ=255)                                                                                                                                                                         tATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTccc                                                                                                                                                >  1:115826/1‑180 (MQ=255)                                                                                                                                                                           ccctgttcagaacccGCCTGACGCATCCAGGAGCCGATGCGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCACACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGTTTCAGGTTCATAAAGTCCATAGCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATAGTCGGGCCGTTGATCca                                                                       <  2:35307/238‑1 (MQ=255)                                                                                                                                                                            ccTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCca                                                                       <  1:218459/250‑1 (MQ=255)                                                                                                                                                                                   agaaaCCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTg                                                                <  1:218903/246‑1 (MQ=255)                                                                                                                                                                                           ccTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCgggg                                                         >  1:22528/1‑249 (MQ=255)                                                                                                                                                                                           ccTGACGCATCCAGGAGACGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGG‑AACCGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGCCCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCgggg                                                         >  1:150669/1‑249 (MQ=255)                                                                                                                                                                                                  cATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGgtgt                                                >  1:85749/1‑251 (MQ=255)                                                                                                                                                                                                         gAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGcc                                          >  2:60719/1‑250 (MQ=255)                                                                                                                                                                                                                                         ccAGTGACCGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAAt         >  1:155035/1‑251 (MQ=255)                                                                                                                                                                                                                                           aGTGACCGGTAGCG‑ACGTCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACGCTCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATtc       <  2:197382/251‑1 (MQ=255)                                                                                                                                                                                                                                                ccGGTAACG‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATtcgcgtc  <  1:167210/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACC‑ACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTC  >  NC_000913/67653‑68142

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