Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	364,607	A→T	7.4%	F567I (TTC→ATC)	lacZ ←	beta‑D‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	364,607	0	A	T	7.4%	153.2 / 3.7	54	F567I (TTC→ATC)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base A (26/24);  new base T (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

ATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATT  >  NC_000913/364461‑364753                                                                                                                                                   |                                                                                                                                                   cctgtatgTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGaa                                                                                                                                                    <  3:114040/140‑1 (MQ=255)       cgTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCATCGACTGATCCACCCAGTCCCAGACGAAGCCGccc                                                                                                                                             <  6:189521/148‑1 (MQ=255)       cgTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCAGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGccc                                                                                                                                             <  6:245705/148‑1 (MQ=255)         tGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTg                                                                                                                                           <  3:113250/148‑1 (MQ=255)          gCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGCGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTACCAGACGAAGCCGACCTg                                                                                                                                           >  6:161990/1‑147 (MQ=255)            ggTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTaa                                                                                                                                        <  5:167885/148‑1 (MQ=255)             gTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAgc                                                                                                                                                  >  1:236177/1‑137 (MQ=255)             gTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAgc                                                                                                                                                  <  2:236177/137‑1 (MQ=255)                ggCAAATAACAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGAAAGCCGACCACTGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACg                                                                                                                                     <  3:312061/147‑1 (MQ=255)                ggCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACg                                                                                                                                     <  7:40124/147‑1 (MQ=255)                 gCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACggg                                                                                                                                   >  4:13729/1‑148 (MQ=255)                  cAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                                                                                                                  <  7:241109/148‑1 (MQ=255)                    aataaaagaccGTTCATAAATAACTGGCGATCGTTCGGCGTATCGACAAAATCACCGCCGTAATCCGACCACGGGTTTCCTTTTTCATCATATTTAATAAGAGACTGATCCAACCAGTCCCAGACGAATCCGCACTGTAAACggc                                                                                                                                   >  6:226731/7‑144 (MQ=255)                          agaccGTTCATACAGAACTGGCGATCGTTCTGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGGTTTCATCATATTTAATCAGCGACTGATCCACACAGTCCCAGACGAAGCCGCCCTGTAAACGGTGATACTGAc                                                                                                                          >  3:215243/1‑148 (MQ=255)                          agaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                                                                                                                  <  2:28484/140‑1 (MQ=255)                          agaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                                                                                                                  >  1:28484/1‑140 (MQ=255)                            accGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGa                                                                                                                        <  3:294670/148‑1 (MQ=255)                             ccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGaa                                                                                                                       <  6:205722/148‑1 (MQ=255)                             ccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGaa                                                                                                                       <  3:13729/148‑1 (MQ=255)                                ttCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGccc                                                                                                                                             <  6:219425/123‑1 (MQ=255)                                ttCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGccc                                                                                                                                             >  5:219425/1‑123 (MQ=255)                                     aCAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGc                                                                                                               >  4:217675/1‑148 (MQ=255)                                      cAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAg                                                                                                                                                   >  7:41042/1‑111 (MQ=255)                                      cAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAg                                                                                                                                                   <  8:41042/111‑1 (MQ=255)                                      cAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAATCGGGGATACTGACGAAACGCCTGc                                                                                                               >  3:215540/1‑147 (MQ=255)                                      cAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGc                                                                                                               >  2:366956/1‑147 (MQ=255)                                         aaCTGGCGATCGTTCGGCGTATCGCCAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGt                                                                                                           <  4:229491/148‑1 (MQ=255)                                            tGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCTTCATATTTTATCAGCGACTGATCCACCCAGTCCCAGACGTAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTTCCAGTAtt                                                                                                        >  1:136688/1‑148 (MQ=255)                                             ggCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACggg                                                                                                                                   >  5:343880/1‑120 (MQ=255)                                             ggCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACggg                                                                                                                                   <  6:343880/120‑1 (MQ=255)                                                                  cAAAATCAGCGCCGTAAGCCGTCCACGGGTTGCCTTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGtt                                                                                  <  4:215540/148‑1 (MQ=255)                                                                     aaTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAttt                                                                                                       <  8:356593/124‑1 (MQ=255)                                                                     aaTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAttt                                                                                                       >  7:356593/1‑124 (MQ=255)                                                                              cgTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGTGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTg                                                                      <  5:226731/148‑1 (MQ=255)                                                                                    ccGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGATGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTa                                                                >  8:386809/1‑148 (MQ=255)                                                                                    ccGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGATGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTa                                                                >  2:235190/1‑148 (MQ=255)                                                                                      gACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGATGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTAtt                                                              <  1:235190/148‑1 (MQ=255)                                                                                      gACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGATGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTAtt                                                              <  7:386809/148‑1 (MQ=255)                                                                                          aCGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGc                                                           >  5:381411/1‑147 (MQ=255)                                                                                           cGGGTTGCCGTTTTCATCATATTTAATCATCGACTTATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCaa                                                         <  3:94775/148‑1 (MQ=255)                                                                                                tGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGAt                                                    >  1:346372/1‑148 (MQ=255)                                                                                                      tttCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTTTTTAGCGCAACCGCCAAGACTGTTCCCCATCGCGTGGGCGTATTCGCAAAGGATCAGCgg                                              >  6:16352/1‑148 (MQ=255)                                                                                                           tcatATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGgcg                                           <  5:160981/146‑1 (MQ=255)                                                                                                           tcatATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGgcg                                           >  6:160981/1‑146 (MQ=255)                                                                                                            catATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGt                                        >  4:273324/1‑148 (MQ=255)                                                                                                              tatTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGtct                                      >  3:163892/1‑148 (MQ=255)                                                                                                               atTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGtctc                                     >  1:342316/1‑148 (MQ=255)                                                                                                                      cAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAgg                               >  2:218873/1‑147 (MQ=255)                                                                                                                       aGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTa                             <  2:346372/148‑1 (MQ=255)                                                                                                                        gggtcTGATGCAGCCAGTCTCTGACGAAGCCGGCCTGTAAACGGTGATACTGACGTAACGGCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGTGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAg                            <  4:215243/144‑1 (MQ=255)                                                                                                                            cTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGa                         >  8:236402/1‑147 (MQ=255)                                                                                                                              gATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAg                      >  8:143113/1‑148 (MQ=255)                                                                                                                                  cACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCAt                  <  7:143113/148‑1 (MQ=255)                                                                                                                                     ccAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCAtttt               >  2:259213/1‑148 (MQ=255)                                                                                                                                          cccAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGAt          >  1:8047/1‑148 (MQ=255)                                                                                                                                                  aaGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCAtt  <  1:259213/148‑1 (MQ=255)                                                                                                                                                  aaGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCAtt  <  1:366956/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   ATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATT  >  NC_000913/364461‑364753

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