Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336,901	0	T	A	41.7%	64.3 / 21.2	49	Y44F (TAT→TTT)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (19/9);  new base A (0/20);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 6.36e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCG  >  NC_000913/4336751‑4337046                                                                                                                                                       |                                                                                                                                                  cagTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAccatccat                                                                                                                                                   <  2:135550/151‑1 (MQ=255) cagTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAccatccat                                                                                                                                                   >  1:69194/1‑151 (MQ=255)  agTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAccatccat                                                                                                                                                   >  2:172342/1‑150 (MQ=255)     agagGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATaaa                                                                                                                                                <  1:120208/150‑1 (MQ=255)           cGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAAt                                                                                                                                         <  1:626543/151‑1 (MQ=255)              tGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATcc                                                                                                                                       <  1:801075/150‑1 (MQ=255)              tGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATcc                                                                                                                                       <  1:644938/150‑1 (MQ=255)                              aaaCGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAgg                                                                                                                      >  1:259382/1‑151 (MQ=255)                              aaaCGGGCAAAGGAAGGGGCGGGCGAAAGCGTAAGAACCCCCAGGACTGGGGCCGAGGAACGACATTTGGGCGTATCCCATCGAGAGCCCCAGTGCCCCGAAAATCGCCACCAACCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCa                                                                                                                        <  2:506991/149‑1 (MQ=255)                                  gggCCAAAGCGGCGGCGGGCGTAAGCTTAAGAACCCCCAGGACTTGGGCCGAGAAACGACTTTTTGGCTTATCCCATCGAGGGCCCAAGTGCCCCGATAATCGTCACCACCCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGttt                                                                                                                   <  1:1002448/150‑1 (MQ=255)                                   ggCAAAAGCCGGGGCGGGAGTAAGCATAAGAACCACCAGGGCTGGGGCGGAGAACCGCCATTTGGGCGTATACCACCGAGGCCCCCAGTGCCCCAAAAATCGCCACCACCCACCCAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTg                                                                                                                  <  2:1167628/150‑1 (MQ=255)                                       aaaGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGCCATTTGGGCTTATACCATCGAGAGCCCCAGGGCCCCGATAATCGTCCCCAACCACCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAgg                                                                                                              <  2:1079069/150‑1 (MQ=255)                                         agcagcGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTa                                                                                                            >  2:1103655/1‑150 (MQ=255)                                            agcGGCGGGCGTAAGCGTAAGAACCACCAGGACTTTGGTCGGGGGACGACATTTTTGCGTATACCATCGAGAGCCCCAGGGCACCGATAATCGTCACCAAACATCCATAAATGGCAATCCCGCCCGTAGAGGCCAGGGTTGCAGGGAACa                                                                                                         >  2:658465/1‑150 (MQ=255)                                               ggcgggcgTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaa                                                                                                      >  1:325134/1‑150 (MQ=255)                                                 cgggcgTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaaa                                                                                                    >  2:1172185/1‑150 (MQ=255)                                                   ggcgTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAac                                                                                                  >  1:458035/1‑150 (MQ=255)                                                   ggcgTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAac                                                                                                  >  2:1167230/1‑150 (MQ=255)                                                        aagcgaaaaAACCCCCAGGCCTTGGGTCGAGGAACGACTTTTTGGCGTATCCCATCGAGACCCCCAGTCCCCCGATAACCGCCCCCACCCATCCAAAAAGGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACAcc                                                                                               <  2:323341/148‑1 (MQ=255)                                                           ataaaaacccaccaGGTTTGGGGTCGAGGAACGACATTTTGGCATATACCACCGGGAGCCCCAGTGCCCGGATAATCGCCACAACCCATCCAAAAATGGAAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAcc                                                                                          <  2:634410/142‑1 (MQ=255)                                                              aagAACCACCAGGACTGGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGACCCCCAGTGCACCAATAATCGCCACCAACCATCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAcccc                                                                                        <  1:823710/149‑1 (MQ=255)                                                               aaacccaccaGGACTTGGGTCGAGGACCGAATTTTTGGCGTATCCCATCGAGAGCCCCAGTGCACCGATAATCGCCACCAACCACCAAAAAATGACAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCata                                                                                     <  2:360481/147‑1 (MQ=255)                                                                   ccaccaGGACTGGGGGCGAGGAACGACATTTTGGCGTATCCCATCGAGACCCCCGGTGCCCCGAAAACCGTCACCACCCATCCAAAAATGCCACTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataata                                                                                  <  2:867010/150‑1 (MQ=255)                                                                         ggACTGGGGTCGAGGAACGACATTTTGGCGTATCCCATCGAGAGCCCCAGGGCACCGATAATCGCCCCCAACCACCCAAAAATGACAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAACATTCCCCg                                                                           <  2:471448/151‑1 (MQ=255)                                                                          gACTTGGGTCGAGAAACGACTTTTTGGCGTATCCCACCGAGGGCCCCAGTGCACCAATAATCGTCACCACCCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGa                                                                          <  1:985644/151‑1 (MQ=255)                                                                             tggggTCGAGGAACATTATTTTGGCATATCCCACCGAGCGCCCAAGTGCCCCAATAATCGTCCCAAACCATCCAGAAGTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac                                                                         <  1:602837/147‑1 (MQ=255)                                                                               gggTCGAGGAACGACATTTTGGCGTATCCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCACCCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac                                                                         <  1:1132692/147‑1 (MQ=255)                                                                                  tCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATc                                                                   >  2:926424/1‑150 (MQ=255)                                                                                    gAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAgg                                                                >  1:66676/1‑151 (MQ=255)                                                                                    gAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGCCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAgg                                                                >  2:1012075/1‑151 (MQ=255)                                                                                    gAGGAACGACATTTGGGCGTATCCCATCGAGAGCCCCAGTGCACCGATAATCGCCCCCAACCACCAAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCa                                                                  <  2:322529/149‑1 (MQ=255)                                                                                     aGGAACGACTTTTGGGCGTATCCCATCGAGACCCCCAGTGCACCGATAACCGTCCCCACCCACCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAgg                                                                <  1:678042/150‑1 (MQ=255)                                                                                       gAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTg                                                             >  2:110684/1‑151 (MQ=255)                                                                                        aaCGACTTTTTGGCGTATCCCATGGAGACCCCCAGTGCCCGGATAACCGCCCCCAACCACCAAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCa                                                                  <  2:1144264/145‑1 (MQ=255)                                                                                                     cGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGcc                                               >  1:1016651/1‑151 (MQ=255)                                                                                                     cGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGcc                                               >  2:155929/1‑151 (MQ=255)                                                                                                       tataCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGccc                                              >  2:587700/1‑150 (MQ=255)                                                                                                        ataCCACCGAGAGCCCCGGGGCCCCGATAATCGCCACCACCCACCCAAAAATGGAAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCa                                             <  2:204487/150‑1 (MQ=255)                                                                                                             accaagagCCCCAGTGCACCGATAATCGTCCCCAACCACCCAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgt                                       <  2:569015/147‑1 (MQ=255)                                                                                                             accGAGACCCCCAGGCCCCCAATAACCGCCCCCACCCACCAAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgt                                       <  2:803803/149‑1 (MQ=255)                                                                                                               cGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                      >  2:689130/1‑150 (MQ=255)                                                                                                                          gTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagcatc                          >  2:1177292/1‑151 (MQ=255)                                                                                                                             cccccaTACTCGCCACCACCCCCCAAAAAATGGAAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGa                       <  2:92012/146‑1 (MQ=255)                                                                                                                             cACCGATAATCGCCCCAAACCATCAAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                        <  2:1039507/150‑1 (MQ=255)                                                                                                                             cACCGAAAATCTCCCCCACCCACCCAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                        <  1:112125/150‑1 (MQ=255)                                                                                                                                 gaaaaTCGTCACCACCCATCAAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagcat                           <  2:544165/140‑1 (MQ=255)                                                                                                                                      tCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAt              >  2:806187/1‑151 (MQ=255)                                                                                                                                       cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtt             >  2:28955/1‑151 (MQ=255)                                                                                                                                          ccccacccaTCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa            <  2:63558/147‑1 (MQ=255)                                                                                                                                          cACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaa           <  1:753033/150‑1 (MQ=255)                                                                                                                                            cccaccaccCAAAAAGGCAAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACt        <  2:552393/148‑1 (MQ=255)                                                                                                                                             caaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc       >  2:692647/1‑151 (MQ=255)                                                                                                                                              acccatccaaAAACGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaaa          <  2:416881/145‑1 (MQ=255)                                                                                                                                              acccatcaaaAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACt        <  2:374665/147‑1 (MQ=255)                                                                                                                                                ccatccaaAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCt     <  1:341329/150‑1 (MQ=255)                                                                                                                                                ccacccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCt     <  1:212096/150‑1 (MQ=255)                                                                                                                                                ccacccaaAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCt     <  1:1030621/150‑1 (MQ=255)                                                                                                                                                   tccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCg  >  1:1206727/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                  CAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCG  >  NC_000913/4336751‑4337046

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