Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	364,526	0	C	A	43.8%	68.0 / 22.1	48	G594C (GGC→TGC)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (24/3);  new base A (0/21);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.81e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTAC  >  NC_000913/364380‑364675                                                                                                                                                       |                                                                                                                                                      ttCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCaaa                                                                                                                                                    >  1:681091/1‑150 (MQ=255)       ggTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCAccg                                                                                                                                             <  1:694074/151‑1 (MQ=255)        gTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCAccgc                                                                                                                                            >  2:466525/1‑151 (MQ=255)         tCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCAccgcc                                                                                                                                           >  1:1004310/1‑151 (MQ=255)             ttCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTa                                                                                                                                        >  1:99405/1‑150 (MQ=255)                gATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGc                                                                                                                                     >  1:294524/1‑150 (MQ=255)                      ttGCCGGGAAAAAGGGAAATGGAAAACCGGCTGCTGGTGTTTTGCTTCC‑GGAGGCGCTGGAGGCGGCGTGCGGGCGG‑AAAACCCAAGACCGTTCAAAAAGAACT‑GGCAATCGTTCG‑GCGTACCGCAAAAACCACCGCCGTAAGCCGACCa                                                                                                                               <  2:923810/150‑1 (MQ=255)                            ggATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑AAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTg                                                                                                                        <  1:36114/151‑1 (MQ=255)                            ggATAAACGGAACTGAAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCGGGATGCGGCGTGCGGGCGG‑CAAAGACCAGCCCGTTCATACAAAACT‑GGCGATCGTTGG‑GCGTATCCCAAAAATCACCGCCGTAAGCCGACCACGGGt                                                                                                                          <  1:959179/149‑1 (MQ=255)                            ggAAAAACGGAACGGGAAAACCGGCTGGTGGTTTTTTGCTTCC‑GCAAGCGGTGGAGGCGGGGTGCGGCCGG‑AAAAGCCCAGCCCTTTCAGAAAGACCT‑GGGGTTGTTTGG‑GCGTAGCCCCAAAATCACCGCCGTAAGCCGACCACGGGtt                                                                                                                         <  2:216619/150‑1 (MQ=255)                                aaaCGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGt                                                                                                                    >  1:1020548/1‑151 (MQ=255)                                    ggAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGACGtttt                                                                                                                 >  2:299741/1‑150 (MQ=255)                                     gAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtc                                                                                                                >  1:541624/1‑150 (MQ=255)                                      aagtGGAAAACCTGCGGGTGGTGTTTTGCTTCC‑GCCAGCGGTGGAGGCGGGGGGCGGGGGG‑AAAAGCCCAGCCGTTTAAAAAAGACCG‑GGGGGTGTTTCG‑GCTTATGACAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtcat                                                                                                              <  2:934469/148‑1 (MQ=255)                                        gggaaaaacCGGCTGTTGGTGTTTTGCTCCG‑TTCAGCCCTGGTGGGGGTTGGCGTGCGG‑GAAAGACAAGCCCGTTATAACAGAGCT‑GGCAATCGTTCG‑GCGTACCGCAAAAAGCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtcatc                                                                                                             <  1:807592/146‑1 (MQ=255)                                          ggAAAACCTGCTGCGGGTTTTTTGCTCCC‑GCCAGCGCTGGAGGGGGCGCGCGGTCGC‑CAAAGACCACCCCGTTCTTACA‑AACTGGGCGGTCGTTCG‑GCGTATCCCAAAAAACGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTcatcat                                                                                                           <  1:57408/150‑1 (MQ=255)                                             aaaaCTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAttt                                                                                                       >  2:255744/1‑151 (MQ=255)                                               aaCTGCTGCTGGTTTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑AAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTaa                                                                                                     <  1:139996/151‑1 (MQ=255)                                                      gctgGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGGCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGGTGCCGTTTTCATCATATTTAATCAGCg                                                                                               >  2:459021/1‑150 (MQ=255)                                                          ggtttttGGTTCCGGTCAGCGTTGGAGGGGGGGTGCGGTCGG‑CAAAGACCGGACCGTTCATAAAGACCG‑GGCGTTGGTTGG‑GCTTATCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTg                                                                                           <  2:334410/148‑1 (MQ=255)                                                              tttGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATcc                                                                                       >  2:421702/1‑150 (MQ=255)                                                              tgtGCTTCC‑GTCAGTGTTGGATGCTTGGTGGGGACGG‑CAAAGACCAGACCGTCCTTACGGAA‑G‑GGCGAGTGTTGGTGGTTCTCGCAAAATTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCa                                                                                                  <  2:847415/137‑1 (MQ=255)                                                                     cc‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTc                                                                               >  1:580347/1‑151 (MQ=255)                                                                     cc‑GTCAGCGCGGGATGGGGCGTGCGGTCGG‑CAAAGCCCAGACCGTTCATACAGACCT‑GGCGATCGTTCG‑GCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCa                                                                                  <  1:388357/148‑1 (MQ=255)                                                                      c‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCACCGACTGATCCACCCAGTcc                                                                              >  2:372609/1‑151 (MQ=255)                                                                        gTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTccc                                                                             >  1:213073/1‑151 (MQ=255)                                                                          cAGCGGGGGATGCGGCGTGGGGCCGG‑CAAAGCACAGCCCGTTCATACAGACCG‑GGCGATCGTTCG‑GTGTACCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCa                                                                            <  2:342505/150‑1 (MQ=255)                                                                            gcgcTGGAGGCGGCGTGGGGTCGG‑CAAAGCCCAGACCTTTCATACAGACCT‑GGCGATCTTTCG‑GCTTACCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAg                                                                           <  2:491605/149‑1 (MQ=255)                                                                                 ggtgggggcgTGC‑GTCGGACAAACCCAAGACCGTCCATAAAGACCG‑GGCGATGGTTCG‑GCGTATCACAAAATTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAg                                                                           <  2:700877/138‑1 (MQ=255)                                                                                 ggAGGCGGCGTGGGGCCGA‑AAAAGCCCAGCCCTTTCATACAGACCT‑GGGGATCGTTCG‑GCTTATCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACg                                                                        <  2:541543/147‑1 (MQ=255)                                                                                            gCGGCCGG‑AAAGGCCCAGACCTTTCAAACAGACCT‑GGCGATCGTCCG‑GCGATCCGCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTa                                                           <  2:710262/149‑1 (MQ=255)                                                                                               gTCGG‑CAAAGCCCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgg                                                      <  1:205276/151‑1 (MQ=255)                                                                                                   g‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGAt                                                  >  1:917100/1‑151 (MQ=255)                                                                                                          accagaccGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGAc                                            >  2:246195/1‑151 (MQ=255)                                                                                                                cctttCAAAAAGAACT‑GGCGATCGTGCG‑GGTTAGCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACg                                      <  2:59717/148‑1 (MQ=255)                                                                                                                 cGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGc                                     >  2:366586/1‑151 (MQ=255)                                                                                                                      aTACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGcc                                >  1:822256/1‑151 (MQ=255)                                                                                                                       aaacatacctgtggATCGTTCG‑GCGTACCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCt                                   <  2:942341/135‑1 (MQ=255)                                                                                                                          aGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGt                             >  1:668468/1‑150 (MQ=255)                                                                                                                           gAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAt                           >  1:107904/1‑151 (MQ=255)                                                                                                                                 ggCGATGGTTGG‑CCGTACCGCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAttt                         <  2:845513/148‑1 (MQ=255)                                                                                                                                 ggCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGc                      >  1:760852/1‑151 (MQ=255)                                                                                                                                  gCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCg                     >  2:81309/1‑151 (MQ=255)                                                                                                                                   ctttCGTTCG‑CCTTTTCGCAAAATTCACCGCCGTAAGCCTACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCg                     <  1:757036/147‑1 (MQ=255)                                                                                                                                    gATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGa                    >  1:115907/1‑150 (MQ=255)                                                                                                                                    gATCGTTCG‑GCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGaa                   <  1:448022/151‑1 (MQ=255)                                                                                                                                         ttCG‑GCTTTCCGCCAAATTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCg               <  2:354385/150‑1 (MQ=255)                                                                                                                                           cG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGcc             >  2:438983/1‑150 (MQ=255)                                                                                                                                                      ccAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTAc  >  1:338651/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      TTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCC‑GTCAGCGCTGGATGCGGCGTGCGGTCGG‑CAAAGACCAGACCGTTCATACAGAACT‑GGCGATCGTTCG‑GCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTAC  >  NC_000913/364380‑364675

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: ‑