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

CTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAA‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTT  >  NC_000913/364384‑364673                                                                                                                                                 |                                                                                                                                                    cTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAA‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCAc                                                                                                                                             >  1:823230/1‑151 (MQ=255)   ggTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAA‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCAc                                                                                                                                             <  2:440082/149‑1 (MQ=255)            gtgggTTTGCCGGGA‑AAAAGGAACCGGAAAAACCT‑GCTGCTGGTTTTTTGCTCCCGTCAGCGCGGGATGGGGCGTGCGGTGGGAAAAGACCAGCCCGTTCATAAAGAACGGGCGATCGTTCGGCTTACCGCCAAAATCACCGCCGTAAGc                                                                                                                                   <  2:471492/147‑1 (MQ=255)                   tGCCCGGATAAACGGAACTGGAAAAA‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACg                                                                                                                           >  1:944247/1‑151 (MQ=255)                    gCCCGGATAACCGGAACTGGAAAAA‑CT‑GCTGCTGGTTTTTTGCTTCCGTCAGCGCTGGATGCGGGGTGCGGTCGGAAAAGCCCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCAc                                                                                                                            <  1:633484/149‑1 (MQ=255)                     cccGGATAAACGGAACTGGAAAAA‑CT‑GCTGCTGGTTTTTTGCTCCCGTCAGCGCGGGAGGCGGCGTGCGGTCGGCAAAGACCAGACCTTTAATACAGACCTGGCGATCTTTCGGCGTATCGCAAAAACCACCGCCGTAAGCCGACCACg                                                                                                                           <  1:353271/149‑1 (MQ=255)                        ggAAAAACGGACCTGAAAAAA‑CT‑GCTGCTGGTGTTTTGCTCCCGTCAGCGCGGGATGCGGCGTGCGGTGGGAAAAGACCAGCCCGTTCATACAAAACTGGCGATCTTTCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGtt                                                                                                                       <  1:449523/150‑1 (MQ=255)                        gaaTAACCGG‑ACTGGAAAAA‑ATGGCTGCGGGTGTTTTGCTTCCGTCAGGGTTGGAGGCGGTGTGCGGTGGACAAAGCCAAGCCCGTCCATACAGAACGGCCGCTGGTTGGGGGTAGCGCAAAAATCCCCGCCGTAAGCCGACCACGGGtt                                                                                                                       <  2:91465/148‑1 (MQ=255)                            aaaCGGAACTGAAAAAC‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGAAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGt                                                                                                                  <  1:719775/151‑1 (MQ=255)                             aaCGGAACTGGAAAAA‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGtt                                                                                                                 >  1:507497/1‑151 (MQ=255)                                   aCTGGAAAAA‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtcatc                                                                                                           >  1:63595/1‑151 (MQ=255)                                     tGGAAAAA‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTcatcat                                                                                                         >  1:1006579/1‑151 (MQ=255)                                      ggAAAAA‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTcatcat                                                                                                         >  2:120595/1‑150 (MQ=255)                                         aaat‑CT‑GCTGCGGTTTTTTTGTTTCCGTCAGCGCGGAGTGCGGCGTGCGGCCGGCAAAGACCAGCCCTTTCATACAGAACTGGCTATCGTTCGGGGTACCGCAAAAATAACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAttt                                                                                                     <  1:365498/151‑1 (MQ=255)                                          aaa‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAttt                                                                                                     >  2:120638/1‑150 (MQ=255)                                                  ctgctgGTTTTTTGCTCCCGTCAGCGCGGGAGGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCa                                                                                                <  1:573613/149‑1 (MQ=255)                                                   tgcggGTGTTTTGCTTCCGTCAGCGCGGGAGGCGGCGTGCGGCGGGCAAAGACAAGCCCGTTCAAACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCa                                                                                                <  1:745098/148‑1 (MQ=255)                                                       ggTGTTTTGCTTCCGTCAGCGCTGGAGGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACt                                                                                          <  1:201209/150‑1 (MQ=255)                                                                   ccGCCAGCGCTGGAGGCGGCGTGCGGTCGGAAAAGACCAGCCCGTTCATACAAAACGGGCGATCGTTCGGCTTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAg                                                                               <  2:50937/149‑1 (MQ=255)                                                                    cgccggcgcGGGATGCGGCGTGCGGTCGGCAAAGACCAGCCCTTTCATACAGAACTGGCGATCGTTCGGCTTATCGCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTcc                                                                            <  1:134982/146‑1 (MQ=255)                                                                     gTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTccc                                                                           >  1:82207/1‑151 (MQ=255)                                                                       cAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTccc                                                                           >  1:223119/1‑149 (MQ=255)                                                                            cTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGa                                                                     >  1:761331/1‑150 (MQ=255)                                                                               gATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAgcc                                                                 >  2:11565/1‑151 (MQ=255)                                                                                       gTGCGGTCGGAAAAGCCCAGCCCGTTCATAAAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGt                                                          <  1:739884/150‑1 (MQ=255)                                                                                           ggTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACg                                                     >  2:26658/1‑151 (MQ=255)                                                                                                  aaaGACCAGCCCGTTCACACAGATCTGGCGATCGTTCGGTTTACCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATa                                               <  2:545595/150‑1 (MQ=255)                                                                                                   aaGACCAGACCGTTCATACAAACCGGGCGTTGGTTGGGCTTACCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACt                                             <  2:722894/151‑1 (MQ=255)                                                                                                        cagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGCATACTGACg                                         >  1:390893/1‑150 (MQ=255)                                                                                                            ccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAAc                                     <  1:128930/150‑1 (MQ=255)                                                                                                               ttCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCt                                 >  2:143033/1‑151 (MQ=255)                                                                                                                  aTACAGAATGGGCGTTCGTTCGGCTTATCCCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCt                                 <  2:507361/148‑1 (MQ=255)                                                                                                                        accgggCGTTCGTTGGGCTTACCGCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAtt                        <  2:944014/147‑1 (MQ=255)                                                                                                                            ggCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTa                      <  1:106198/149‑1 (MQ=255)                                                                                                                             gCGGTCGTTCGGGGTCCCGGAAAAACAACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCGGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGc                    <  1:425844/150‑1 (MQ=255)                                                                                                                                  cGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAAcc              >  2:598239/1‑151 (MQ=255)                                                                                                                                    ttCGGCGTATCGCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAAcc              <  1:616636/149‑1 (MQ=255)                                                                                                                                    tcgggCTTACCGCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAAcc              <  2:63583/146‑1 (MQ=255)                                                                                                                                      cGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGcc           >  1:483144/1‑150 (MQ=255)                                                                                                                                       ggCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCa          >  2:256774/1‑150 (MQ=255)                                                                                                                                         cGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGcc           >  2:205879/1‑147 (MQ=255)                                                                                                                                          gTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGAc      >  1:962971/1‑151 (MQ=255)                                                                                                                                              cGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGtt  >  2:832683/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAA‑CT‑GCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTT  >  NC_000913/364384‑364673

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