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

AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984293‑2984591                                                                                                                                                     |                                                                                                                                                       aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                      >  2:867723/1‑151 (MQ=255)   tatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt                                                                                                                                                    >  1:671940/1‑151 (MQ=255)    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                                   <  2:1113125/151‑1 (MQ=255)     tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                  >  1:1195743/1‑151 (MQ=255)     tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                  >  2:534000/1‑151 (MQ=255)      gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                  <  2:538875/150‑1 (MQ=255)           gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccg                                                                                                                                            >  1:871644/1‑151 (MQ=255)              aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg                                                                                                                                         >  2:856318/1‑151 (MQ=255)               aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGa                                                                                                                                        <  2:399674/151‑1 (MQ=255)                  gAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATAc                                                                                                                                     >  1:308369/1‑151 (MQ=255)                  gAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATAc                                                                                                                                     >  1:371462/1‑151 (MQ=255)                         gggggACCCCAAAAAATGTATGTCGCTTTGCCACCAATCATTATTAAACATATAAATGACGGCAAATATTTATCGGTAATGAGAAAATGCAGATCTATTCTTACCGTACAATGGTAAATGCCCCCCCCTGGCCCCCGCCGATACAAAGc                                                                                                                                <  1:168573/147‑1 (MQ=255)                           gggAGCCAAAAAATTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAAATATTTATCGTAAGGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGtt                                                                                                                             <  1:217004/150‑1 (MQ=255)                                   aaaaaaTTTATGTCGACTTCCCGCCCATCATTATTCAACAGATAAAGAGCGGCAGATTTTTATCGGTAAGGAGAAAATGCAGATCTATTCGTCACTTTAAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                    <  1:936066/149‑1 (MQ=255)                                      aaaGGTATGTCGACTTGCCAGCAATCTTTATTCAAAAGATAAATGACGGAAGATATTTATGGGTAAGGAGAAAATGCAGAGTTTTTCTCCACTTTAAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                                 <  2:394109/151‑1 (MQ=255)                                           tatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                            >  2:983302/1‑151 (MQ=255)                                                cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGAAGATATTTATCGGTAAGGAGAAAAGGCAGAGCTTTTCGTCACGTTAAATGGTCAAGGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                         <  2:466186/149‑1 (MQ=255)                                                      gccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                                  >  2:49309/1‑150 (MQ=255)                                                       ccagccaTTATTATTCAAAAAATTAATGGCGGCAAATTTTTTTCGGTAAGGAGAAAAAGCAAAACTTTTTGTCCCTTTAAATGGTAAACCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                    <  2:22058/147‑1 (MQ=255)                                                       ccagccaTCATTATTCAACAGATTAATGAGGGCAGATTTTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAATGGTAAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                                  <  2:1020345/149‑1 (MQ=255)                                                       ccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAc                                                                                                >  2:512921/1‑151 (MQ=255)                                                       ccagccaTCATTATCAAACAGATTAAGGACGGCAGATATTTATGGGAAAGGAGGAAAGGCGGGTCTTTTCGCAACGTTAAAGGGTCAACGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                                  <  2:794103/149‑1 (MQ=255)                                                       ccaccaatcatTATACAAAAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTAATTGTTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                                  <  1:511587/143‑1 (MQ=255)                                                           ccttctttatcCAACAGATTAAAGAGGGAAGATTTTTAGGGTGAAGGAGGAAATCCAGTGCTTTTCGTGTGTTTAAGGGGCCAACGCCCCGCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               <  2:764249/142‑1 (MQ=255)                                                           ccaTAATTATTAAACAGATTAATGACGGCAATTATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCCCTTTAAATGGTAAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                             <  2:863988/150‑1 (MQ=255)                                                           ccaTAATTATCAAACAGATTAATGCCGGCAGATATTTATCGGAAGGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                             <  2:982489/150‑1 (MQ=255)                                                                 ttattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                       >  2:622310/1‑150 (MQ=255)                                                                    ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                   >  2:443552/1‑151 (MQ=255)                                                                    tacaaaaaaTTAATGACGGAATATTTTTATGGGTAGGGAGAAAATGCAGAGTTATTCGTCACGTTAAATGGTCATTCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                    <  2:784762/142‑1 (MQ=255)                                                                       aacaattaaaTGCCGGCAGATATTTAGGGTTAGGGAGGAAATGCAGCGTTTTTCGGCTTTTTAAATGGTCAATCCCCCTCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                     <  1:602052/138‑1 (MQ=255)                                                                       aaCAGATAAATGACGGAAAATATTTATGGGTAAGGAAGAAAGGCAGAGCTTTTCGTCACGTTCAATGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                      <  1:433965/145‑1 (MQ=255)                                                                       aaCAAATAAAGGACGGAAAATATTTATGGGAAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAAGGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                                <  2:710631/151‑1 (MQ=255)                                                                          aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                              >  2:238342/1‑150 (MQ=255)                                                                           gATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                             >  2:287393/1‑150 (MQ=255)                                                                             ttAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCCTGAACCAGAGAAAcc                                                                          >  1:561764/1‑151 (MQ=255)                                                                               aaaggcGGCAGTTTTTTAGAGGGAAGGAGGAATTGCAGGGCTATTCGATTTGTTAAGTGGCGACCGCCACCCCCGGGCACCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg                                                                        <  1:346352/146‑1 (MQ=255)                                                                                aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                       >  1:213053/1‑151 (MQ=255)                                                                                aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                       >  2:987668/1‑151 (MQ=255)                                                                                       cAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                >  2:1044084/1‑151 (MQ=255)                                                                                       cAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                >  2:1026829/1‑151 (MQ=255)                                                                                        aaatttttATCGGTAAGAAGAAAATGCAAACCTTTTCGCCCCTTTCATTGGTCAAGCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                                 <  1:1113982/144‑1 (MQ=255)                                                                                        aGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                >  1:446841/1‑150 (MQ=255)                                                                                               tATGGGAAAGGAGGAAATGCAGATCTATTCGTCACTTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                         <  1:529385/150‑1 (MQ=255)                                                                                                  cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                     >  1:1029442/1‑151 (MQ=255)                                                                                                    gTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                                   >  2:1128335/1‑151 (MQ=255)                                                                                                      aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                  >  1:170163/1‑150 (MQ=255)                                                                                                             aaaGGCAGAGCTTTTCGTCACGTTCAATGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt                                           <  2:895449/150‑1 (MQ=255)                                                                                                                   agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  1:455747/1‑149 (MQ=255)                                                                                                                        ttttCGTCCCGTTAAAGGGACAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                                <  2:667158/148‑1 (MQ=255)                                                                                                                            cGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                           >  2:287561/1‑151 (MQ=255)                                                                                                                               ccctttCAATGGTCAACCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                            <  1:1063570/143‑1 (MQ=255)                                                                                                                                actttCAATGGTAAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             <  1:366456/142‑1 (MQ=255)                                                                                                                                  gTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     >  1:540206/1‑151 (MQ=255)                                                                                                                                      aaTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg                 >  1:862014/1‑151 (MQ=255)                                                                                                                                             aaTGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt            <  1:480744/149‑1 (MQ=255)                                                                                                                                                     acCCTGGCCCCCACCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  1:1160873/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                       AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984293‑2984591

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