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

TATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAA  >  NC_000913/364378‑364659                                                                                                                                                     |                                                                                                                                      tATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCa                                                                                                                                      >  1:222889/1‑150 (MQ=255)      gCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATc                                                                                                                                 >  2:401189/1‑150 (MQ=255)         ggTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCTccg                                                                                                                             <  2:962951/151‑1 (MQ=255)         ggTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCAc                                                                                                                               <  1:371055/149‑1 (MQ=255)                tCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTaa                                                                                                                       >  1:433238/1‑150 (MQ=255)                tCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTa                                                                                                                        >  1:18242/1‑149 (MQ=255)                   aTGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCg                                                                                                                   >  1:23838/1‑151 (MQ=255)                     gttttGCCCGGATAAACGGAACGGAAAAACCTGCTGCTGGTTTTTTGCTCCCGTCAGCGCTGGAGGGGGCGTGCGGTCGGAAAAGACCAGACCGTTCATACAGAACGGGCGATCGTTCGGGGTACCGCAAAAATCCCCGCCGTaa                                                                                                                       <  1:214152/143‑1 (MQ=255)                              gaaaaaaCGGAACGGAAAAACCTGCTGCGGTTGTTTTTCTCCCGTCAGCGCGGGATGGGGCGTGCGGCCGGAAAAGACCAGCCCTTTCATACAGACCGGGCGATCGTGCGGCGTACCCCAAAAATCACCGCCGTAAGCCGACCACGGGtt                                                                                                         <  1:428827/146‑1 (MQ=255)                                  aaggggAACGGAAAAACCTTCTGCTGGTTTTTTGCTCCCGTCAGCGCTGGATGCGGCGTGCGGGCGGAAAAGCCCAGCCCGTTCATACAGAACGGGCGATCGTTCGGCTTACCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGt                                                                                                    <  1:489652/147‑1 (MQ=255)                                  aaaCGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGt                                                                                                    >  2:992193/1‑151 (MQ=255)                                         aCTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtcatc                                                                                             >  2:165813/1‑151 (MQ=255)                                          cTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtcatc                                                                                             >  1:545441/1‑150 (MQ=255)                                          cTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtcatc                                                                                             >  2:296785/1‑150 (MQ=255)                                            ggAAAAACTGCTGCTGGTGTTTTGCTCCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCTTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTcatcat                                                                                           <  2:667364/150‑1 (MQ=255)                                              aaaaaCTGCTGCTGTTTTTTTGCTCCCGTCAGCGCGGGAGGCGGCGTGCGGTCGGCAAAGACCAGACCTTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAtt                                                                                        <  2:786979/151‑1 (MQ=255)                                               aaatCTGCTGGTGTTTTTTTGCTTCCGTCAGCGCGGGGTGCGGCGGGCGGACGGCAAAGACCAGCCTGTTCATAAAGACCGGGCGATTGTTCGGCTTCTCGCCAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAttt                                                                                       <  2:357407/151‑1 (MQ=255)                                               aaacCTGCTGCTGGTGTTTTGCTTCCGTCAGCGCGGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATAAAGAACGGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAttt                                                                                       <  2:806646/151‑1 (MQ=255)                                                aaaCTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAtt                                                                                        >  2:437204/1‑149 (MQ=255)                                                      ctgctgGTGTTTTCCTCCCGTCAGGGCTGGTTGGGGCGTGCGGCCGGAAAAGACCACACCGTCCATACAGAACGGGCGATCGTCCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCa                                                                                  <  2:243792/149‑1 (MQ=255)                                                         ctgGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGAc                                                                             >  1:859342/1‑151 (MQ=255)                                                          tgGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACt                                                                            >  1:486206/1‑151 (MQ=255)                                                          tgGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACt                                                                            >  2:359378/1‑151 (MQ=255)                                                           ggTTTTTTGCTTCCGGCAGCCCTGGATGTGGCGTGCGGTCGGCAACGCCCAGACCTTTCACACAGACGTGGCGATCGTTGGGGGTACCGCAAAAATCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTg                                                                           <  2:251611/151‑1 (MQ=255)                                                                tttGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATcc                                                                       >  2:36213/1‑150 (MQ=255)                                                                  tGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCTACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAcc                                                                    >  1:699267/1‑151 (MQ=255)                                                                  tGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAcc                                                                    >  2:418992/1‑151 (MQ=255)                                                                       ccGTCAGCGCTGGATGCGGGGTGCGGTCGGCAAAGACCAGACCTTTCATACAGACCGGGCGATCGTTCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCa                                                                  <  1:712552/148‑1 (MQ=255)                                                                       ccGCCAGCGCTGGATGCGGCTTGCGGCCGGCAAAGCCCAGCCCTTTAAACAAGAACTGGCGATCGTTCGGTTTACCGCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAg                                                                 <  2:96760/149‑1 (MQ=255)                                                                           cAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTccc                                                             >  1:696833/1‑149 (MQ=255)                                                                            accgcTGGAGGCGGCGGGCGGTCGGCAAAGCCCAGCCCTTTCATAAAGACCTGGCGATGGTTGGGCGTTCCGCAAAAACCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGa                                                          <  1:3418/149‑1 (MQ=255)                                                                             gcggggaaTGCGGCTTGCGGTCGACAAACCCCAAACCTTGCATACAAAAGGGGCGATCTTGCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTcc                                                              <  1:990287/139‑1 (MQ=255)                                                                              cgcTGGATGCGGCGTGCGGCCGGCAAAGACCAGACCGTTCATACACAACTGGCGATCGTTCGCCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGGTGCCGTTCTCTTCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                         >  1:420047/1‑150 (MQ=255)                                                                               gcTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACg                                                        >  2:214064/1‑150 (MQ=255)                                                                                  ggATGCGGCGTGCGGTCGGAAAAGCCCAGCCCGTTCATACAAAACTGGCGATCGTCCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAgc                                                    <  2:68343/151‑1 (MQ=255)                                                                                            tGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTa                                           >  1:874226/1‑150 (MQ=255)                                                                                               ggTGGGAAAAGACCGGACCGTTCATACAGAACGGGCGATCGTCGGGCGTTCCCCCAAATTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACg                                       <  1:550159/151‑1 (MQ=255)                                                                                                   ggCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                    >  2:331133/1‑150 (MQ=255)                                                                                                     aaaaGACCAGCCCTTACATAAAAAACTGGCGATCGTTCGGCTTACCGCAAAATTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATa                                 <  1:383023/150‑1 (MQ=255)                                                                                                      aaaGACCAGCCCGTTCATACAGACCTGGCGATGGTTCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATa                                 <  1:414140/150‑1 (MQ=255)                                                                                                           caagaccTTTCATACAGAACTGGCGTTCGTTCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTg                              <  1:533060/146‑1 (MQ=255)                                                                                                                ccGTTCATACAGACCGGGCGATCGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGaaa                        <  1:165747/149‑1 (MQ=255)                                                                                                                      aTACAAAACGGGCGATCGTTGGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGc                 <  2:18250/150‑1 (MQ=255)                                                                                                                            aaCTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAtt          >  2:691513/1‑151 (MQ=255)                                                                                                                             aCTGGCGATATTGCGTCTTATCCCAATAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAttt         <  1:741612/151‑1 (MQ=255)                                                                                                                               tGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAttt         >  1:900297/1‑149 (MQ=255)                                                                                                                                agCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGc      >  2:443406/2‑151 (MQ=255)                                                                                                                                 gCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCg     >  2:118944/1‑151 (MQ=255)                                                                                                                                    aTCGTTCGCGGTACCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATGCTGACGAAACGCCTGCCAGTATTTAGCg     <  2:180273/148‑1 (MQ=255)                                                                                                                                       gTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGaaa  <  2:428459/148‑1 (MQ=255)                                                                                                                                           ggCGTATCGCAAAAATCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTa        <  2:113633/138‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                      TATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAA  >  NC_000913/364378‑364659

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