Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	43.9%	86.5 / 17.1	57	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/4);  new base C (0/25);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 2.72e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCG  >  NC_000913/2984293‑2984587                                                                                                                                                     |                                                                                                                                                   aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                  <  2:129372/151‑1 (MQ=255)  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                                 >  1:845631/1‑151 (MQ=255)    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt                                                                                                                                                <  1:619186/150‑1 (MQ=255)    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                               >  2:605635/1‑151 (MQ=255)        cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccc                                                                                                                                           <  1:972111/151‑1 (MQ=255)         aTGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcccc                                                                                                                                          >  1:930874/1‑151 (MQ=255)                  gAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATAc                                                                                                                                 >  2:770642/1‑151 (MQ=255)                       aTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                                                                            >  2:536750/1‑151 (MQ=255)                               gCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                                                                      >  1:197747/1‑149 (MQ=255)                               gCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAg                                                                                                                    <  2:1131146/151‑1 (MQ=255)                                   aaaaaaTGTATGGCGACTTGCCACCCAAAATTTTTCAACAAATTAAGGCCGGAAGATATTTATGGTTAAGGAGGAAAGGCAGACCTTTTCGTCACGTAAAATGGTAAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                <  1:1148872/149‑1 (MQ=255)                                    aaaaaTGTATGTCGCTTTCCCCCCCATATTTATTAAAAAAATTAAGGCCGGAATATTTTTATCGGAAAGGAGAAAATGCAGAGTTTTTCGTCCCTTTCATTGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                               <  1:965399/150‑1 (MQ=255)                                     aaaaTGTATGTCGTCTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATTTTTATGGGTAAGGAGAAAAGGCAGAGTTTTTCGTCACTTTCAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                               <  2:903643/150‑1 (MQ=255)                                          gtatgtCGACTTCCCACCCATCATTATTAAACAAATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGGGTTTTTCGTCACTTTCAAGGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                           <  2:123162/149‑1 (MQ=255)                                           tatgtCGACTTGCCAGCCATCTTTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                          <  2:703006/149‑1 (MQ=255)                                                 gACTTGCCAGCTATCATTATTCAACAATTTAATAGCGGCAAATATTTACCGGAAAGGAGAAAAGGCAGAGCTTTTCGTCACGTTAAATGGCCAACGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                                      <  2:65827/147‑1 (MQ=255)                                                    ttccaagccaTCTTTATTAAAAAAATAAATGACGGAAGATATTTATGGGTAGGGAGAAAAGGCAGAGCTTTTCGTCCGTTTAAATGGTAATCGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                                                                 <  1:531013/144‑1 (MQ=255)                                                         agccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                           >  1:131543/1‑150 (MQ=255)                                                           ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                        >  1:165595/1‑151 (MQ=255)                                                             atcatTATCAAACAAATTAATGCCGGCAGATATTTATCGGAAAGGAGAAAATGAAGACTTTTTCGCCACGTTCAATGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                        <  2:138452/149‑1 (MQ=255)                                                                       aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                               >  1:1055999/1‑148 (MQ=255)                                                                       aaCAAATAAATGACGGCAGATTTTTATCGGTAGGGAGAAAATGCAGAGTTTTTCCTCACTTTCAATGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                             <  2:900651/150‑1 (MQ=255)                                                                        aCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                             >  1:314153/1‑149 (MQ=255)                                                                          aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                          >  1:800435/1‑150 (MQ=255)                                                                            ataaaTGACGGAAGATATTTATGGGTAAGGAGAAAATGCAGAGCTATTCGTCATTTTAAATGGTAAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa                                                                        <  2:64407/147‑1 (MQ=255)                                                                              tAATGACGGCAAATATTTATCGGTAAGGAGAAAAGGCAGAGCTTTTCGTCCCGTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                       <  1:606577/149‑1 (MQ=255)                                                                               aaggACGGCAGATTTTTATGGGAAAGGAGAAAATGCAGAGTTTTTCCTCACTTTAAATGGTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                     <  1:979040/147‑1 (MQ=255)                                                                                aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                   >  2:338049/1‑151 (MQ=255)                                                                                 tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                   >  1:964983/1‑150 (MQ=255)                                                                                    cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                               >  1:492447/1‑151 (MQ=255)                                                                                    cGGAAGATATTTATCGGTAGGGAGAAAATGCAGACCTTTTCGTCACGTTCAAGGGTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                                <  2:854089/150‑1 (MQ=255)                                                                                     ggCAGATATTTATCGGTAAGGAGAAAATGCAGACCTTTTCGTCACGTTAAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                  <  2:998685/147‑1 (MQ=255)                                                                                        aGATATTTATCGGCAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                            >  1:787872/1‑150 (MQ=255)                                                                                                atgggTAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAAGGGTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa                                                    <  2:148537/147‑1 (MQ=255)                                                                                                  cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                 >  2:132416/1‑151 (MQ=255)                                                                                                   ggTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                                >  1:662767/1‑151 (MQ=255)                                                                                                    gTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                               >  1:462979/1‑151 (MQ=255)                                                                                                       aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                            >  2:207414/1‑151 (MQ=255)                                                                                                       aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACAGGAAGCGCCTa                                             >  2:108236/1‑150 (MQ=255)                                                                                                        ggagaaAATGCAGAGCTTTTCGCCCCTTTCAATGGTCACTCCCCCCCCCGCCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                               <  1:620931/147‑1 (MQ=255)                                                                                                         gaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                          >  2:191195/1‑151 (MQ=255)                                                                                                          agaaaaTGAAGCGCTTTTCTTCACTTTCAATGGACAAGGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                      <  1:214161/135‑1 (MQ=255)                                                                                                           gaaaaTGCAGATCTTTTCGTCCCGTTCAATGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                         <  1:78170/148‑1 (MQ=255)                                                                                                             aaaTGCAGAGCTTTTCGTCACGTTAAATGGTAAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                         <  2:238402/148‑1 (MQ=255)                                                                                                                 gCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                  >  1:652270/1‑151 (MQ=255)                                                                                                                   agagTTTTTCGTCCCGTTAAATGGAAAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                   <  2:110302/148‑1 (MQ=255)                                                                                                                        ttttcttcACGTTCAATGGGCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                             <  1:1087738/147‑1 (MQ=255)                                                                                                                          ttCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                         >  1:801680/1‑151 (MQ=255)                                                                                                                           tCGTCCCGTTAAAGGGTAAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                         <  2:988545/150‑1 (MQ=255)                                                                                                                             gtcagtttaaaTGTTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                        <  1:912741/140‑1 (MQ=255)                                                                                                                              tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                     >  2:118482/1‑151 (MQ=255)                                                                                                                               cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGAGACCGAGTGCGATcgcg                      >  2:892653/1‑149 (MQ=255)                                                                                                                                     cAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt              >  2:674339/1‑151 (MQ=255)                                                                                                                                      aaTGTTAAAGGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat               <  2:444892/149‑1 (MQ=255)                                                                                                                                       aTGGTCACTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg             >  1:304282/1‑150 (MQ=255)                                                                                                                                         ggTCAACCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat               <  2:426510/146‑1 (MQ=255)                                                                                                                                                 ccACACCCCCGCCCCCGCCGATCCAAAGCCTCCCCAGCCCTCTGCGGGCATCACGCTTTACCCTTCCACGAACCCGTGACACCAGGATTCGGCAACCGGAAGTGCCTACCGGGTGCCCGAGCGCGCTCGCGCCACCATTCACACTGACTCg  >  1:908980/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCG  >  NC_000913/2984293‑2984587

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