Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,752	T→A	G99G (GGA→GGT)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,752	0	T	A	97.1%	107.2 / NA	34	G99G (GGA→GGT)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (11/22);  minor base C (0/1);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAG  >  NC_000913/69507‑69989                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              cacacctcaacccacgcgggggCATATTGCGCCACGGCGCTGTCCGGGCGAGTCACATGAAGGTTTTTTGCCCAGAACCATTCGCTGGAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTc                                                                                                                                                                                                                                          <  2:278539/231‑1 (MQ=255) cacaGCACAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTc                                                                                                                                                                                                                                          <  1:68198/251‑1 (MQ=255)  acaGCTCAACCCACGATGCGGCAGATGGCCCCCCGGCGCTGTCCGGGCGAGTCACATGCAGGATTTTTGCCCAGACCCATTCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCg                                                                                                                                                                                                                                         <  1:113877/251‑1 (MQ=255)      cacaagccccgaggCGGCAGATGGCGCCGCGGCGCTGTCCGGGGGAGTCAAAGGCAGGTTTTTTGCCAAGAACCATTCGTGGGAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCAt                                                                                                                                                                                                                                     <  2:194320/238‑1 (MQ=255)         aaTCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCa                                                                                                                                                                                                                                   >  1:321564/1‑250 (MQ=255)                         gctggcgcCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGAGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCCGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCg                                                                                                                                                                                                                   >  2:90245/5‑250 (MQ=255)                                      cgcTTTCCTGGCGAGTCAAAGGCAGGTTTTTTGCCCAGAACCTTTCGCGGAAATAAAGACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCa                                                                                                                                                                                                     <  1:261349/251‑1 (MQ=255)                                       ggtgtccgggcgagtaaaaggaagtaTTTTGGCCCAAAACCATCCGCGGAAAAAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCaa                                                                                                                                                                                                    <  2:236684/226‑1 (MQ=255)                                               ggctcttcacttgctgtttttttCCCCAGACCCTTTCGCTGGAATAATTACCGCCAATATAGCGGTAGTAGTCATCATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCATTGTGTTCTTTCCTCAATACGAACATCGCTTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCATTCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATc                                                                                                                                                                                            <  1:220687/233‑1 (MQ=255)                                                             aGGATTTTTGCCAAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTgc                                                                                                                                                                              >  1:279493/1‑251 (MQ=255)                                                                    ttcccaaGAACCTTGCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCTCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTg                                                                                                                                                                        <  1:150910/244‑1 (MQ=255)                                                                        ccAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCCCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCg                                                                                                                                                                    >  2:294091/1‑250 (MQ=255)                                                                               catccGCTGGAATAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAg                                                                                                                                                            <  1:176116/247‑1 (MQ=255)                                                                               cATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACACACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTTCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTa                                                                                                                                                              >  1:201244/1‑249 (MQ=255)                                                                                                                         acattGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAAt                                                                                                                  <  1:156250/248‑1 (MQ=255)                                                                                                                           attGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGt                                                                                                                >  1:157442/2‑251 (MQ=255)                                                                                                                             tGCCCGGCGCGTGGCACCAACGGGTACTCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGt                                                                                                                >  1:99085/1‑249 (MQ=255)                                                                                                                             tGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTGTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGGTACCGTCGGCATCAAGCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCGGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTa                                                                                                               <  1:145473/250‑1 (MQ=255)                                                                                                                                  ggCGCGTGGCCCAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTc                                                                                                            >  2:195097/1‑248 (MQ=255)                                                                                                                                  ggCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGGCAACGCCAATCCCGACCACAGCTGCGCGCTGTGCGACGCTAAGCTCTTCAAGCACGGTTTGCAGTGCCGCTTCCATTGACTCAATGTAGTc                                                                                                            <  1:195097/248‑1 (MQ=255)                                                                                                                                                   ggTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACt                                                                                        <  2:112163/251‑1 (MQ=255)                                                                                                                                                            cttcGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTAt                                                                                  <  1:56335/246‑1 (MQ=255)                                                                                                                                                               cGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCggg                                                                             <  1:2122/251‑1 (MQ=255)                                                                                                                                                                gCTTCTTTCACCCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCgggg                                                                            <  2:265946/251‑1 (MQ=255)                                                                                                                                                                gCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCgggg                                                                            <  2:237129/251‑1 (MQ=255)                                                                                                                                                                cttcttt‑cAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGc                                                                           >  1:274612/6‑251 (MQ=255)                                                                                                                                                                     tt‑cAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATc                                                                        <  1:117598/249‑1 (MQ=255)                                                                                                                                                                               aGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTCCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGcc                                                             <  2:24294/251‑1 (MQ=255)                                                                                                                                                                                  gtgGTCTTTCCACAATCCGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGTACTGGTTATTCGGGGCATCACAAAATTGCCCt                                                           <  2:153617/250‑1 (MQ=255)                                                                                                                                                                                                  aCGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGAt                                           <  1:305819/250‑1 (MQ=255)                                                                                                                                                                                                         aagaggtCGGGTGTTGGGCAAACTCCGGGGGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGGGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTGCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCActct                                   <  2:321564/245‑1 (MQ=255)                                                                                                                                                                                                              ttCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCt                              >  2:22803/1‑251 (MQ=255)                                                                                                                                                                                                                        cGGCAAACTCCGGTCGCAGAGCCAGCACGTTACCGTCGGCATCAATCGGAGCGGGCGTCGAGCCGGTACTGTCAACGCCACTCCCGACCACAGCCGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCAGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGAtc                    >  2:214417/1‑251 (MQ=255)                                                                                                                                                                                                                                          gcgcCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAg  <  1:294091/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              CACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAG  >  NC_000913/69507‑69989

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