Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	364,526	0	C	A	43.6%	49.1 / 16.2	39	G594C (GGC→TGC)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (16/6);  new base A (0/17);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 2.06e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACC  >  NC_000913/364382‑364676                                                                                                                                                 |                                                                                                                                                       cGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATc                                                                                                                                                  <  1:394960/151‑1 (MQ=255)     ggTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCAccg                                                                                                                                              <  1:498272/151‑1 (MQ=255)        cACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCAccgcc                                                                                                                                            <  1:305814/150‑1 (MQ=255)         aCTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGt                                                                                                                                          >  2:200473/1‑151 (MQ=255)                 gttttGCCCGGATAAACGGAACTGAAAAAACTGCTGCTGGTGTTTTGCTCCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCg                                                                                                                                    <  1:468155/147‑1 (MQ=255)                    ttGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCAc                                                                                                                               >  1:386021/1‑151 (MQ=255)                      gCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACgg                                                                                                                             >  1:665652/1‑151 (MQ=255)                                 cGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGtttt                                                                                                                  >  1:328651/1‑151 (MQ=255)                                  ggaatggaaaaaaCTGCTGGTGGTGTTTTGCTCCCGTCAGCGCTGGATGCGGCGTGCGGTCGGAAAAGACCAGCCCGTTCATACAGAACTGGCGATCGTTCGGCGTATCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGttt                                                                                                                   <  1:485260/141‑1 (MQ=255)                                        ggAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTcatcat                                                                                                            >  1:757825/1‑150 (MQ=255)                                          aaaaaCTGCTGCTGGTGTTTTGCTTCCTTCGGCGCTGGATGCGGCTGGCGGTCGGCAAAGCCCAGCCCGTCCTAACAGAACGGGCGATCGTTCGGCTTTTCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtcatca                                                                                                             <  2:662948/147‑1 (MQ=255)                                           aaacCTGCTGCGGGTTTTTTGTTCCCTTCAGCGCTGGATGCGGGGTGCGGTCGGCAAAGCCCAGACCGTTCATACAGAACTGGCGATCGTTGGGCGTACCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAttt                                                                                                        <  1:753616/151‑1 (MQ=255)                                             aaCTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTa                                                                                                       >  1:159361/1‑150 (MQ=255)                                             aaCTGCGGTTGGTGTTTTGCTTCCGTCGGCGCGGGAGGGGGCGTGCGGTCGGCAAAGACCGGACCGTTCATAAAAAAGTGGGGTTCTTTCGGCTTTTCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTaa                                                                                                      <  2:540200/151‑1 (MQ=255)                                              aCTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAAt                                                                                                     >  1:303848/1‑151 (MQ=255)                                                 gctgctgGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATc                                                                                                    >  2:765930/1‑149 (MQ=255)                                                      tgGTGTTTTGATTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACt                                                                                             >  2:694105/1‑151 (MQ=255)                                                       ggTTTTTTGCTCCCGTCAGCGCGGGATGCGGCGTGCGGTCGGCAAAGCCCAGACCTTTCATACAGACCTGGCGATCGTTGGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTg                                                                                            <  1:773563/151‑1 (MQ=255)                                                       ggTTTTTTGCTCCCGTCACCGCTGGAGGCGGCGGCCGGTCGCCAAAGACCAGCCCTTCCATCCAGACCTGGCGATCGTTCGGCGTATCCCAAAATTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACt                                                                                             <  2:611389/150‑1 (MQ=255)                                                           ttttGCTCCCGTCAGCGCTGGAGGCGGCGTGCGGTGGGAAAAGACCAGACCGTTCATACAGAACTGGCGTTCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGAt                                                                                          <  1:702724/149‑1 (MQ=255)                                                                     gTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTccc                                                                              >  1:621759/1‑151 (MQ=255)                                                                              ggATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCAAACAGAACTGGCGATCGTGCGGCGTATCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAg                                                                      <  2:215181/150‑1 (MQ=255)                                                                              ggAGGCGGCGGGCGGTCGGCAAAGCCCAGCCCGTCCAAACAGACCGGGCGATCGTTGGGCGTCTCCCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAGgc                                                                     <  2:271158/151‑1 (MQ=255)                                                                                         gggggcGGCAAAGCCCAGCCCGTTCATAAAGAACTGGCCATCTGTGGGCGTATCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTaa                                                           <  2:452985/145‑1 (MQ=255)                                                                                               ggCAAAGCCCAGACCGTTCATACAAAACGGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgg                                                       <  2:52143/148‑1 (MQ=255)                                                                                                 cAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGa                                                    >  2:471726/1‑149 (MQ=255)                                                                                                      accagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGAc                                             >  1:275466/1‑151 (MQ=255)                                                                                                       ccagcccTTTCAAAAAGACGTCGCGATCGTCCGGCGTACCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACg                                            <  2:271968/151‑1 (MQ=255)                                                                                                       ccagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGAc                                             <  1:241650/150‑1 (MQ=255)                                                                                                            ccGTTCATACAGAACTGGCGATCGTTCGGGGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGaa                                          <  1:665245/148‑1 (MQ=255)                                                                                                            ccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGCCTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACg                                       >  2:190784/1‑151 (MQ=255)                                                                                                              gTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACAGCCGTAAGCCGACAACGGGTTGCCGTTTTCATCATATTTAATCAGCGACCGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGcc                                     >  2:79183/1‑151 (MQ=255)                                                                                                                  aTACAGACCTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGc                                  <  1:446970/150‑1 (MQ=255)                                                                                                                  aTACAGACCTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGc                                  <  1:417697/150‑1 (MQ=255)                                                                                                                      agagCTGCCGTTCGTCCGGGGTATCCAAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCa                                <  1:772144/148‑1 (MQ=255)                                                                                                                         aCTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAttt                          >  1:467274/1‑151 (MQ=255)                                                                                                                            ggCGATCGTCCGTCGTTTCGCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCa                                <  2:592505/142‑1 (MQ=255)                                                                                                                                         cGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCa             >  2:20357/1‑148 (MQ=255)                                                                                                                                              ccacAAAAGCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCACCATATTTAATCAGCGACTGTTCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAACCGCCTGCCAGTATTTAGCGAAGCCGCCAAGACTg       <  2:200139/146‑1 (MQ=255)                                                                                                                                              cGCAAAATTCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCCCCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGAc         <  1:200515/147‑1 (MQ=255)                                                                                                                                                 cAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTAcc  >  2:198903/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                       CGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACC  >  NC_000913/364382‑364676

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