Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	364,433	G→A	Q625* (CAG→TAG)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	364,433	0	G	A	93.8%	97.8 / NA	32	Q625* (CAG→TAG)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/12);  minor base C (0/2);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGAC  >  NC_000913/364194‑364669                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             cGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTg                                                                                                                                                                                                                                   <  2:92027/251‑1 (MQ=255)         cTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTc                                                                                                                                                                                                                           >  1:212791/1‑251 (MQ=255)                 aCTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCGTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTg                                                                                                                                                                                                                    <  1:168999/250‑1 (MQ=255)                       aGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATgcgg                                                                                                                                                                                                             <  1:352588/251‑1 (MQ=255)                                 gCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCg                                                                                                                                                                                                   <  1:89213/251‑1 (MQ=255)                                  cccGCCGCTCTCCGGCTGCGGTAGTTCAGGCAGCTCAATCACCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACAGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGCTAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCAGTCAGCCCTGGATGCGGCGTGCGGTCgg                                                                                                                                                                                                  >  2:224470/1‑251 (MQ=255)                                    cGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCa                                                                                                                                                                                                <  1:198295/251‑1 (MQ=255)                                           ctcCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTCTTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGCAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAagacc                                                                                                                                                                                          >  1:57803/1‑250 (MQ=255)                                                   gCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGtt                                                                                                                                                                                  <  2:302338/250‑1 (MQ=255)                                                        aGTTCAGGCAGTTCCATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGCATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATAc                                                                                                                                                                             >  2:177652/1‑250 (MQ=255)                                                        aGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATAc                                                                                                                                                                             >  1:100451/1‑250 (MQ=255)                                                         gTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATAc                                                                                                                                                                             >  1:160957/1‑249 (MQ=255)                                                                 aGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCg                                                                                                                                                                   >  1:70881/1‑251 (MQ=255)                                                                                   ccTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGcc                                                                                                                                                 >  2:63069/1‑251 (MQ=255)                                                                                   ccTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGcc                                                                                                                                                 <  1:78796/251‑1 (MQ=255)                                                                                       gtgGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCaaa                                                                                                                                              <  1:239583/250‑1 (MQ=255)                                                                                         ggAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATc                                                                                                                                           >  2:407569/1‑251 (MQ=255)                                                                                                           aCTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGAcc                                                                                                                          >  1:198680/1‑250 (MQ=255)                                                                                                                         cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGtttt                                                                                                           >  2:73161/1‑251 (MQ=255)                                                                                                                            cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtcat                                                                                                        >  1:257749/1‑251 (MQ=255)                                                                                                                                        caccgcCACCACCCGGTGCAGGGGCTCGTTATCGCTAGGACGGAACAGGTATTCGATGGTCACTTCGATGTTTTGCCCGGATAAACGGAACTGAAAAACCTGCTACTGGTGTTTTGCTTCCGTCAGCGGTGGATGCGGCGTGCGGTCGGAAAAGACCAGCCCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCTCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCa                                                                                            <  2:8070/248‑1 (MQ=255)                                                                                                                                                          cacgAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGAGTTTTGCTTCCGTCAGCGCTGGATGCTGCGTGCGGTCGGCCAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGAGTATCGCCAAAATAACCGCCGTAAGCCGTCCACGGGTTGCCGTTTTCATCAAATTTAACCAGCGACTGATCCACCCAGt                                                                          >  2:230195/4‑251 (MQ=255)                                                                                                                                                                  cgttatcgttatGACGGAACAGGTTTTCGCGGGTCACTGCGATGGTTTGCCCGGATAAACGAAACTGAAAAAACTGCTCCTGGTGTTTTGCTCCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGCCCGTTCATACAGAACTGGCGATCTTTCGGCTTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                   <  1:299155/250‑1 (MQ=255)                                                                                                                                                                     tATCGCTATGACGGAACAGGTATACGCGGGTCCCTTCGAGGGTTTGCCCGGATAAACGGAACGGGAAAAACTGCTACGGGTGTTTTGCTTCCGTCAGCGCTGGAGGCGGCGTGCGGGCGGCAAAGACCAGACCGTTCATACAGAACGGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGaa                                                                <  2:100466/250‑1 (MQ=255)                                                                                                                                                                              gACGGAACAGTTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAACCTGCTACTGGTGTTTTGCTCCCGTCAGCGCTGGATGCGGCGGGCGGTCGGCAAAGACCAGCCCTTTCATACAGAACGGGCGATCGTTGGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGt                                                      <  1:63069/251‑1 (MQ=255)                                                                                                                                                                              gACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGt                                                      >  1:228836/1‑251 (MQ=255)                                                                                                                                                                                         tATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAACATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGAt                                            >  1:29701/1‑250 (MQ=255)                                                                                                                                                                                          aTTCGCTGGTCACTTCGATGGTTTGCCCGGATCAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATAc                                          >  1:255399/1‑251 (MQ=255)                                                                                                                                                                                               cggttcatttCGATGTTTTGCCGGGAAAAACGGAACTGGAAAATCTTCTCCGGTTGTTTTTCGTCCGTCAGCGCGGGATGCGGCGTGCGGTCGGAAAAGACCCGCCCGTTCATCCAGAACTGGCGGTCGTTCGGCGTGCCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACg                                     <  1:352660/243‑1 (MQ=255)                                                                                                                                                                                                        tCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACGGGTGTTTTGCTTCCGCCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGCCCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCt                             <  1:90001/250‑1 (MQ=255)                                                                                                                                                                                                           aTGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTACTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAACCGCCTGcc                          >  2:79092/1‑250 (MQ=255)                                                                                                                                                                                                                                  cTGGAAAAACTGCTACTGGTGTTTTGCTTACGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGAc  >  2:187021/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             CGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGAC  >  NC_000913/364194‑364669

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