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

AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCG  >  NC_000913/2984293‑2984587                                                                                                                                                     |                                                                                                                                                   aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                  >  2:36603/1‑151 (MQ=255) aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                  >  1:844230/1‑151 (MQ=255)  atatGTCATGAGAATGGAAAAATGAGAGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                                 >  2:47564/1‑151 (MQ=255)     tGTCATGAGAATGGAAAAATGATGGAGCCCATAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGTAGGAAATGCAGAGCTATTCGTCACGTTCAATAGTCAATGCCACAccc                                                                                                                                                 >  2:381973/1‑148 (MQ=255)           gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc                                                                                                                                         <  1:921861/150‑1 (MQ=255)                    aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGACAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACa                                                                                                                                >  2:842494/1‑150 (MQ=255)                           gggAGCCCAAAAAATGTTTTTCGATTTGCCAGCTTTCATTATACAACAGATTAATGGGGGAAGATATTTACCGGTAGGGAGGAAATGCAGAGTTTTTCGTCACGTTCAATGGTTAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTg                                                                                                                        <  1:769883/151‑1 (MQ=255)                               gCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                                                                      >  1:24317/1‑149 (MQ=255)                                  cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                 >  1:885582/1‑151 (MQ=255)                                        atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                           >  1:926227/1‑151 (MQ=255)                                         tgtttgcCGACTTGCCACCCATAATTATCAAAAAGATAAATAACGGCATATTTTTCTCGGTAAGGAGAAAATGCAGGGCTTTTCGTCCCGTTCAATGCTCAACCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                            <  2:664977/149‑1 (MQ=255)                                               tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                     >  1:794640/1‑150 (MQ=255)                                                       ccagcaatatttttccAACAGTTTAAGGCCGGAAGATATTTATCGGTAAGAAGAAAATCCAGAGTTTTTCGTCACTTTCATTGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                             <  2:783089/135‑1 (MQ=255)                                                        cacccaTCATTATTAAACAAATTAATGCCGGCAGATATTTAGGGGAAAGGAGGAAATGCAGGCTTTTTCGTCACTTTAAATGGTAAAGGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                     <  1:673580/138‑1 (MQ=255)                                                         agccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                           >  1:611894/1‑150 (MQ=255)                                                               catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                    >  2:176279/1‑151 (MQ=255)                                                                 ttattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                   >  2:755565/1‑150 (MQ=255)                                                                 ttattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                   >  1:769884/1‑150 (MQ=255)                                                                   attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGAAACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                >  2:518694/1‑151 (MQ=255)                                                                   attAAACATATAAATGACGGCAGATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACTTTCAATGGTAAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                  <  1:870618/149‑1 (MQ=255)                                                                      cAACAGATTAGTGGCGGCAGATATTTTTCGGAAGGGAGAAAATCCAGAGCTATTCGTCCCTTACAATGTTAAACGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                               <  1:532728/149‑1 (MQ=255)                                                                      cAACAGATTAATGACGGCAGATATTTATAGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCAACACCCTGGCACCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTGTCACCATTTCATGAAc                                                                                >  1:678208/1‑148 (MQ=255)                                                                        aCATATTAATGCCGGAAGATATTTATGGGAAAGGAGAAAAGGCAGAGCTTTTCGTCACTTCCAATGGTAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                           <  2:76308/151‑1 (MQ=255)                                                                         cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                             >  1:918752/1‑148 (MQ=255)                                                                          aGATTACTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGCGCTATTCGGCACGTTCATTGGTCAATGCCTAACCCTGATCCCCGCCGATACAAAGCGTCGCTGGTCCTTTGCGGGCATTACGTTTCACCATTTCTTGAACCAgaga                                                                          >  2:633027/1‑150 (MQ=255)                                                                          aGATTAATGCCGGCAGATATTTATCGGAAAGGAGAAAATGCAGAGCTATTCGTCACTTTCAATGGCCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                          <  2:3803/150‑1 (MQ=255)                                                                          aGATTAAGGACGGAAGATATTTATGGGAAGGGAGAAAAGGCAGAGCTATTGTCCCCTTTAAAGGGCCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                          <  2:40770/150‑1 (MQ=255)                                                                             ttAATGACGGCAGATATTTATCGGAAAGGAGAAAATGAAGAGCTTTTCGTCACGTTAAAGGGTAAAGGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa                                                                        <  1:546189/149‑1 (MQ=255)                                                                               aaTGACGGAAGATATTTATAGGTAAGGAGAAAATGAAGAGTTTTTCGTCACTTTCATTGTTCAATCCCCCACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                      <  1:504953/149‑1 (MQ=255)                                                                                 tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                   >  1:792311/1‑150 (MQ=255)                                                                                   aCGGCATATATTTATGGGTAAGGAGGAAATGCAGACTTTTTGGTCCCGTTCAGTGGTCAATTCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                   <  1:183388/148‑1 (MQ=255)                                                                                    cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                               >  1:651460/1‑151 (MQ=255)                                                                                          atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTAGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa                                                          >  2:747138/1‑150 (MQ=255)                                                                                             tttagtggTAAGGAGAAAATGTAGACTTATTCGTCATGTTCAAGGGTAAACGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                        <  2:194026/143‑1 (MQ=255)                                                                                             tttATGGGTAAGGAGAAAATGCAGAGTTTTTCGTCACGTCAAGTGGTAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                        <  2:779297/149‑1 (MQ=255)                                                                                             tttATCGGTAAGGAGGAAATGCAGGGCTATTCGTCACGTTAAATGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                        <  2:456042/149‑1 (MQ=255)                                                                                              ttaggggaaaGGAGGAAAGGAAGAGTTTTTCGCCCCGTTAAATGGCAAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                      <  2:515659/142‑1 (MQ=255)                                                                                              ttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                      >  2:484409/1‑150 (MQ=255)                                                                                              ttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                       >  1:859079/1‑149 (MQ=255)                                                                                                  cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                 >  1:117160/1‑151 (MQ=255)                                                                                                   ggTAAGGAGAAAATGAAGATTTTTTGGTCGCTTTAAATGATAAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                      <  1:170621/145‑1 (MQ=255)                                                                                                     tAAGGAGAAAATGCAGAGCTATTCGTCGCGTTCAATGGACAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                               <  2:103026/150‑1 (MQ=255)                                                                                                      aaGGAGAAAATGAAGAGCTATTCTCCACTTTAAATGGTAAAGGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                              <  1:327382/150‑1 (MQ=255)                                                                                                         gagaaAAGGCAGAGTTATTCGTCACGTTCAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                          <  2:68334/151‑1 (MQ=255)                                                                                                               aTGCAGACCTTTTCGTCATTTTAAATGTTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                                    <  1:693517/151‑1 (MQ=255)                                                                                                                 gCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                  >  1:136400/1‑151 (MQ=255)                                                                                                                     agCTTTTCTTCACGTTCAATGGAAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                               <  2:170624/150‑1 (MQ=255)                                                                                                                        ttttCGTCCCTTTCAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                            <  1:306028/148‑1 (MQ=255)                                                                                                                           tCGTCATTTTCAAGGTTCACCGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                         <  1:883433/150‑1 (MQ=255)                                                                                                                            cGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                       >  2:701351/1‑151 (MQ=255)                                                                                                                              tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                     >  2:471466/1‑151 (MQ=255)                                                                                                                               cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                    >  1:771852/1‑151 (MQ=255)                                                                                                                                                  cacaCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg  >  2:798591/1‑150 (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: ‑