Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	370,876	0	C	G	40.4%	81.7 / 22.8	47	S200R (AGC→AGG)	mhpB	2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (24/4);  new base G (0/19);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.27e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

ACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCAT  >  NC_000913/370727‑371025                                                                                                                                                          |                                                                                                                                                        aCCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCg                                                                                                                                                        >  1:535928/1‑151 (MQ=255)       aTGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAg                                                                                                                                                  <  1:324398/151‑1 (MQ=255)                       ttGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTg                                                                                                                                 >  2:422527/1‑150 (MQ=255)                       ttGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTg                                                                                                                                 >  2:75655/1‑150 (MQ=255)                                 cACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAAtt                                                                                                                       >  2:147297/1‑150 (MQ=255)                                  aCCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGc                                                                                                                     >  1:3076/1‑151 (MQ=255)                                    ccgCACTCTAAAAAACCGCGTCTTTTTCCTGTGTCCGGGTGGTTTTC‑CCCCCCGGCCGCGGGTCCCGGAAG‑GGGAGAAA‑GCCGATCCCA‑ATATCGGCGCCGGTTTG‑TGGGGGGGGGGAAAGTT‑TTTGCCCGCCAGTGAGCGCGAATTGCGt                                                                                                                   <  1:131329/149‑1 (MQ=255)                                       ccctctAAAAAACCGGGTGTTG‑TCCCGGGTTCCGGGGGGGTTTCCCCACCAGCCGCGGGGGCCCGAAC‑GGGAAAAA‑GCCGATGCCA‑ATATGCGCGCCCGTTGGTGGGGGGGGGGAAAAG‑T‑TTTACCCGCCAGTGAGCGCGAATTGCGt                                                                                                                   <  2:922465/146‑1 (MQ=255)                                              aaTAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACggg                                                                                                         >  2:737329/1‑151 (MQ=255)                                              aaTAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACggg                                                                                                         >  1:341013/1‑151 (MQ=255)                                               aaaaacgggggtttttCCTGGGGTCCGGGGGGGTTTCCCCACCGGCCGCGGGGGCCCGAGC‑GGGAAAAG‑GCCGGTCCCA‑ATGGGCGCGCCGGTGTGTGGGGGGGGGGAAAAG‑T‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCa                                                                                                              <  1:805471/132‑1 (MQ=255)                                                 aacggcgTGTTGTTCGTGGGTTCGGGGGGTCTTC‑CCCCTCCGCGGCGGGTGCCCGAGC‑GGGAGAAA‑GCCGCTGCCC‑ATATGCGCGCCCGTTTGTGGGGGGGGGGAAAAT‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGa                                                                                                      <  2:292094/147‑1 (MQ=255)                                                 aaaCGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTg                                                                                                       >  1:401884/1‑150 (MQ=255)                                                      cgTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCTAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAgc                                                                                                 >  1:282033/1‑151 (MQ=255)                                                         gCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgcc                                                                                              >  2:650537/1‑151 (MQ=255)                                                           tGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgccgc                                                                                            >  2:808729/1‑151 (MQ=255)                                                            tttCCTGTGTGCGGGGGGGCTTC‑CCCCTCGGCCGCGGGTGCCCGAGC‑TGGCGAAA‑GCGGTTCCAA‑ATTGGCGCGACCGTCTGTTGGGGGGCGGGAAGG‑T‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgcc                                                                                              <  2:300759/147‑1 (MQ=255)                                                                  ggttcccGGGGGGCTTC‑CCCATCGGCCGCGGGCGCCCGAAC‑GGGCAAAG‑CCCGAGCCCC‑GGGTGCGCGCCCGTCTGTTGGGGGGGGGGAAAA‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAg                                                                                     <  2:376122/146‑1 (MQ=255)                                                                  gggTTGGGGTGGTCTTC‑CCCATCCCCCGCGGGGCCCCGACC‑GGGAAAAA‑GCGGAGCCCC‑ATTGGCGCGCCGGTTTGTTGGGGGGGGGAAAAG‑T‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgccg                                                                                             <  2:884622/143‑1 (MQ=255)                                                                  gggTTCCGGTGGTTTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGCCCGTCTGTGGGGGGGGGGAAAAG‑A‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAg                                                                                     <  2:599628/151‑1 (MQ=255)                                                                  gggTTCCGGGGGGCCTT‑CCCCCCCGCCCGCGGGGCCCGGACGGGGCAAAA‑GCGGACGCCA‑ATAGGCGCGACGGTTTGTGGGGGGGGGGGAAAT‑T‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAtt                                                                                                    <  1:95495/137‑1 (MQ=255)                                                                   ggTTGCGGGGGGTTTT‑CCCAACCGCCCGCGGGGGCCGAAC‑TGGGGGAA‑GGCGGTTCCCAATGTCCCCGCCGGTGTGTGGGGGGGGGGAAAAT‑T‑TTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTgaga                                                                                       <  1:56173/149‑1 (MQ=255)                                                                    gTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGtt                                                                                   >  2:300892/1‑151 (MQ=255)                                                                      tCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttg                                                                                 >  2:412794/1‑151 (MQ=255)                                                                       ccGGGGGGTTTC‑CCCACCGGCCGCGGGTGCCCGAAC‑GGGAAAAA‑GCGGATCCCC‑ATGTGCGCGCCCGTCTGGGGGGGGGGGGAAAAG‑T‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAg                                                                                     <  2:852463/146‑1 (MQ=255)                                                                        gggggggCTTC‑CCCACCGGCCGCGGGTCCCCGAGC‑GGGAAAAA‑GCGGAGGCCA‑ATAGGCGCGCCGGTGTGTGGGGGGGGGGAAAAG‑A‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                               <  2:365493/147‑1 (MQ=255)                                                                        cGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                               >  1:89207/1‑151 (MQ=255)                                                                         ggTGTGTTTT‑CCCTCCCGCCGCCGTTGCCCGAAC‑TGGAGAAG‑GACGATCCCC‑ATACGCGCGACTGTCTGTGGGGGAGGGGGAAAT‑A‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGtt                                                                                   <  1:872436/146‑1 (MQ=255)                                                                         ggTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgttg                                                                              >  2:493067/1‑151 (MQ=255)                                                                         ggTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgttg                                                                              >  1:862661/1‑151 (MQ=255)                                                                                     ccaccccccgcggGTGCCGAATC‑TGGCAAAG‑GCCGATCCCA‑ATGGGCGGGCCTGTCTGTGGGGGGGGGGGAAAG‑T‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGa                                                                          <  1:965436/137‑1 (MQ=255)                                                                                     ccATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCag                                                                      >  2:405704/1‑148 (MQ=255)                                                                                      ctgcccccgccgGCG‑CCGAAT‑CGGGGAAAGGCCGATGCGC‑AAGGGCGCGCCTGTTTGTGGGGGAGGGGAAAAT‑T‑TTTCGCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAgaga                                                                   <  2:755996/144‑1 (MQ=255)                                                                                          aGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCt                                                              >  1:679386/1‑151 (MQ=255)                                                                                                  gTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGc                                                      >  2:102066/1‑151 (MQ=255)                                                                                                        gAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAcc                                                >  1:847538/1‑151 (MQ=255)                                                                                                                 gAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTgg                                        >  1:336840/1‑151 (MQ=255)                                                                                                                  aaa‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTggg                                       >  1:544703/1‑151 (MQ=255)                                                                                                                             ccc‑ATGTGGGCGACCGTGTGGGGGGGGGGGGAAAAG‑T‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAg                               <  1:249907/149‑1 (MQ=255)                                                                                                                                 atatGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCAt                          >  1:269682/1‑151 (MQ=255)                                                                                                                                  tatGCGCGCCCGTTTGGGGGGGAGCGGAAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCACTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCa                           <  1:301429/149‑1 (MQ=255)                                                                                                                                   atGCGCGCCCGTCTGTGGGGGGGGGGGAAAG‑T‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCa                           <  2:882746/148‑1 (MQ=255)                                                                                                                                        cgACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTg                   >  1:650083/1‑151 (MQ=255)                                                                                                                                             cgttggtGGGGGGGGGGAAAA‑ATTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACt                      <  2:948689/138‑1 (MQ=255)                                                                                                                                                tgtgGGGGGGGGGGAAAG‑A‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGc            <  1:112010/150‑1 (MQ=255)                                                                                                                                                tgtgGGGGCGGGGAAAAT‑T‑TTCACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGa              <  1:459319/148‑1 (MQ=255)                                                                                                                                                   gggggggggGAAAAG‑T‑TTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGa              <  2:368053/144‑1 (MQ=255)                                                                                                                                                         gCGGAAAAT‑T‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCa   <  2:108743/150‑1 (MQ=255)                                                                                                                                                          cGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCAt  >  2:665459/1‑150 (MQ=255)                                                                                                                                                          |                                                                                                                                                        ACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCATCAGCCGCCGGTGCCCGAAC‑TGGCGAAA‑GCCGATGCCC‑ATATGCGCGACCGTCTGTTGGGGAGCGGGAAAG‑A‑TTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCAT  >  NC_000913/370727‑371025

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