Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,326,450	0	T	G	47.9%	64.6 / 17.3	48	G114G (GGT→GGG)	atoB	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (3/22);  new base G (23/0);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 8.48e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAA  >  NC_000913/2326301‑2326597                                                                                                                                                      |                                                                                                                                                    ccGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGt                                                                                                                                                     <  2:666053/150‑1 (MQ=255) ccGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGt                                                                                                                                                     <  1:296092/150‑1 (MQ=255) ccGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTa                                                                                                                                                    <  2:872241/151‑1 (MQ=255)          cAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAata                                                                                                                                           <  2:863685/151‑1 (MQ=255)              cACTGTTAAAAAGCGGGCTGGCAGAAACGGTTTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGa                                                                                                                                        <  1:195010/150‑1 (MQ=255)                cTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGGAAAAATGAGtt                                                                                                                                     >  1:622580/1‑151 (MQ=255)                  gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                                     <  1:936122/149‑1 (MQ=255)                  gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                                                   <  2:92867/151‑1 (MQ=255)                         aGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGTTTAGccccc                                                                                                                             >  2:548122/1‑150 (MQ=255)                         aGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGGTTAGccccc                                                                                                                             >  1:632053/1‑150 (MQ=255)                                 ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttac                                                                                                                      <  1:605701/149‑1 (MQ=255)                                    aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGAt                                                                                                                 <  1:813810/151‑1 (MQ=255)                                             gtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGCTTAGCCCCCCCCTTTCTCGATGCcaaa                                                                                                           >  1:229705/1‑148 (MQ=255)                                                 ggATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCCACTTACTCCATGCAAAAGCACg                                                                                                      >  1:229158/1‑149 (MQ=255)                                                    ttCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGccc                                                                                                   >  1:106185/1‑147 (MQ=255)                                                     tCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTgg                                                                                                >  2:356001/1‑151 (MQ=255)                                                      cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGTTGGGAAAAATGAGTTTAGCCCCCCACTTTCTCGACGCAAAAGAAAGCTCTGGt                                                                                               >  1:121995/1‑151 (MQ=255)                                                                aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCtt                                                                                      <  1:202738/150‑1 (MQ=255)                                                                aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCtt                                                                                      <  1:732035/150‑1 (MQ=255)                                                                aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCCACCTACTCGGTGCAAAAACACGCCCTGGGTATTGTCtt                                                                                      >  1:744015/1‑150 (MQ=255)                                                                 aaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTgg                                                                                    <  1:892353/151‑1 (MQ=255)                                                                 aaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGTTTAGCCCCCTACCTACTCGATGGCAAAACACGCCCTGGTTATCGTCTTgg                                                                                    >  1:118258/1‑151 (MQ=255)                                                                   gTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATTGCAAACATGCACTTTTCCCCCTACTTACTCGTTCTGAAAACACACTCTGGCTACTGTCt                                                                                       >  1:627738/1‑146 (MQ=255)                                                                   gTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTgg                                                                                    >  2:187842/1‑149 (MQ=255)                                                                       gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACg                                                                               <  1:70734/150‑1 (MQ=255)                                                                           ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACgg                                                                              >  2:348149/1‑147 (MQ=255)                                                                                gTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGGCAGGTTTa                                                                     <  2:620433/151‑1 (MQ=255)                                                                                  cTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTACTTAATCGATGCAAAAGCACGCTCTGGTTATTGTCTTTGGGAAGGACAGGgtt                                                                      >  2:611315/1‑148 (MQ=255)                                                                                   ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGGGACGGACAGGGTTAt                                                                    >  1:462099/1‑149 (MQ=255)                                                                                    tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGGACGCTCTGGGTATCGGCTTGGGGACGGACAGGgtt                                                                      >  2:135929/1‑146 (MQ=255)                                                                                     aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGGAATATGAGTTTAGCCCCCTTCTTACTCGATGGAAAAGCACGCTCTGGTTATCGGCTTGGAGACGGGCAGGTTTATGAc                                                                 >  1:896177/1‑150 (MQ=255)                                                                                       aaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCCACTTACTCGATGCAAAAGCACGCCCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGt                                                               >  1:164543/1‑150 (MQ=255)                                                                                        aGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTaa                                                             <  1:61157/151‑1 (MQ=255)                                                                                             ggCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCt                                                         <  2:452643/150‑1 (MQ=255)                                                                                              gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgc                                                       <  2:586258/151‑1 (MQ=255)                                                                                                 cTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGACCCCCACTTTCTCGATGCAAAAGCACGCTCTGGGTTTCGTCTTGGGGGCGGACAGGGTTTTGACGGAAAcc                                                          >  1:533148/1‑145 (MQ=255)                                                                                                       gccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATTAGTTTAGCCCCCTACTTAATCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGGTgg                                                >  1:8156/1‑149 (MQ=255)                                                                                                           aGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATTACGTAATACTTCGCGATGGCCTGAt                                          >  2:654146/1‑151 (MQ=255)                                                                                                            ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGAt                                          <  2:699624/150‑1 (MQ=255)                                                                                                                aTTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTgcg                                     <  2:649018/151‑1 (MQ=255)                                                                                                                 ttCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATTAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTgcg                                     >  2:424756/1‑150 (MQ=255)                                                                                                                    aggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCCACTTACTCGGTGCAAAAGCACGCCCTGGTTATTGGCTTGGGGACGGACAGGGTTATGACGTAATCCTGCGCGATGGGCTGATGGGGGCCAc                                 >  1:695565/1‑151 (MQ=255)                                                                                                                        aggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCag                             <  1:803447/151‑2 (MQ=255)                                                                                                                               gcgcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATc                      <  2:799160/151‑1 (MQ=255)                                                                                                                                cgcAGAGCATTGTGGCGGGGGGGATGGAAAATAAGAGTTTAGCCCCCTACTTATTCGAGGCAAAAGCACGCTCTGGGTCATGTCTTGGGGACGGAAAGGTTTATTACGTAATTCTTCGCGGTGGCCCGACGGGCGCCACCCATGGttttt                      >  1:535706/1‑147 (MQ=255)                                                                                                                                    gagCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATg                 <  1:37043/151‑1 (MQ=255)                                                                                                                                              gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccg       <  2:467379/151‑1 (MQ=255)                                                                                                                                               cGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgc      >  1:849243/1‑151 (MQ=255)                                                                                                                                               cGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgc      >  1:533476/1‑151 (MQ=255)                                                                                                                                                   gggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaa  <  1:187838/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAA  >  NC_000913/2326301‑2326597

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