Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,286,478	0	A	C	41.8%	87.8 / 10.3	56	G298G (GGT→GGG)	acs	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (23/9);  new base C (1/22);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 2.90e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTG  >  NC_000913/4286335‑4286627                                                                                                                                                |                                                                                                                                                       accaccTGCGCCATACGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGgat                                                                                                                                                 >  2:426701/1‑151 (MQ=255)  ccaccTGCGCCATACGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGgat                                                                                                                                                 <  2:184391/150‑1 (MQ=255)     ccTGCGCCATACGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATa                                                                                                                                             <  2:226421/151‑1 (MQ=255)        gcgcCATACGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATaa                                                                                                                                            <  1:525448/149‑1 (MQ=255)        gcgcCATACGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATc                                                                                                                                          >  2:480072/1‑151 (MQ=255)                cGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACat                                                                                                                                  <  1:425466/151‑1 (MQ=255)                   gCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATAtt                                                                                                                               >  1:664777/1‑151 (MQ=255)                     aGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTa                                                                                                                             <  1:458247/151‑1 (MQ=255)                         gTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAg                                                                                                                          >  1:502782/1‑150 (MQ=255)                               cAGTTGGGTACGCCTTCAAACATCAGCGTGGGCGCACCGCAGTCCAGCGGGCCGTCCAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCCCCAGTAGATATTACCCAGATGATAATCAAAGACATATTTAAAGGTCa                                                                                                                      >  1:538323/1‑148 (MQ=255)                                aGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑gcgc                                                                                                                 >  2:732171/1‑151 (MQ=255)                                 gTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑Gcgcc                                                                                                                <  2:700515/151‑1 (MQ=255)                                   tGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑Gcgccgc                                                                                                              >  2:649715/1‑151 (MQ=255)                                   tGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑Gcgccg                                                                                                               >  2:153478/1‑150 (MQ=255)                                               cAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATa                                                                                                  <  2:642470/151‑1 (MQ=255)                                                aaaCATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATa                                                                                                  >  2:766015/1‑150 (MQ=255)                                                      cAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAAccgcc                                                                                            <  2:359961/150‑1 (MQ=255)                                                                  aCCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATCTCCCCCGGATGATAATCAAAGACATATTTAAAGGCCA‑GCGCCGCGTACAACAGATAACCGCCGGTAGTATgcag                                                                                >  1:125024/1‑150 (MQ=255)                                                                   ccGCGGCCCAGGGGCCCTTACAGCAAGTACCTGTGCCCGGCCCCCCACCCCACAGCGGCGGTCCACCAGGAGATACCCCCCGGATGAAAATCAAAGACATATTAAAAGGTCA‑GCGCCGCTTCCCCCAGATAACCGCCGGTAGTAt                                                                                    <  1:783832/145‑1 (MQ=255)                                                                     gCAGCCCCGCGGACCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACCTCGGCGGTGGTCCAGGAGTTCTCTCCCACCTGGTTAACCAAGACATATTTCGAGGACG‑GCGCCGGCGCCACCAGGCAACCGCCGGTGCGATgcagca                                                                              >  1:728126/1‑149 (MQ=255)                                                                        gcccggggggCCGAACAGCAGGTAGTTGTGGCGGGTCCCCCAGCCCACACCGGGGGTGCACGAGTAGATACCCCCGGGATGAAAACCAAAGACATATTAAAAGGTCA‑GCGCCGCGTCCACCAGATAACCGCCGGTAGTATgcag                                                                                <  1:631793/137‑1 (MQ=255)                                                                          ccAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGGAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGGACACCAGATAACCGCCGGTAGGAAGCAGCACAACtt                                                                        >  1:650161/1‑150 (MQ=255)                                                                            aGCGGGCCGTACAGCAAGACACTGTGTCCGGTCACCCAGCCCACATCGGCGGCGCACCAGTAAATATCAGCCGGATTATAAACAAAAAAATAATTAAAGGACA‑GCGCCGTGTAAAACAGATAACCGGCGGTAGGATgc                                                                                  >  1:663323/1‑138 (MQ=255)                                                                               gggCCGTCCAGCAGGTATCTGTGCCCGGTCACCCCCCGCACCTGGGCGGTCCCCCGGCAAATACCCCCCGGAGGATAACCAAGGACATTTTAAAACGTCG‑GCGCGCCGTCCGCCAGCTAACCGCCGGTAGTATGCAGCACTCCt                                                                         <  1:342884/144‑1 (MQ=255)                                                                                 gCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGACATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCtt                                                                 >  2:223234/1‑150 (MQ=255)                                                                                  ccGTACAGCAAGTAACTGTGTCCGGTCCCCCAGCCCACATCGGCGGTGCACCAGTAGACACCCCCCGGATGAAAATCAAAGACATATTAAAGGGTCA‑GCGCCGCGTACCCCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTa                                                                <  1:223212/150‑1 (MQ=255)                                                                                   cGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTAcc                                                              >  1:694378/1‑151 (MQ=255)                                                                                     tACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCgg                                                            >  2:674078/1‑151 (MQ=255)                                                                                     tACAGCAAGTAACTGTGCCCGGTCCCCCAGCCCACATCGGGGGTGCACCAGTAGATACCACCCGGATAAAAATCAAAGAAATATTTAAAGTTCA‑GCGCGGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTa                                                                <  1:156704/147‑1 (MQ=255)                                                                                           aaGTAACTGTGTCCGGTCCCCCAGCCCACACCGGCGGTGCACCAGTAGATATCCCCCGGATAATAATCAAAGACATATTAAAAGGTCA‑GCGCCGCGTACCCAAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTa                                                          <  2:718236/147‑1 (MQ=255)                                                                                           aaGTAACTGGTCCCGTTCACCCAGCCCACACCGGCGGTGCCCCAGTAGATATCCCCCGGATAATAATCAAAGACATTTTTAAAGGTCG‑GCGCCGCGTCCCCAAAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGa                                                        <  1:695428/149‑1 (MQ=255)                                                                                                aCTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCgg                                                   >  2:757267/1‑149 (MQ=255)                                                                                                aCTGGGCCCGGTCCCCCACCCCACACCGGGGTGGCACCAGTAAATACCCCCCGGATGAAAATCAAAGAAATATTAAAAGGCCG‑GCGCCGCGTACCCAAAAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAg                                                 <  1:566139/151‑1 (MQ=255)                                                                                                  tgtgTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAgg                                                >  1:221784/1‑150 (MQ=255)                                                                                                    tgtCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGATATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGgtgt                                             >  2:766879/1‑151 (MQ=255)                                                                                                    tgtCCGCCCACCCACCCCACACCGGCGGTGCACCAGTAGAACCCCCCCGAATGATAAACAAAGAAATA‑TTAAAGGGCAGGCGCCGCCTACCCAAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGgt                                               <  1:104054/149‑1 (MQ=255)                                                                                                       ccGGTCCCCCACCCCACACCGGCGGTGCCCCAGTAGATACCCCCCGGATGATAATAAAAGACATATTAAAAGGCCG‑GCGCCGCGTACACAAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTag                                           <  1:5394/150‑1 (MQ=255)                                                                                                       ccGGTCCCCAAGCCCACAGCGGCGGTGCCCCAGTAGATACCCCCCGGATAATAATAAAAGACATATTAAAGGGCCG‑GCGCGGCGTACCCAAAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTaga                                          <  1:567150/151‑1 (MQ=255)                                                                                                            cccccAGCCCACATCGGCGGTGCACCAGTAGATATCCCCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATa                                     <  2:617868/149‑1 (MQ=255)                                                                                                            cacccaccCCACATCGGCGGTGCACCAGTAGATATCCCCCGTATGATAATCAAAGACATATTAAAAGGTCA‑GCGCCGCGTACCCAAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGaa                                       <  1:92345/149‑1 (MQ=255)                                                                                                              cccaccccacAGCGGCGGTGCACCATTAGATACCCCCGGTATAATAATCAAAGACATATTAAAAGGTCA‑GCGCCGCTTACACCAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATaaa                                   <  2:702320/151‑1 (MQ=255)                                                                                                              cccAGCCCCCATGGGGGGGGCCCCAGTAGATACCCCCCGAATGAAATCCAAAGACATATTAAAAGGTCA‑GCGCCGCCTACACCAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATaaa                                   <  2:139898/151‑1 (MQ=255)                                                                                                                cggCCCACACGGGGGGTGAACAAGTAGAAATCCCCCGGAGGAAAACAAAAGAAATATTAAAAGGTCG‑GCGCCGCGTACACCAAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACa                                 <  1:295924/149‑1 (MQ=255)                                                                                                                cAGCCCCCAGGGGCGGTGCCCCAGTAGATACCCCCCGGAAGATAACCAAAGACATTTTTAAGGGCCA‑GCGCCGCGTACACCAGATACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACa                                 <  1:531812/151‑1 (MQ=255)                                                                                                                   cccACATCGGCGGTGCACCAGTAGATATCCCCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACCCCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAg                                <  2:66832/149‑1 (MQ=255)                                                                                                                   cccACACCGGCGGTGCACCAGTAGATATCCCCCGGATGAAAATCAAAGACATATTTAAAGGCCA‑GCGCCGGGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAg                                <  2:188333/149‑1 (MQ=255)                                                                                                                     cacaTCGGCGGTCCACCAGTAGATATCACCCGGATGTTATTCACAGACCTATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCCCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGGAGAGAATAAACAGCgg                             >  1:198651/1‑150 (MQ=255)                                                                                                                        accggcgtTCCCCCAGTAGAAACCCCCCGGATGAAAATAAAAAACATATTAAAAGGCCA‑GCGCGGCTTCCACAAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGAt                           <  1:425085/147‑1 (MQ=255)                                                                                                                           gccggTCCACCATAAAATATCCCCCGGATGATAACCAAAGCAATATTAAAAGGTCA‑GCGCCGCTTCCACAAAAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCtt                        <  1:245692/147‑1 (MQ=255)                                                                                                                            gcggTCCACCAGTAGATATCCCCCGGATAATAATCAAAGACATATTAAAGGGTCA‑GCGCCGCGTCCACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCgg                     <  2:138026/151‑1 (MQ=255)                                                                                                                               gTGCCCCAATAGATATCCCCCGAATGATAAAAAAAGACATATTAAAAGGTCG‑GCGCCGCGTACCCCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGc                    <  1:490764/149‑1 (MQ=255)                                                                                                                                    ccAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGAGTAGAGAATAAACAGCGGATATTCGGCGTTCAtc             >  1:785706/1‑151 (MQ=255)                                                                                                                                      aGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCAtctc           >  2:589983/1‑151 (MQ=255)                                                                                                                                         aGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTc        >  1:114276/1‑151 (MQ=255)                                                                                                                                              tCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCt   >  1:434086/1‑151 (MQ=255)                                                                                                                                               cACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTg  >  2:740204/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                       ACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCA‑GCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTG  >  NC_000913/4286335‑4286627

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