Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	172,911	0	A	G	46.5%	56.4 / 13.7	43	M484V (ATG→GTG)	fhuB	iron(3+)‑hydroxamate import ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/23);  new base G (20/0);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.04e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGGTGTTGATGCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATT  >  NC_000913/172762‑173060                                                                                                                                                       |                                                                                                                                                        gCGTGGTGTTGATGCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGAt                                                                                                                                                        <  1:801618/151‑1 (MQ=255)                tGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGG‑G‑GTTAGCACcgc                                                                                                                                       >  1:188530/1‑151 (MQ=255)                 gTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACcg                                                                                                                                        <  1:125963/149‑1 (MQ=255)                        gTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGGTTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGG‑G‑GTTAGCAcccccgtgcacc                                                                                                                               >  2:372758/1‑141 (MQ=255)                            cGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCt                                                                                                                           <  1:183453/151‑1 (MQ=255)                            cGGGTAATGCCTTTGGCTGACTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGTCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCt                                                                                                                           <  1:852867/151‑1 (MQ=255)                               gTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCtttt                                                                                                                        <  2:169198/151‑1 (MQ=255)                                    gCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTtgatga                                                                                                                   <  1:910941/151‑1 (MQ=255)                                    gCCTGTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTtgatg                                                                                                                    <  2:10799/150‑1 (MQ=255)                                      cTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGG‑G‑TTGAGCACCCCCGTCCCCCCGCTTTtggtgtt                                                                                                                  >  2:110894/1‑148 (MQ=255)                                                  gTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAgt                                                                                                     <  1:411059/151‑1 (MQ=255)                                                       cTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACTGTATGTTACTGGC‑GGGGGGGGG‑GTTTAAGC‑CCCCGGTCCCCCTGCTTTTTTTGATGTTggcgggaaatgg                                                                                                   >  1:116304/1‑137 (MQ=255)                                                           aggcagTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGGGGGGG‑G‑TTAAG‑ACCCCGGTCACCCTGCTTTTTTTGGTGTTTCAGGCAAgggggggaccc                                                                                            >  1:83382/1‑144 (MQ=255)                                                             gcagTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGG‑G‑GTAAGCAACGCGGTTACCATGCTTTTTTTGATGTTTTAGGCAAgggggg                                                                                                 >  1:446192/1‑142 (MQ=255)                                                               agTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGG‑G‑GTTAGCACCCCGCTCACCCATCTTTTTGTGATTTTTCCAGGAAGTGGGGACCcccg                                                                                          >  2:331226/1‑149 (MQ=255)                                                               agTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGC‑G‑GTTAACACCCCGTTCCCCCTGCTTTTGTTGATTTTTCAGGCCAGAGGGGGCCcgc                                                                                           >  1:128941/1‑148 (MQ=255)                                                                    cggcgcggcgGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATg                                                                                      <  2:891479/148‑1 (MQ=255)                                                                         cggcggTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGt                                                                              <  1:905948/151‑1 (MQ=255)                                                                           gcggTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGG‑G‑GTAAGCACCCCGTTCACCATGCTTTTTATGATGTTGCAGGCAAGAGGTGACCCGCGAATGGGGCAAGTg                                                                             >  1:409210/1‑150 (MQ=255)                                                                           gcggTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGc                                                                            <  1:731252/151‑1 (MQ=255)                                                                             ggTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGC‑G‑GTTAGCACCGCGGTCACCATGCTTTTTGTGATGTTGCAGGCAAGTGGGGACCCCCGAAAGGGGCAAGTg                                                                             >  1:882825/1‑148 (MQ=255)                                                                                   gCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTgg                                                                    <  2:503687/151‑1 (MQ=255)                                                                                      gTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGAtt                                                                 <  1:813257/151‑1 (MQ=255)                                                                                        tGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGC‑G‑TTAAGCACCGCGCTCCCCCTGCTTTTGTTGATGTTGCAGGCAAGTGGGGACCCGCGAATGGCGCAAGGGCTGACCTGGAttt                                                                >  1:870867/1‑150 (MQ=255)                                                                                         gATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCTCCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTcc                                                              <  2:916447/151‑1 (MQ=255)                                                                                          agcATTATGGTCGCCGCCGGCCGCGGTGGATTTTCCCCACGCCGTATGTTACTGGC‑GGGGAGGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTGGATGATGTTGCAGGCAAGTGGTGACCGGCGAATGGCGCAAGTGCTGACCTGGATTTc                                                               <  2:408724/147‑1 (MQ=255)                                                                                            cATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGC‑G‑TTAAGCACCGCCTTCCCCCTGCTTTTTATTATGTTGCAGGCAAGGGGTGGCCCGCGGATGGCGCAAGTGCCGACCCGGGttttc                                                              >  2:860542/1‑146 (MQ=255)                                                                                             aTTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGG‑G‑TTAAGCACCCCGTTCACCCTGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCCCGAAAGGCGCAAGAGCTGGCCTGGGTTTCCgg                                                            >  1:402933/1‑149 (MQ=255)                                                                                                 tGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGCGG‑TTAAG‑ACCGCGGTCACCATGCTTTTTGTGATGTTGCAGGCAAGTGGGGACCCGCGAAAGGGGCAAGTGCTGACCTGGGTTTTCGGTTCg                                                        >  1:141734/1‑149 (MQ=255)                                                                                                   aTCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGGGGG‑G‑GTTAACAACGCGGTCACCATGCTTTTTATGGTGTTGGAGGCCAGTGGTGAACCCCGAAAGGCGCAAGGGCTGGCCTGGGTTTTCGGGTCGAcc                                                     >  1:317729/1‑150 (MQ=255)                                                                                                         gccggccgCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAAcg                                              <  1:308571/151‑1 (MQ=255)                                                                                                          ccggccgCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGGGGG‑G‑GTTAACAACCCGCTCACCCTCCTTTTTTTGTTGTTGGAGGGAAGTGGTGGACCCCGCAAGGGGCCAGTGGTGGTCTGGGTTTTCgggctgacctccca                                                >  2:603097/1‑136 (MQ=255)                                                                                                          ccggccgCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGGGGC‑G‑GTTAGCACCGCGTTTACCATGCTTTTTGTGATGTTTCAGGCAAGGGGTGACCCGCGGATGGCGCAAGTGCTGACCCGGGTTTTCGGGTCGACCTACaa                                                >  1:940446/1‑148 (MQ=255)                                                                                                               cgcgGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGC‑G‑TTTAACACCGCGGTCACCATGCTTTTTTTGATGTTGCAGGCAAGTGGTGACCCGCGAAAGGCGCAAGTGCTGACCTGGGTTTCCGGTTCGACCCACAAAGCGAcc                                         >  1:964569/1‑150 (MQ=255)                                                                                                                  ggTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGC‑G‑GTAAGCACCCCGTTCACCATGCTTTTTGTGATGTTGCAGGCAAGTGGGGACCCGCGAATGGGGCCAGTGCTGACCTGGGTTTTCGGTTCGACCTACAACGCGACCGATg                                     >  1:621877/1‑151 (MQ=255)                                                                                                                   gTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATgc                                    <  2:99825/151‑1 (MQ=255)                                                                                                                   gTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATgc                                    <  2:649559/151‑1 (MQ=255)                                                                                                                       aTTTTCCCCACACCGTATGTTACTGGC‑GGGGGTGGG‑G‑TTAAACACCGCGGTCACCATGCTTTTTTTGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTTCGGTTCGACCTACAACGCGACCGATgc                                    >  1:231386/1‑147 (MQ=255)                                                                                                                         tttCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAgg                               <  1:755408/150‑1 (MQ=255)                                                                                                                               cacaCCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCGAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGgcg                        <  2:377546/151‑1 (MQ=255)                                                                                                                                         gTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAAtt              <  1:471987/151‑1 (MQ=255)                                                                                                                                                   ggggATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTg     <  1:77245/151‑1 (MQ=255)                                                                                                                                                      gATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGAtt  <  1:397362/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                        GCGTGGTGTTGATGCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGC‑GGGGATGGC‑G‑TTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATT  >  NC_000913/172762‑173060

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