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

GGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTT  >  NC_000913/364376‑364673                                                                                                                                                       |                                                                                                                                                    ggTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGcc                                                                                                                                                     >  2:587831/1‑151 (MQ=255)       cGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAAt                                                                                                                                                <  1:834413/150‑1 (MQ=255)        gCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATc                                                                                                                                               >  1:830404/1‑150 (MQ=255)                       ggTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGcc                                                                                                                                  <  1:68183/148‑1 (MQ=255)                          ttGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTTTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGACCTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCa                                                                                                                             <  2:69669/150‑1 (MQ=255)                               cGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGtt                                                                                                                       >  2:928813/1‑151 (MQ=255)                                  aTAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGcc                                                                                                                    >  1:851353/1‑151 (MQ=255)                                   tAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCg                                                                                                                   >  1:504268/1‑151 (MQ=255)                                        ggAACTGAAAAACCTGCGGCGGGTGTTTTGCTTCCGTCGGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACGGGCGATCGTTCGGCGTACCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGtttt                                                                                                               <  1:716652/150‑1 (MQ=255)                                        gaaaCTGAAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGGGTGCGGTCGGCAAAGACCAGACCGTTCATCCAGAACGGGGGATCGTTGGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGttt                                                                                                                <  1:851351/147‑1 (MQ=255)                                                aaaaggtgctgttgGTGTTTTGCTGCCTACAGCGCTGGGTGCGGGGGGCGGTCGGCAAAGACCACACCGTCCATACAAACCTGGGGATCGTTCGGCGTTTCGCAAAAAGCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTcatcat                                                                                                         <  1:406360/142‑1 (MQ=255)                                                 aaaaCTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAttt                                                                                                     >  1:426831/1‑151 (MQ=255)                                                       gctgctgGTTTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTGGGCAAAGCCCAGCCCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAAt                                                                                                  <  2:882468/148‑1 (MQ=255)                                                       gctgcggGGTTTTTGCTCCCGTCAGCTGGGGATGCGGCGTGCGGGCGGGAAAGCCCAGCCCGTTCATACAGACCTGGAGCTCGTCCGCGTTGTCCCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATc                                                                                                 <  2:778674/149‑1 (MQ=255)                                                        ctgctgTTGTTTTGCTCCCGTCAGCGCGGGATGCGTCGTGCGGCCGGAAAAGACCAGACCTTTCATCAAGACCGGGCGATCGTTCGGCGTATCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAg                                                                                               <  2:601153/150‑1 (MQ=255)                                                         tgctgGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCg                                                                                             >  1:175804/1‑151 (MQ=255)                                                          gctgGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCg                                                                                             <  1:248133/150‑1 (MQ=255)                                                             ggTTTTTTGCTTCCGTCGGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAAt                                                                                                  <  1:522591/142‑1 (MQ=255)                                                                    tGCTTCCGTCAGCGCTGGATGGGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAc                                                                                   <  1:684531/150‑1 (MQ=255)                                                                           gTCAGCGCTGGATGGGGCGTGCGGTCGGCAAAGCCCAGACCGTTCATACAGAACGGGCGATCGTTCGGCGTATCCCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTccc                                                                           <  2:176847/151‑1 (MQ=255)                                                                           gTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTccc                                                                           >  2:19894/1‑151 (MQ=255)                                                                              aGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTccc                                                                           >  2:96084/1‑148 (MQ=255)                                                                               gcgcTGGATGCGGCGTGCGGCGGGAAAAGCCCAGCCCGTTCATACAGACCGGGCGATCGTTGGGCGTTTCCCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGa                                                                        <  2:447101/150‑1 (MQ=255)                                                                                   tGGTTGCGGCGGGCGGTCGGCAAAGACCAGCCCGTTCATACAGACCTGGCGACCGTTCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACg                                                                      <  2:449428/148‑1 (MQ=255)                                                                                          ggcgTGCGGGCGGAAAAGCCCAGACCGTTCATACAGAACGGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCt                                                            <  2:538981/151‑1 (MQ=255)                                                                                          ggcgTGCGGCCGGCAAAGCCCAGCCCGTTCATCCAGACCTGGCGACCGTTCGGCGTACCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAgccg                                                                <  2:667054/147‑1 (MQ=255)                                                                                           gcgTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTg                                                           >  2:801604/1‑151 (MQ=255)                                                                                             gTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTCTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGt                                                          >  1:95512/1‑150 (MQ=255)                                                                                                    cGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                                  >  2:719042/1‑151 (MQ=255)                                                                                                     ggCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                                  >  2:818005/1‑150 (MQ=255)                                                                                                     ggCAAAGACCAAACCTTTAATACAAAACGGGCGATCGTCCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                                  <  1:872410/150‑1 (MQ=255)                                                                                                        aaaGCCCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGAt                                                <  1:326688/149‑1 (MQ=255)                                                                                                        aaaGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGAt                                                <  1:57266/149‑1 (MQ=255)                                                                                                               agaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGa                                        >  2:459189/1‑150 (MQ=255)                                                                                                                  ccGTTCATCCAGAACTGGGGATCGTTCGGCGTATCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACg                                    <  1:474602/151‑1 (MQ=255)                                                                                                                          aCAGATCTGGCGATCGTGCGGCGTACCGCCAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAg                            <  1:756769/151‑1 (MQ=255)                                                                                                                          aCAGAACTGGCGATCGTTCGGGGTACCGCAAAAATCTCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCa                             <  1:214176/150‑1 (MQ=255)                                                                                                                           cAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGt                           >  2:205028/1‑151 (MQ=255)                                                                                                                            aGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGt                           >  1:767296/1‑150 (MQ=255)                                                                                                                                cTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAttt                       >  2:708388/1‑150 (MQ=255)                                                                                                                                 gggcgttgttggggCGTGCCGCAAAAACCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTa                      <  2:23270/138‑1 (MQ=255)                                                                                                                                  ggtgatcgatcGGCTTATCGCAAAAATCGCCGCCGTAAGCCGACCGCGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGACCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTa                      <  1:223466/146‑1 (MQ=255)                                                                                                                                    cGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCg                   >  2:126780/1‑150 (MQ=255)                                                                                                                                     gaccGTGCGGGGTGTCGCAAAACTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGa                  <  2:763036/147‑1 (MQ=255)                                                                                                                                         gTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCg             >  1:208261/1‑151 (MQ=255)                                                                                                                                             ggCGTATCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCa          <  1:849544/150‑1 (MQ=255)                                                                                                                                              gCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCa          >  2:493606/1‑149 (MQ=255)                                                                                                                                                    cGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGtt  >  2:431711/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                    GGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTT  >  NC_000913/364376‑364673

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