Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,286,504	0	T	A	44.4%	71.9 / 30.0	54	K290* (AAA→TAA)	acs	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base T (20/10);  new base A (0/24);  total (20/34)
Fisher's exact test for biased strand distribution p-value = 9.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACC  >  NC_000913/4286354‑4286652                                                                                                                                                        |                                                                                                                                                     cAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATAttt                                                                                                                                                      <  1:292753/151‑1 (MQ=255)         ggCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCa                                                                                                                                              <  1:982565/151‑1 (MQ=255)             aGTTGGGTACGCCTTCAAACATCAGCGTGGGCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAgcgc                                                                                                                                          <  1:674638/151‑1 (MQ=255)              gTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGcgcc                                                                                                                                         <  2:911678/151‑1 (MQ=255)                 gggTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGcgccgc                                                                                                                                       <  1:183608/150‑1 (MQ=255)                 gggTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGcgcc                                                                                                                                         >  2:77244/1‑148 (MQ=255)                     aCGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTaca                                                                                                                                  <  1:1022386/151‑1 (MQ=255)                     aCGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTac                                                                                                                                   <  2:13885/150‑1 (MQ=255)                     aCGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTac                                                                                                                                   <  2:374051/150‑1 (MQ=255)                          ttCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACAcc                                                                                                                                >  1:993880/1‑148 (MQ=255)                            cAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGGAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGAt                                                                                                                            >  2:786184/1‑150 (MQ=255)                               aCATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCAACAGTAGATATTACCCGG‑ATGATAATCAAAGACATATTTAAAGGATCGCGCCGCGTACAACAGATAAcc                                                                                                                        >  2:1009314/1‑151 (MQ=255)                                    aGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCgg                                                                                                                   >  1:1070082/1‑151 (MQ=255)                                     gCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCgg                                                                                                                   <  2:779319/150‑1 (MQ=255)                                         ggTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGggagt                                                                                                               >  2:821344/1‑150 (MQ=255)                                            cGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTATGTCCGGGCACCCAGCCCACATCGGCGGCGCACCAGGAGGAATCACCCGG‑ATGATAATCAAAGACATAATTTAAGGGCCGCGCCGCGTAAACCAGATAACCGCCGGGAGAATgc                                                                                                           >  1:950052/1‑151 (MQ=255)                                             gCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATgc                                                                                                           >  1:163069/1‑150 (MQ=255)                                                ccGCGGGCCAGCGGGCCGTACAAAAAGTAACGGTGCCCGGCCACCCACCCCCCATCGGCGGTTCACCAGAAGATAGCACCCGG‑AAAATAAACAAAGACATATTAAAAGGTCACCGCCGCGTACCCCAGATAACCGCCGGTAGTATgcag                                                                                                         <  2:459089/149‑1 (MQ=255)                                                ccGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATgcagc                                                                                                        >  1:199289/1‑150 (MQ=255)                                                        cAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCtttt                                                                                               >  1:582582/1‑151 (MQ=255)                                                          gggggCCTTACAGAAAGTAACTGTGTCCGTCCCCCCACCCCACAGCGGCGGTGCACAAGTAGAATCCCCCCGG‑AGAATAAAAAAAGAAATATTAAAAGGGCGGCCCCGCGTACCCCAGAAAACCGCCGGTAGTATGCAGCACACCTTTTg                                                                                              <  1:374685/148‑1 (MQ=255)                                                          gCGGGCCGTACAGCAATTAACTGTTCCCGGTCACCCACCCCACATCGGCGGTGCACCAGTAAATACCCCCCGG‑TTGATAATCAAAGACATATTAAAAGGCCGGCGCCGCGTACCCCAAAAAACCGCCGGTAGTATGCAGCACACCtt                                                                                                 <  1:886863/147‑1 (MQ=255)                                                          gCGGGCCGAACAGCAAGTAACTGTGTCCGGTCCCCCACCCCACATCGGCGGTGCACCAGTAGATATCCCCGGG‑ATGATAATCAAAGAAATATTTAAAGGCCGGCGCCCCGTCCCCAAGATAACCGCCGGTAGTATGCAGCACACCtttt                                                                                               <  1:825397/149‑1 (MQ=255)                                                                  tACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCgg                                                                                     >  1:541454/1‑151 (MQ=255)                                                                     aGCAAGTAACTGTGTCCGGTCCCCCACCCCACATCGGCGGTGCACCAGTAGATATCCCCCGG‑ATAATAATAAAAGACATATTAAAAGGTCAGCGCCGCGTACACCAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCgg                                                                                     <  1:8134/148‑1 (MQ=255)                                                                      gaaggTAACTGCGCCGGCTCCCCCACCCCACAGCGGCGGTGCACCGGAAAATCCCCCCCAA‑ATAATACACAAAAAAATTTTAAAGGGATTGGGCCACGTCCACGAAATACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAg                                                                                  <  1:121323/146‑1 (MQ=255)                                                                        aattaccTGTGTCCGGCCCCCCACCCCCCAGCGGCGGTCCACCAGTAAATACCACCCGG‑ATAATAATAAAAGACATATTAAAAGGCCGGCGCCGCGTACCCCAAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAAc                                                                               <  2:795036/145‑1 (MQ=255)                                                                        aaGTAACTGGGTCCGGTCCCCCAGCCCACACGGGCGGTGCCCCAGTAGATACCCCCGGG‑ATAAAAAAAAAATACTTATTTAAGGGCCAGCCCCGCGTACAACAAAAACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAAc                                                                               <  1:1112964/151‑1 (MQ=255)                                                                        aaGTAACTGGGTCCGGTCACCCAGCCCACAGCGGGGGTCCCCCAGTAGATACCCCCCGG‑ATAATAATAAAAGACATATTAAAGGGCCAGCCCCGCTTACACCAAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAg                                                                                  <  1:790075/148‑1 (MQ=255)                                                                        aaGAAACTGTGCCCGGGCCCCCCGCCCACATGGGGGGCTCACCATTAGATCCCCCCGGG‑GAGAAAACAAAAGACATTTTAAAGGGGCGGCGCCGCGCCCCCCAAATACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAAc                                                                               <  2:844290/151‑1 (MQ=255)                                                                           tAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAAATAACCGCCGGTAGTATGCAGCACAACCTTTGGCTTACCGGTAGAACCgg                                                                            >  2:303169/1‑151 (MQ=255)                                                                               tgtgTCCGGTCCCCCACCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATAAAAGACATATTAAAAGGTCAGCGCCGCGTACCCCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAg                                                                          <  2:931153/149‑1 (MQ=255)                                                                                    ccgtccccccaccccacACCGGCGGTGCACCAGTAGATACCCCCCGG‑ATAATAATCAAAGACATATTAAAAGGCCACCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGgtgt                                                                      <  2:793792/141‑1 (MQ=255)                                                                                        ccccccaccccacAGCGGCGGTGCCCCAGATGATACCCCCGGG‑ATAATAAAAAAAGACATTTTAAAAGGTCAGCGCCGCGTACACCAGATACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAAt                                                               <  1:660998/148‑1 (MQ=255)                                                                                           cccgcccaacaTGGGGGGTGCACCAG‑GGATTACCCCCGGAATAATAAACAAAAAAATATTAAAAGGTCCGCGCCGCGTCTACGAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAAt                                                               <  2:434640/140‑1 (MQ=255)                                                                                            ccACCCCACACCGGCGGTGCACAAGTAGATACCCCCCGG‑AAGAAAATAAAAGACATATTAAAAGGCCGGCGCCGCGTACACCAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAAc                                                           <  2:363264/151‑1 (MQ=255)                                                                                              aGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAAc                                                           >  1:850417/1‑149 (MQ=255)                                                                                               ccccACCTCGGGGGTCCCCCAGTAGATACCCCCCGA‑ATAAAAACAAAAGACATATTAAAAGCTCACCCCCGCTAACCCCAAAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGc                                                        <  2:925388/150‑1 (MQ=255)                                                                                                    ccgcggcggTCCCCAAGAAAATACCCCCCGG‑ATGAAAAAAAAAAAAATATTAAAGGGCCAGCGCCGCTTACACAAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATc                                                   <  1:22783/148‑1 (MQ=255)                                                                                                     agcgggggTGCACCAGTAGATATCCCCCGG‑ATGATAACAAAAGAAATATTAAAAGGTCAGCGCCGCGTCCACCAAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGAt                                                    <  1:714809/147‑1 (MQ=255)                                                                                                        ggcggTGCACCAGTAGATATCCCCCGG‑ATGAAAATAAAAGACATATTAAAAGGTCAGCGCCGCGTACCCCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCt                                                  <  1:634244/148‑1 (MQ=255)                                                                                                           ggTGCCCCAGTAGATACCCCCGGT‑ATAAAAAAAAAAGCCATATTAAAGGGTCAGCGCCGCGTCCCCCAGATATCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCgg                                              <  1:446410/149‑1 (MQ=255)                                                                                                               ccagaGTAGATATCCCCCGG‑ATGAAAATAAAAGAAATTTTAAAAGGCCAGCGCCGCGTACCCCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTc                                         <  1:533190/146‑1 (MQ=255)                                                                                                               cACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTc                                         >  1:215750/1‑150 (MQ=255)                                                                                                                caaagaaGATACCCCCGGG‑ATAAAAACAAAAAACATATTAAAAGGTCGGCCCCGCGAACCCAAAATACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGtt                                          <  2:878146/142‑1 (MQ=255)                                                                                                                        atatCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCg                               >  1:1065133/1‑151 (MQ=255)                                                                                                                        ataccccccGG‑AAAATAATCAAACACATATTAAAGGGTCGGCGCCGCTTACCCCAAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCg                               <  1:49855/145‑1 (MQ=255)                                                                                                                            cccccGG‑ATGATAATCAAAGACATATTAAAGGGCCAGCGCCGCGTACACCAAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTcc                                <  1:74941/144‑1 (MQ=255)                                                                                                                              cccGG‑ATAATAAACAAAAAAATATTAAAAGGCCAGCCCCGCTTACCACAAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGt                         <  2:910019/151‑1 (MQ=255)                                                                                                                                cGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGc                       >  2:22758/1‑151 (MQ=255)                                                                                                                                        tAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGAtcg                 >  2:522585/1‑150 (MQ=255)                                                                                                                                         aataaaaacattatTAAAAGTCCAGCGCCGCTTCCCCCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGAtcgct               <  1:1077341/140‑1 (MQ=255)                                                                                                                                               aGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTg         >  2:539803/1‑151 (MQ=255)                                                                                                                                                 aCATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCt       >  2:137449/1‑151 (MQ=255)                                                                                                                                                 aCATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTCCCGCCTGGTGCTGATCGCTCGCTTGCt       >  2:579007/1‑151 (MQ=255)                                                                                                                                                       ttAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAAcc  >  2:627913/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                     CAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGG‑ATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACC  >  NC_000913/4286354‑4286652

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